A model for stability of the semi-implicit backward differentiation formulas

摘要

A model is presented for stability for an extension of linear multistep methods for stiff ordinary differential equations. The method is based on a prediction followed by a fixed number of corrections obtained by a Newton scheme with inexact Jacobian matrix. The impact on stability of error in the matrix over a broad range of linear, constant coefficient equations is modeled. The model provides practical guidance for implementation of software for stiff equations.