Error Propagation in the RK5GL3 Method

Authors: J.S.C. Prentice

Abstract:

The RK5GL3 method is a numerical method for solving initial value problems in ordinary differential equations, and is based on a combination of a fifth-order Runge-Kutta method and 3-point Gauss-Legendre quadrature. In this paper we describe the propagation of local errors in this method, and show that the global order of RK5GL3 is expected to be six, one better than the underlying Runge- Kutta method.

Keywords: RK5GL3, RKrGLm, Runge-Kutta, Gauss-Legendre, initial value problem, order, local error, global error.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1076954

References:

[1] J.S.C. Prentice, "The RKGL method for the numerical solution of initialvalue problems", Journal of Computational and Applied Mathematics, 213, 2 (2008) 477.
[2] E. Hairer, S.P. Norsett, and G. Wanner, Solving ordinary differential equations I: Nonstiff problems, Berlin: Springer-Verlag, 2000, p177.
[3] J.C. Butcher, Numerical methods for ordinary differential equations, Chichester: Wiley, 2003, pp151 − 155.
[4] D. Kincaid andW. Cheney, Numerical Analysis: Mathematics of Scientific Computing, 3rd ed., Pacific Grove: Brooks/Cole, 2002, pp492 − 498.
[5] T.E. Hull, W.H. Enright, B.M Fellen, and A.E. Sedgwick, "Comparing numerical methods for ordinary differential equations", SIAM Journal of Numerical Analysis, 9, 4 (1972) 603.