SciPy needs an adaptive 4th order Runge-Kutta method...