Second derivative of high-order accuracy methods for the numerical integration of stiff initial value problems

Abstract

In Mitsui (Sci Eng Rev Doshisha Univ Jpn 51(3):181–190, 2010) the author introduced a new class of discrete variable methods known as look-ahead linear multistep methods (LALMMs) which consist of a pair of predictor-corrector (PC), including a function value at one more step beyond the present step for the numerical solution of ordinary differential equations. Two-step family of Look-Ahead linear multistep methods of fourth-order pair were derived and shown to be A( $$\theta $$ )-stable in Mitsui and Yakubu (Sci Eng Rev Doshisha Univ Jpn 52(3):181–188, 2011). The derived integration methods are of low orders and unfortunately cannot cope with stiff systems of ordinary differential equations. In this paper we extend the concept adopted in Mitsui and Yakubu (Sci Eng Rev Doshisha Univ Jpn 52(3):181–188, 2011) to construct second-derivative of high-order accuracy methods with off-step points which behave essentially like one-step methods. The resulting integration methods are A-stable, convergent, with large regions of absolute stability, suitable for stiff systems of ordinary differential equations. Numerical comparisons of the new methods have been made and enormous gains in efficiency are achieved.