Diagonally implicit Runge-Kutta-Nyström methods for oscillatory problems. (English) Zbl 0676.65072

A number of second order two and three stage diagonally implicit Runge- Kutta-Nyström (DIRKN) methods with low phase error are developed for the initial value problem \((i)\quad y''=f(t,y),\) \((ii)\quad y(0)=y_ 0,\) \(y'(0)=y'_ 0.\) The definition of phase error and amplitude error for the Runge-Kutta-Nyström algorithm associated with (i), (ii) has been developed by the authors in an earlier paper [ibid. 24, 595-617 (1987; Zbl 0624.65058)]. The phase error \(\phi\) (\(\omega\) h) and amplification error \(\alpha\) (\(\omega\) h) are defined as extensions of the explicitly known formulas for deviation of the phase and amplitude of the RKN solution of (i), (ii) from the true solution when \(f=-\omega^ 2y\). An RKN method is said to be dispersive and dissipative of order q and r if \((iii)\quad \phi (\omega h)=O(h^{q+1}),\) \((iv)\quad \alpha (\omega h)=O(h^{r+1}).\) The authors construct two-stage second order DIRKN algorithms in which \(q=4,6\), one of which is unconditionally strongly stable. Also three-stage algorithms are constructed with \(q=6,8,10\), one of which is P-stable. Results of computation of actual dispersion error are presented using the new algorithms on five test problems. The new algorithms are better than standard DIRKN methods in case the dispersion error grows larger than the truncation error.
Reviewer: J.B.Butler jun.


65L05 Numerical methods for initial value problems involving ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
34C10 Oscillation theory, zeros, disconjugacy and comparison theory for ordinary differential equations


Zbl 0624.65058
Full Text: DOI