Search

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.

We study error propagation in both an explicit and an implicit method for solving Volterra integro-differential equations. We determine the relationship between local and global errors. We derive upper bounds for the global error, and show that the global order for both methods is expected to be first-order. A few numerical examples illustrate our results., Comment: 14p, 5 figs

We present implicit and explicit versions of a numerical algorithm for solving a Volterra integro-differential equation. These algorithms are an extension of our previous work, and cater for a kernel of general form. We use an appropriate test equation to study the stability of both algorithms, numerically deriving stability regions. The region for the implicit method appears to be unbounded, while the explicit has a bounded region close to the origin. We perform a few calculations to demonstrate our results., Comment: 10 pages, 1 Figure

We show how the error term for the Trapezium Rule can be estimated, by solving an initial value problem using a Runge-Kutta method. The error term can then be added to the Trapezium approximation, yielding a much more accurate result. We also show how the risk of singularities in the relevant initial value problem can be mitigated.

We describe an algorithm, based on Euler's method, for solving Volterra integro-differential equations. The algorithm approximates the relevant integral by means of the composite Trapezium Rule, using the discrete nodes of the independent variable as the required nodes for the integration variable. We have developed an error control device, using Richardson extrapolation, and we have achieved accuracy better than 1e-12 for all numerical examples considered., Comment: 11 pages

We transform a double integral into a second-order initial value problem, which we solve using Euler's method and Richardson extrapolation. For an example we consider, we achieve accuracy close to machine precision (1e-15). We also use the algorithm to determine the error curve for a Simpson cubature rule., Comment: 17p, 4 figs

We report on a novel algorithm for controlling global error in a step-by-step (stepwise) sense, in the numerical solution of a scalar, autonomous, nonstiff or weakly stiff problem. The algorithm exploits the remainder term of a Taylor expansion of the solution. It requires the use of the DP853 triple to solve an auxiliary problem which, in turn, enables the remainder term to be determined. A quenching process then allows the solution generated by Euler's method to be controlled. We have achieved tolerances on the relative global error as strict as 1e-10., Comment: 32 pages, 11 figures

We determine the pointwise error in Hermite interpolation by numerically solving an appropriate differential equation, derived from the error term itself. We use this knowledge to approximate the error term by means of a polynomial, which is then added to the original Hermite polynomial to form a more accurate approximation. An example demonstrates that improvements in accuracy are significant., Comment: 10 pages, 3 figures

We study the efficiency of the multisection method for univariate nonlinear equations, relative to that for the well-known bisection method. We show that there is a minimal effort algorithm that uses more sections than the bisection method, although this optimal algorithm is problem dependent. The number of sections required for optimality is determined by means of a Lambert W function., Comment: 12 pages, 4 figures

We determine the Lagrange function in Taylor polynomial approximation by solving an appropriate initial-value problem. Hence, we determine the remainder term which we then approximate by means of a natural cubic spline. This results in a significant improvement in the quality of the Taylor approximation. We observe improvements in the accuracy of the approximation of many orders of magnitude, including a case when the independent variable x lies beyond the relevant radius of convergence., Comment: 14 pages, 6 figures

We determine the Rolle function in Lagrange polynomial approximation using a suitable differential equation. We then propose a device for improving the Lagrange approximation by exploiting our knowledge of the Rolle function., Comment: 9 pages including 3 pages of figures

We study the effect of global error control in the numerical solution of Hamiltonian systems. In particular, we apply the RKQ algorithm in the numerical solution of a Hamiltonian system. This algorithm is designed to provide stepwise control of both local and global error. A test problem demonstrates the error control features of RKQ. Good results are obtained, despite the fact that explicit Runge-Kutta methods have been used in RKQ, rather than symplectic Runge-Kutta methods. This simply emphasizes the value of stepwise global error control, as per the RKQ algorithm.

We incorporate explicit Nystrom methods into the RKQ algorithm for stepwise global error control in numerical solutions of initial-value problems. The initial-value problem is transformed into an explicitly second-order problem, so as to be suitable for Nystrom integration. The Nystrom methods used are fourth-order, fifth-order and 10th-order. Two examples demonstrate the effectiveness of the algorithm., Comment: This is an extension of ideas published in J. Math. Res. (open access); see refs [1] and [2]

We study the relationship between local and global error in Runge-Kutta methods for initial-value problems in ordinary differential equations. We show that local error control by means of local extrapolation does not equate to global error control. Our analysis shows that the global error of the higher-order solution is propagated under iteration, and this can cause an uncontrolled increase in the global error of the lower-order solution. We find conditions under which global error control occurs during the initial stages of the RK integration, but even in such a case the global error is likely to eventually exceed the user-defined tolerance., Comment: Corrections to some absolute values signs in Section 3.2