ODE
Shampine and Gordon ODE Solver

ODE is a FORTRAN90 library which implements the Shampine and Gordon ODE solver.

Given a system of ordinary differential equations of the form

        Y' = F(T,Y)
        Y(T0) = Y0
      

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.

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 implements the Runge-Kutta-Fehlberg ODE solver.

TEST_ODE is a FORTRAN90 library which defines test problems for ODE solvers.

Reference:

1. Lawrence Shampine, Marilyn Gordon,
   Computer Solution of Ordinary Differential Equations: The Initial Value Problem,
   Freeman, 1975,
   ISBN: 0716704617,
   LC: QA372.S416.

Source Code:

• ode.f90, the source code.
• ode.csh, commands to compile the source code.

Examples and Tests:

• ode_prb.f90, a sample calling program.
• ode_prb.csh, commands to compile and run the sample program.
• ode_prb_output.txt, the output file.

List of Routines:

• ODE is the user interface to an ordinary differential equation solver.
• DE carries out the ODE solution algorithm.
• STEP integrates the system of ODE's one step, from X to X+H.
• INTRP approximates the solution at XOUT by polynomial interpolation.
• TIMESTAMP prints the current YMDHMS date as a time stamp.

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