> | restart:with(DEtools):with(plottools):with(plots): |
An animation of a vector field under going a birfucation
I will include a disiption soon.
> | PP:=display(seq(phaseportrait([D(x)(t)=x(t)*(y(t)^2+4),
D(y)(t)=(x(t)-y(t)^2+alpha/10)], [x(t),y(t)],t=0..1, [[x(0)=0,y(0)=0]], x=-5..5,y=-5..5,color=blue,linecolor=red),alpha=-30..30),insequence=true): |
> | parabola:=display(seq(plot([t^2-alpha/10,t,t=-5..5],x=-5..5,y=-5..5,color=green),alpha=-30..30),insequence=true): |
> | display(PP,parabola); |
> |