Solve stiff differential equations and DAEs; variable order method - MATLAB ode15s

This MATLAB function with tspan = [t0 tf] integrates the system of differential equations y′ = f(t,y) from time t0 to tf with initial conditions y0. ... If you specify an output function as the value of the OutputFcn property, the solver calls it with the

www.mathworks.com