Question: Plotting the first and second derivative of the numerical solution of an ODE

Hello, I am solving an ODE numerically:

restart: deq1 := diff(y(t), t, t)+y(t)*abs(y(t)) = 0; ic1 := y(0) = 1.0, (D(y))(0) = 0.; 
dsol1 := dsolve({deq1, ic1}, numeric, range = 0 .. 10); 
odeplot(dsol1);

I need to get the plots of the first and second derivatives of the solution. Ok, I can get the first derivating by transforming the equation to system with two unknown functions and reducing the order, but I also need the second and possibly higher derivatives. 

 

Please Wait...