Differential equations

Runge-Kutta-Cash-Karp method
Adaptive stepsize, moderate accuracy