feny.mws

>    deq1:=diff(x(t),t,t)=g*(1+mu)*x(t)/L-mu*g;

deq1 := diff(x(t),`$`(t,2)) = 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;

g := 10

mu := .2

L := 1

d := .2

>    DEplot(deq1,x(t),t=0..1.5,
[[x(0)=d,D(x)(0)=0]],scene=[t,x(t)],linecolor=t,stepsize=0.01);

[Maple Plot]

>