An attempt to avoid exact Jacobian and nonlinear equations in the numerical solution of stiff differential equations

Abstract

A class of linear implicit methods for numerical solution of stiff ODE’s is presented. These require <italic>only</italic> occasional calculation of the Jacobian matrix while maintaining stability. Especially, an effective second order stable algorithm with automatic stepsize control is designed and tested.