ppanimation.mw

> 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);

[Plot]

>