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 an effective local error control algorithm for RK5GL3, which uses local extrapolation with an eighth-order Runge-Kutta method in tandem with RK5GL3, and a Hermite interpolating polynomial for solution estimation at the Gauss-Legendre quadrature nodes. Keywords—RK5GL3, RKrGLm, Runge-Kutta, Gauss-Legendre, Hermite interpolating polynomial, initial value problem, local error.
[1]
J. Butcher.
Numerical methods for ordinary differential equations
,
2003
.
[2]
Ernst Hairer,et al.
Solving Ordinary Differential Equations I: Nonstiff Problems
,
2009
.
[3]
J. Dormand,et al.
A family of embedded Runge-Kutta formulae
,
1980
.
[4]
W. Cheney,et al.
Numerical analysis: mathematics of scientific computing (2nd ed)
,
1991
.
[5]
J. Douglas Faires,et al.
Numerical Analysis
,
1981
.
[6]
T. E. Hull,et al.
Comparing Numerical Methods for Ordinary Differential Equations
,
1972
.