Analysis of error propagation in the RK3GL2 method

J. S. C. Prentice
2024-08-14
Abstract:The RK3GL2 method is a numerical method for solving initial value problems in ordinary differential equations, and is a hybrid of a third-order Runge-Kutta method and two-point Gauss-Legendre quadrature. In this paper we present an analytical study of the propagation of local errors in this method, and show that the global order of RK3GL2 is expected to be four.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to analyze the propagation of local errors in the RK3GL2 method and show why its global error can reach fourth - order accuracy. ### Problem Background The RK3GL2 method is a numerical method for solving initial - value problems of ordinary differential equations. It combines the third - order Runge - Kutta (RK) method and the two - point Gauss - Legendre (GL) quadrature method. Generally, the local error of the Runge - Kutta method is \(O(h^{r + 1})\), and the global error is \(O(h^r)\), where \(h\) is the step size. However, by introducing the GL quadrature method, the RK3GL2 method can increase the order of the global error under certain conditions. ### Core Content of the Paper The main contributions of the paper are as follows: 1. **Analysis of Local Error Propagation**: The author analyzes in detail how local errors in the RK3GL2 method propagate and shows the cumulative effects of these errors at different nodes. 2. **Improvement of Global Error Order**: By introducing the GL quadrature method, the author proves that the global error of the RK3GL2 method can reach fourth - order accuracy, that is, \(O(h^4)\). This is an improvement over the global error \(O(h^3)\) of the traditional third - order Runge - Kutta method. ### Key Formulas - **Definition of Local Error**: \[ \varepsilon_{i + 1}=\left[y_i+h F(x_i, y_i)\right]-y_{i + 1} \] where \(y_i\) is the exact solution, \(w_i\) is the numerical solution, and \(F(x, y)\) is a function related to a specific RK method. - **Definition of Global Error**: \[ \Delta_i = w_i - y_i \] - **Error Propagation Formulas**: \[ \Delta_3=\varepsilon_3+A_{1,2}h \] \[ \Delta_6 = (\varepsilon_6+\varepsilon_3)+(A_{4,5}h+A_{1,2}h)+B_6\Delta_3h \] - **Summary of Global Error**: \[ \Delta_{3N}=O(h^4) \] ### Conclusion Through the above analysis, the author shows the propagation mechanism of local errors in the RK3GL2 method and proves that its global error can reach fourth - order accuracy. This result is of great significance for improving methods of numerically solving ordinary differential equations, especially in cases where high - precision solutions are required. ### References The paper cites multiple references on classical numerical analysis and solving ordinary differential equations, such as "Solving Ordinary Differential Equations I: Nonstiff Problems" by E. Hairer et al. and "Numerical Methods for Ordinary Differential Equations" by J.C. Butcher, ensuring the theoretical basis of the research and the reliability of the method.