Added a fourth order runge-kutta solver.