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

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