cheers,
Sanjay.
What's wrong with the built-in function NDSolve? Although it is a
general purpose integrator which uses the Adams method or the Gear
method depending on stiffness, you can choose to use a RKF method if
you want, with the option
Method -> RungeKutta
which uses a Fehlberg order 4-5 Runge-Kutta algorithm.
Carl Woll
Dept of Physics
U of Washington