nm

5404 Reputation

17 Badges

8 years, 229 days

MaplePrimes Activity


These are answers submitted by nm

one way might be to use solve or PDEtools:-Solve

restart;
eq:=5*x^3-x^2+x-1=0;
PDEtools:-Solve(eq,x)

second order ODE can be written as two first order ode's. So convert your second order ODE by introducting two state variables, x1(t) and x2(t). One for "position" and second for "speed".  

Since you did not give an example of your second order ode's and its IC, I made one up.

restart;
ode := diff(x(t),t$2) +1/2*diff(x(t),t)+ x(t) = 0;
eq1:=diff(x1(t),t)=x2(t);
eq2:=subs([diff(x(t),t)=x2(t),x(t)=x1(t)],ode);
DEtools:-DEplot([eq1,eq2],[x1(t),x2(t)],t=0..35,[[x1(0)=1,x2(0)=1]],x1=-2..2,x2=-2..2,
        numpoints=200, linecolor=black, axes=boxed);