> | with(DEtools): |
> | eq1:=diff(y(x),x)=1/5*y(x)*(1-y(x)/10); |
> | DEplot(eq1,y(x),x=0..10,[[y(0)=0],[y(0)=2],[y(0)=5],[y(0)=10],[y(0)=15]],y=0..15,title="Direction Field",color=black,dirgrid=[30,30],linecolor=blue); |
> | eq2:=diff(y(x),x)=1/5*y(x)*(y(x)-2)*(1-y(x)/10); |
> | DEplot(eq2,y(x),x=0..10,[[y(0)=0],[y(0)=1],[y(0)=2],[y(0)=4],[y(0)=10],[y(0)=15]],y=0..15,title="Direction Field",color=black,dirgrid=[30,30],linecolor=blue); |
> | eq3:=diff(p(t),t)=p(t)*(1-p(t))-2*p(t); |
> | DEplot(eq3,p(t),t=0..10,[[p(0)=2],[p(0)=4],[p(0)=15]],p=0..15,title="Direction Field",color=black,dirgrid=[30,30],linecolor=blue); |
> | eq4:=diff(p(t),t)=2*p(t)*(1-p(t))-p(t); |
> | DEplot(eq4,p(t),t=0..10,[[p(0)=0.1],[p(0)=0.5],[p(0)=4]],p=0..5,title="Direction Field",color=black,dirgrid=[30,30],linecolor=blue); |
> |