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 ...

Sugiura, H. and T. Torii, A method for constructing generalized Runge-Kutta methods, Journal of Computational and Applied Mathematics 38 (1991) 399-410. In the implementation of an implicit Runge-Kutta formula, we need to solve systems of nonlinear equations. In this paper, we analyze the Newton iteration process and a modified Newton iteration process for solving these equations. Then we propo...

A set of validated numerical integration methods based on explicit and implicit Runge-Kutta schemes is presented to solve, in a guaranteed way, initial value problems of ordinary differential equations. Runge-Kutta methods are well-known to have strong stability properties, which make them appealing to be the basis of validated numerical integration methods. A new approach to bound the local tr...

In this paper a new embedded Singly Diagonally Implicit Runge-Kutta Nystrom fourth order in fifth order method for solving special second order initial value problems is derived. A standard set of test problems are tested upon and comparisons on the numerical results are made when the same set of test problems are reduced to first order systems and solved using the existing embedded diagonally ...

In this paper, we study the preservation of quadratic conservation laws of Runge-Kutta methods and partitioned Runge-Kutta methods for Hamiltonian PDEs and establish the relation between multi-symplecticity of Runge-Kutta method and its quadratic conservation laws. For Schrödinger equations and Dirac equations, the relation implies that multi-sympletic RungeKutta methods applied to equations wi...

System theory for numerical analysis has recently become a focus of research. In this paper we regard dynamics of Newton’s method as a nonlinear feedback system and derive convergence conditions, based on the internal model principle and systems of Lur’e type. We then focus our attention on the analysis of the region of absolute stability of Runge-Kutta type methods. We derive a linear matrix i...

8. First-Order Equations: Numerical Methods 8.1. Numerical Approximations 2 8.2. Explicit and Implicit Euler Methods 3 8.3. Explicit One-Step Methods Based on Taylor Approximation 4 8.3.1. Explicit Euler Method Revisited 4 8.3.2. Local and Global Errors 4 8.3.3. Higher-Order Taylor-Based Methods (not covered) 5 8.4. Explicit One-Step Methods Based on Quadrature 6 8.4.1. Explicit Euler Method Re...

Based on reasonable testing model problems, we study the preservation by symplectic Runge-Kutta method (SRK) and symplectic partitioned Runge-Kutta method (SPRK) of structures for fixed points of linear Hamiltonian systems. The structure-preservation region provides a practical criterion for choosing step-size in symplectic computation. Examples are given to justify the investigation.

This paper presents a family of Runge{Kutta type integration schemes of arbitrarily high order for di erential equations evolving on manifolds. We prove that any classical Runge{Kutta method can be turned into an invariant method of the same order on a general homogeneous manifold, and present a family of algorithms that are relatively simple to implement.