RKF45 - Runge-Kutta-Fehlberg ODE Solver

RKF45 is the Watt and Shampine implementation of the Runge-Kutta-Fehlberg ordinary differential equation solver, with automatic error estimation. Both single and double precision versions are available.

Files you may copy include:

• rkf45.f90, the source code;
• rkf45_prb.f90, a sample calling program;
• rkf45_prb.out, the output from the sample calling program;

The list of routines includes:

• FEHL takes one Fehlberg fourth-fifth order step (single precision).
• FEHL_D takes one Fehlberg fourth-fifth order step (double precision).
• RKF45 carries out the Runge-Kutta-Fehlberg method (single precision).
• RKF45_D carries out the Runge-Kutta-Fehlberg method (double precision).
• RKFS implements the Runge-Kutta-Fehlberg method (single precision).
• RKFS_D implements the Runge-Kutta-Fehlberg method (double precision).
• TIMESTAMP prints the current YMDHMS date as a time stamp.

Return to the FORTRAN software page.