Efficient split linear multistep methods for stiff ordinary differential equations

A new family of predictor-corrector schemes is designed for the numerical solution of stiff differential systems. Based on split Adams–Moulton formulas through sixth order, members of the new family achieve higher order and possess smaller error constants than corresponding split backward differentiation formulas of the same stepnumber, while maintaining similar stability properties. Some confirmation of this is obtained using a variable step implementation on test problems from the literature.