Numerical Solution of Stiff Ordinary Differential Equations.

An analysis is presented for alternate numerical techniques for solving stiff ordinary differential equations. These techniques include a singular perturbation or pseudo-steady-state method and an imbedded, error-monitoring semi-implicit Runge-Kutta method. Extensive numerical experience on equations which are linearnonlinear, smalllarge dimensional, and moderatelystrongly stiff reveals that the singular perturbation method is most efficient for very stiff problems while the imbedded Runge-Kutta method is superb over a wide range of stiffness. Author