QPPS

Solve systems of ordinary differential equations with polynomial right-hand sides.

We present a general purpose solver for non-stiff initial value problems in ordinary differential equations using the Parker-Sochacki method. The transformation of a system of ODEs into a system with polynomial right-hand sides is illustrated. Step size and order are selected automatically to satisfy the user's error tolerance. User-supplied reporting times, event functions, and output functions are accommodated.

For details, see qpps.pdf in the .zip file. (For some reason this Description panel provides no "attachment" link.)