(1) Mathematisches Institut, Universität Tübingen, Auf der Morgenstelle 10, D-72076 Tübingen, Germany;(2) Institut für Mathematik II, BioComputing Group, Freie Universität Berlin, Arnimallee 2-6, D-14195 Berlin, Germany
Abstract:
Numerical integrators for second-order differential equations with time-dependent high frequencies are proposed and analysed. We derive two such methods, called the adiabatic midpoint rule and the adiabatic Magnus method. The integrators are based on a transformation of the problem to adiabatic variables and an expansion technique for the oscillatory integrals. They can be used with far larger step sizes than those required by traditional schemes, as is illustrated by numerical experiments. We prove second-order error bounds with step sizes significantly larger than the almost-period of the fastest oscillations.AMS subject classification (2000) 65L05, 65L70.Received February 2004. Accepted February 2005. Communicated by Syvert Nørsett.
