Runge–Kutta integration