hidro.mws

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

deq1 := diff(x(t),`$`(t,2)) = C/m/x(t)^2/y(t)

deq2 := diff(y(t),`$`(t,2)) = C/m/x(t)/y(t)^2

varlist := {x(t), y(t)}

>    m:=1;C:=20;

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

[Maple Plot]

>