Midpoint, Heun and Runge–Kutta Methods — Cheat sheet
On this page
Core rule
Watch for
Higher order does not remove stability limits or repair discontinuous forcing automatically. Split a step at known discontinuities. Adaptive methods estimate local error and adjust steps, but the requested tolerance is not an unconditional bound on every global error.