Is there any documentation some where that talks about the algorithm used to control the step size of the Matlab Ode solvers? How do they calculate the truncation error how do they change the step size... Etc
Thanks
I am not an expert on the topic, but a safe assumption is: it depends. AFAIK there are Runge-Kutta (ode23, ode45) and Adams-Bashforth-Moulton solvers in use. Information about these algorithms can e.g. be found in the references given in http://www.mathworks.de/de/help/matlab/ref/ode45.html
Good luck!
Most of the information regarding specific solvers is under the "Algorithms" section indicated by Marc:
http://www.mathworks.co.uk/help/matlab/ref/ode45.html#bti6n8p-40
For instance, ode45 uses a Dormand-Prince pair, ode113 is a variable order Adams-Bashforth-Moulton PECE solver, ode15s is based on numerical differentiation formulas, etc.
Maybe this will be of interest to you:
http://www.cs.berkeley.edu/~wkahan/Math128/ODEsuite.pdf
In any case, ode functions have open codes, so one can find the most precise information there.
- Badges
- Science topic