Higher order implicit multistep methods for matrix differential equations

This paper proposes implicit multistep matrix methods for the numerical solution of stiff initial value matrix problems. The study of matrix difference equations involving the matrix coefficients of the multistep method permits one to obtain convergence results, as well as bounds for the global discretization error in terms of the data. An illustrative example is included.