> | deq1:=diff(x(t),t,t)=g*(1+mu)*x(t)/L-mu*g; |
> | with(DEtools): |
> | g:=10; mu:=0.2; L:=1; d:=L*0.2; if d*(1+mu)-mu*L>0 then "parameters ok" else "parameters not ok!!!" fi; |
> | DEplot(deq1,x(t),t=0..1.5, [[x(0)=d,D(x)(0)=0]],scene=[t,x(t)],linecolor=t,stepsize=0.01); |
> |