How to take the derivative for the retained value from ode45 solver?

[t,y]= ode45(@my_function,[0,10],[5,0]);

I'm getting angular velocity and angular position from the ODE solver above. Now I need to find the angular acceleration, I'm suppose to take the derivative, but I'm not sure how to do that?

[Jan]

Please do not post a question twice. You can and should delete this duplicate question as long as nobody had posted an answer.