I can't understand how the deval function works in matlab, it's hard for me to understand the code, are there any books or descriptions of numerical methods where I can read about numerical algorithms, how to evaluate solutions of differential equations at the right points? I need to implement this feature myself. Thanks for the help!
$\endgroup$ 11 Answer
$\begingroup$The solvers in last mode mentioned in the documentation of them produce as part of the returned object a "dense output" solution, that is, a piecewise polynomial interpolation, represented by the coefficients for every segment. Many of the standard solvers can produce this interpolation directly from the values of the ODE function produced during the solver step, or with a few additional evaluations.
So deval has just to determine the correct segment, the offset inside that segment, and evaluate the polynomial.
