> | deq1:=diff(x(t),t,t)=C/m/(x(t)^2*y(t)); deq2:=diff(y(t),t,t)=C/m/(x(t)*y(t)^2); varlist:={y(t),x(t)}; |
> | m:=1;C:=20; |
> | with(DEtools): |
> | DEplot({deq1,deq2},{x(t),y(t)},t=0..2, [[x(0)=10,y(0)=8,D(x)(0)=-10,D(y)(0)=-10]],scene=[x(t),y(t)],linecolor=t,stepsize=0.01); |
> |