I'm solving an ODE using ODE15i. As per the usual use case of 15i, the function I'm integrating is a function of t,y,and prime (the time derivative of y):
function sol = myFunc
sol = ode15i(@nested_integrand,tspan,y0,yp0)
function f = nested_integrand(t,y,yp)
However, the solution sol only contains values of t and y at the integration points. Not yp! This seems like an oversight to me, since the consistent values of yp have already been calculated within ode15i.
So my question is in two parts:
1. What is the most elegant way of getting the yp values that correspond with sol.t and sol.y?
2. Does anyone else have a use for this output? If so I'll make a feature request.