CWG_ODE
Three ODE Solvers by William Gear.

CWG_ODE is a FORTRAN77 library, using single precision arithmetic, which includes three solvers for systems of ordinary differential equations, written by C William Gear.

Related Data and Programs:

LSODI is a FORTRAN77 library which implements the Livermore Solver for Ordinary Differential Equations, Implicit case.

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

ODE is a FORTRAN90 library which contain the Shampine and Allen ODE solver.

POLKING is a pair of MATLAB routines, DFIELD and PPLANE, for solving and displaying the solution of differential equations and their phase planes.

RKF45 is a FORTRAN90 library which applies the Runge-Kutta-Fehlberg method to solve a system of ODE's.

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

Reference:

Charles Gear,
Numerical Initial Value Problems in Ordinary Differential Equations,
Prentice-Hall, 1971,
ISBN: 0136266061,
LC: QA372.G4.

Source Code:

cwg_ode.f, the source code.
cwg_ode.csh, commands to compile the source code.

Examples and Tests:

cwg_ode_prb.f, a sample calling program.
cwg_ode_prb.csh, commands to compile and run the sample program.
cwg_ode_prb_output.txt, the output file.

List of Routines:

SRKODE applies a Runge-Kutta method to differential equations.
SRK1 performs one Runge Kutta step.
SPRODE uses polynomial or rational approximation to solve ODE's.
SMVODE solves an ODE with a multi-value method.
PDUMMY is a dummy derivative routine for calls to SMVODE with mf.ne.1.
SGEFA factors a matrix by gaussian elimination.
SGESL solves a linear system that was factored by SGEFA.
ISAMAX finds the index of element having maximum absolute value.
SAXPY computes constant times a vector plus a vector.
SDOT forms the dot product of two vectors.
SSCAL scales a vector by a constant.
TIMESTAMP prints out the current YMDHMS date as a timestamp.

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