I've adapted what I thought were the best ideas for an ordinary
differential equation integrator and placed them on the internet at:
http://forums.devshed.com/programming-languages-139/numerical-recipes-odeint-961791.html
I'd like to see the code installed into the j addons/math/misc
distribution. I believe I have permission without the "know how".
The code is adaptive, taken from the Patrick Harrington's Numerical
Recipes translation;
uses a gerund for the derivative verbs from the EelVex code;
and includes a generalized version of my stepper for which you supply
the Butcher tableau.
The file includes the common 3 RK45 steppers.
Thank you all!
(And since I haven't got Dr. Iverson's cleverness, I'll let the compiler
allocate registers and delay lines and relays. Thank you for showing
his doctoral thesis.)
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm