An ecient family of second derivative Runge-Kutta collocation methods for oscillatory systems
Abstract
An ecient family of high-order second derivative Runge-Kutta collocation methods is de-
rived for the numerical solution of oscillatory systems. The approach uses polynomial inter-
polation and collocation techniques to construct continuous schemes which were evaluated
at both step and o-step points to obtain hybrid formulae. The hybrid formulae can be
applied simultaneously as block methods for moving the integration process forward at a
time, if desired. The block methods based on hybrid formulation can also be converted to
second derivative Runge-Kutta collocation methods. The stability properties and order of
accuracy of the methods are studied. They can also be implemented easily since they are
collocation methods and provide a high order of accuracy. The methods were illustrated by
the applications to some test problems of oscillatory system found in the literature and the
numerical results obtained conrm the accuracy and eciency of the methods.
References
K.Burrage, J.C.Butcher, Nonlinear stability for a general class of dierential equation methods.
BIT Numer. Math. 20(1979) 185 {203
S.O. Fatunla, Numerical integrators for sti and highly oscillatory dierential equations. J.
Math. Comp. 34 (150) (1980)373{390.
L.R. Petzold, L.O. Jay, J.Yen, Numerical solution of highly oscillatory ordinary dierential
equations. Acta Numerica. 6 (1997)437{483.
F. Costabile, A. Napoli, Stability of Chebychev collocation methods, Comput. Math. Appl. 47
(2004) 659-666.
F. Costabile, A. Napoli, A class of collocation methods for numerical integration of initial
value problems, Comput. Math. Appl. 62 (2011)3221{3235.
M. Urabe, An implicit one-step method of high-order accuracy for the numerical integration
of ordinary dierential equations. J. Numer. Math. 15(1970) 151 { 164.
H. Shantani, On one-step methods utilizing the second derivative. Hiroshima Math. 1 (1971)
{ 372.
H. Shantani, On explicit one-step methods utilizing the second derivative. Hiroshima Math. 2
(1972) 353 { 368.
J.R.Cash, High order methods for numerical integration of ordinary dierential equations.
Numer. Math. 30 (1978) 385 { 409.
T. Mitsui, Runge-Kutta type integration formulas including the evaluation of the second
derivative Part I Publ. RIMS Kyoto University ,Japan 18(1982) 325 { 364.
S.D. Capper, J.R. Cash, D. R. Moore, Lobatto-Obrechko formulae for 2nd , order two-points
boundary value problems. J. Numer. Anal. Indus. and Appl. Math. JNAIAM 314(2006) 13 {
R.P.K. Chan, A.Y.J. Tsai, On explicit two-derivative Runge-Kutta methods. J. Numer. Algor.
(2010) 171 { 194.
D.G. Yakubu, S. Markus, Second derivative of high order accuracy methods for the numerical
integration of sti initial value problems. Afrika Math. 27 (5-6)(2016) 963{977.
J.C. Butcher, Numerical Methods for Ordinary Dierential Equations. 2nd edn. Wiley Chichester,
T.Mitsui, D.G. Yakubu, Two-step family of look-ahead linear multistep methods for ODEs.
The science and engineering review Doshisha University Japan. 52(3)(2011) 181 { 188.
P. Onumanyi, D.O. Awoyemi, S.N. Jator, U.W. Sirisena, New linear multistep methods with
continuous coecients for rst order initial value problems. J. Nig. Math. Soci. 13(1994) 37{51.
I. Lie, S. P. Nrsett, Superconvergence for multistep collocation. Math. Comp. 52 (1989) 65{
D.G. Yakubu, M. Aminu, A. Aminu, The numerical integration of sti systems using stable
multi-step multi-derivative methods, J. Mod. Meth. Numer. Math., 8(1-2)(2017) 99{117.
