The
Cash-Karp formula is a 5(4) Runge-Kutta formula for
solving non-stiff systems of initial value ordinary differential equations. It
is fully described in the book

*Numerical
Recipes in C (and also Numerical Recipes in Fortran):
The Art of Scientific Computing (second edition)*

*by** Press, Flannery, Teukolsky and Vetterling.*

It
is also described in

*A
variable order Runge-Kutta method for initial value
problems with rapidly varying right hand sides*

*J.R.
Cash and A.H. Karp, ACM Transactions on Mathematical Software, Vol 16, No.3, Pages 201-222*

The
MATLAB code can be found in

*Numerical
Methods in Engineering with Matlab*

*Jaan
Kiusalaas, CUP, August 2005, pages 280-285*

Programs
which implement this method can be found in the ODE section of the above
Numerical Recipes books.

Then
starting from compute and define

A
solution at an off step point can be computed at
virtually zero extra cost using

,

,

,

,

,

,

,

5947 people have visited this page since May 2008.