Runge-Kutta Integration (Algorithm 9 [D2])
CACM April, 1966
Thacher Jr., H. C.
