New second derivative multistep methods for stiff systems

In this paper we present details of a new class of second derivative multistep methods. Stability analysis is discussed and an improvement in stability region is obtained. With a simple modification we take advantage of calling for the solution of algebraic equations with the same coefficient matrix in each step. Moreover, using IOM to solve the resulting linear systems, the coefficient matrix is not needed explicitly. Some numerical experiments and comparison with several popular codes are given, showing strong superiority of this new class of methods. (c) 2005 Elsevier Inc. All rights reserved.