Integration of ODE


Adaptive Stepsize for R-K

Modified Midpoint Method

\[\begin{split}z_0 &= y(x) \\ z_1 &= z_0 + h f(x,z_0) \\ z_{m+1} &= z_{m-1} + 2h f(x+mh,z_m) \\ y(x+H) &\approx y_n = \frac{1}{2} \left( z_n + z_{n-1} + h f(x+H,z_n) \right) .\end{split}\]

This method contains only the even powers of \(h\) thus we can gain two orders of precision at a time by calculating one more correction.

Refs & Notes

Back to top

© 2016-2018, Lei Ma | Created with Sphinx and . | On GitHub | Physics Notebook Statistical Mechanics Notebook Neutrino Physics Notes Intelligence | Index | Page Source