One-leg variable-coefficient formulas for ordinary differential equations and local–global step size control

In this paper we discuss a class of numerical algorithms termed one-leg methods. This concept was introduced by Dahlquist in 1975 with the purpose of studying nonlinear stability properties of multistep methods for ordinary differential equations. Later, it was found out that these methods are themselves suitable for numerical integration because of good stability. Here, we investigate one-leg formulas on nonuniform grids. We prove that there exist zero-stable one-leg variable-coefficient methods at least up to order 11 and give examples of two-step methods of orders 2 and 3. In this paper we also develop local and global error estimation techniques for one-leg methods and implement them with the local–global step size selection suggested by Kulikov and Shindin in 1999. The goal of this error control is to obtain automatically numerical solutions for any reasonable accuracy set by the user. We show that the error control is more complicated in one-leg methods, especially when applied to stiff problems. Thus, we adapt our local–global step size selection strategy to one-leg methods.

[1]  Werner Liniger,et al.  Stability of Two-Step Methods for Variable Integration Steps , 1983 .

[2]  Werner Liniger,et al.  Contractive methods for stiff differential equations part I , 1978 .

[3]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[4]  G. Dahlquist.  Error analysis for a class of methods for stiff non-linear initial value problems , 1976 .

[5]  Rolf Dieter Grigorieff,et al.  Stability of multistep-methods on variable grids , 1983 .

[6]  W. Hundsdorfer,et al.  Convergence of linear multistep and one-leg methods for stiff nonlinear initial value problems , 1991 .

[7]  G. Dahlquist.  On One-Leg Multistep Methods , 1983 .

[8]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[9]  J. Butcher.  Numerical Methods for Ordinary Differential Equations: Butcher/Numerical Methods , 2005 .

[10]  Germund Dahlquist,et al.  G-stability is equivalent toA-stability , 1978 .

[11]  E. Hairer,et al.  Stiff and differential-algebraic problems , 1991 .

[12]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[13]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[14]  G. Yu. Kulikov,et al.  A technique for controlling the global error in multistep methods , 2000 .

[15]  M. Crouzeix,et al.  The Convergence of Variable-Stepsize, Variable-Formula, Multistep Methods , 1984 .

[16]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[17]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[18]  G. Yu. Kulikov,et al.  Global Error Estimation and Extrapolated Multistep Methods For Index 1 Differential-Algebraic Systems , 2005 .

[19]  G. Kulikov,et al.  One-leg Integration of Ordinary Differential Equations with Global Error Control , 2005 .