From: Yihan on
I am using the ode solvers.

I want to extract the derivatives from the ode solver too.

For example:

function dy = fun(t,y)
dy = zeros(2,1);
dy(1) = y(1)*y(2);
dy(2) = -0.5*y(1)*y(2);
end

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

I obtain 'y' by solving this.
However, I require 'dy' at the various 't' too.
Anyone know how to extract this?