RK4
Runge-Kutta 4th Order ODE Solver

RK4 is a FORTRAN90 library which implements a simple Runge-Kutta solver for an initial value problem.

The rk4() routine does not include any error estimator. It takes a single step at a time, requiring the user to control the step size.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

RK4 is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version and a Python version.

Related Data and Programs:

NMS, a FORTRAN90 library which includes the DDRIV package of ODE solvers.

ODE, a FORTRAN90 library which solves a system of ordinary differential equations, by Shampine and Gordon.

RKF45, a FORTRAN90 library which implements the Runge-Kutta-Fehlberg ODE solver.

TEST_ODE, a FORTRAN90 library which contains routines which define some test problems for ODE solvers.

Reference:

1. Erwin Fehlberg,
   Low-order Classical Runge-Kutta Formulas with Stepsize Control,
   NASA Technical Report R-315, 1969.
2. Lawrence Shampine, Herman Watts, S Davenport,
   Solving Non-stiff Ordinary Differential Equations - The State of the Art,
   SIAM Review,
   Volume 18, pages 376-411, 1976.

Source Code:

• rk4.f90, the source code.

Examples and Tests:

• rk4_test.f90, a sample calling program.
• rk4_test.txt, the output file.
• rk4_test.png, an image of the solution.

List of Routines:

• RK4 takes one Runge-Kutta step for a scalar ODE.
• RK4VEC takes one Runge-Kutta step for a vector ODE.

You can go up one level to the FORTRAN90 source codes.