Runge-Kutta-Fehlberg ODE Solver

RKF45 is a FORTRAN77 library which implements the Watt and Shampine RKF45 ODE solver.

The RKF45 ODE solver is a Runge-Kutta-Fehlberg algorithm for solving an ordinary differential equation, with automatic error estimation using rules of order 4 and 5.


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


RKF45 is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version and a PYTHON version.

Related Data and Programs:

CWG_ODE, a FORTRAN77 library which implements three ODE solvers by C William Gear.

DRIV, a FORTRAN77 library which solves real or complex systems of ordinary differential equations;

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

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

ODEPACK, a FORTRAN77 library which contains nine ODE solvers, including LSODE, LSODES, LSODA, LSODAR, LSODPK, LSODKR, LSODI, LSOIBT, and LSODIS, by Alan Hindmarsh.

RK4, a FORTRAN77 library which applies the fourth order Runge-Kutta algorithm to estimate the solution of an ordinary differential equation at the next time step.

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


  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.
  3. The source code for Shampine and Watt's original FORTRAN77 routine is available at http://www.netlib.org/ode/ the NETLIB ODE web site.

Source Code:

Examples and Tests:

List of Routines:

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

Last revised on 03 March 2010.