Fortran Runge-Kutta numerical solver