Trigonometrically-fitted Scheifele two-step methods for perturbed oscillators

In this paper, a new family of trigonometrically-fitted Scheifele two-step (TFSTS) methods for the numerical integration of perturbed oscillators is proposed and investigated. An essential feature of TFSTS methods is that they are exact in both the internal stages and the updates when solving the unperturbed harmonic oscillator y″=−ω2y for known frequency ω. Based on the linear operator theory, the necessary and sufficient conditions for TFSTS methods of up to order five are derived. Two specific TFSTS methods of orders four and five respectively are constructed and their stability and phase properties are examined. In the five numerical experiments carried out the new integrators are shown to be more efficient and competent than some well-known methods in the literature.