For all of the routines below, a user can very easily specify the differential equation to be integrated as a octave function and the actual integration is performed by the Netlib FORTRAN code.

As these routines stand they are a bit primative in the way that they input auxillary variables like absolute and relative tolerances. Changing the code to input these parameters is not difficult and was not done because at some point in the future it might be nice to have a common interface for all the routines.

*ode.f**by:*Shampine and Gordon*for:*ordinary differential equation initial-value problem solver*alg:*Adam's methods*prec:*double*rel:*excellent*age:*old*gams:*I1a1b*octave gateway:*ODE

*rkf45.f**by:*Watts and Shampine*for:*ordinary differential equation initial-value problem solver*alg:*Runge-Kutta Fehlberg fourth-fifth order*prec:*double*rel:*good*age:*old*gams:*I1a1a*octave gateway:*RKF45