High order adaptive methods of Nyström-Cowell type (Q1362354)

A class of unconditionally stable multistep methods for second-order periodic or quasiperiodic initial value problems is developed. Problems of this type occur very frequently in celestial mechanics, nonlinear oscillations and in other situations. The methods belong to the category of methods with nonconstant coefficients. The coefficients of the methods depend upon the ``frequency'' of the problem, i.e. upon a parameter \(\omega > 0\). The methods integrate exactly trigonometric functions and algebraic polynomials. It is known that these methods for \(\omega = 0\) are reduced to classical methods, i.e. methods with constant coefficients. In this paper adaptive Nyström-Cowell methods of arbitrarily high order of accuracy are constructed. These methods are reduced to classical Nyström-Cowell methods for \(\omega = 0\). Numerical illustrations indicate the efficiency of the new methods.