Initial Runge-Kutta method for ODE, plus interpolation.