A Family of Fifth-Order Runge-Kutta Pairs

Authors:
S. N. Papakostas and G. Papageorgiou

Journal:
Math. Comp. **65** (1996), 1165-1181

MSC (1991):
Primary 65L05

MathSciNet review:
1333323

Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: The construction of a Runge-Kutta pair of order with the minimal number of stages requires the solution of a nonlinear system of order conditions in unknowns. We define a new family of pairs which includes pairs using function evaluations per integration step as well as pairs which additionally use the first function evaluation from the next step. This is achieved by making use of Kutta's simplifying assumption on the original system of the order conditions, i.e., that all the internal nodes of a method contributing to the estimation of the endpoint solution provide, at these nodes, cost-free second-order approximations to the true solution of any differential equation. In both cases the solution of the resulting system of nonlinear equations is completely classified and described in terms of five free parameters. Optimal Runge-Kutta pairs with respect to minimized truncation error coefficients, maximal phase-lag order and various stability characteristics are presented. These pairs were selected under the assumption that they are used in Local Extrapolation Mode (the propagated solution of a problem is the one provided by the fifth-order formula of the pair). Numerical results obtained by testing the new pairs over a standard set of test problems suggest a significant improvement in efficiency when using a specific pair of the new family with minimized truncation error coefficients, instead of some other existing pairs.

**1.**J. C. Butcher,*On Runge-Kutta processes of high order*, J. Austral. Math. Soc.**4**(1964), 179–194. MR**0165692****2.**J. C. Butcher,*The numerical analysis of ordinary differential equations*, A Wiley-Interscience Publication, John Wiley & Sons, Ltd., Chichester, 1987. Runge\mhy Kutta and general linear methods. MR**878564****3.**C. R. Cassity,*The complete solution of the fifth order Runge-Kutta equations*, SIAM J. Numer. Anal.**6**(1969), 432–436. MR**0255061****4.**A. R. Curtis,*High-order explicit Runge-Kutta formulae, their uses, and limitations*, J. Inst. Math. Appl.**16**(1975), no. 1, 35–55. MR**0383750****5.**J. R. Dormand and P. J. Prince,*A family of embedded Runge-Kutta formulae*, J. Comput. Appl. Math.**6**(1980), no. 1, 19–26. MR**568599**, 10.1016/0771-050X(80)90013-3**6.**W. H. Enright and J. D. Pryce,*Two FORTRAN packages for assessing initial value methods*, ACM Trans. Math. Software**13**(1987), 1--27.**7.**E. Fehlberg,*Classical fifth, sixth, seventh, and eighth order Runge-Kutta formulas with stepsize control*, TR R--287, NASA, 1968.**8.**------,*Low order classical Runge-Kutta formulas with stepsize control and their application to some heat-transfer problems*, TR R--315, NASA, 1969.**9.**E. Hairer, S. P. Nørsett, and G. Wanner,*Solving ordinary differential equations. I*, Springer Series in Computational Mathematics, vol. 8, Springer-Verlag, Berlin, 1987. Nonstiff problems. MR**868663****10.**P. J. van der Houwen and B. P. Sommeijer,*Explicit Runge-Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions*, SIAM J. Numer. Anal.**24**(1987), no. 3, 595–617. MR**888752**, 10.1137/0724041**11.**T. E. Hull, W. H. Enright, B. M. Fellen, and A. E. Sedgwick,*Comparing numerical methods for ordinary differential equations*, SIAM J. Numer. Anal.**9**(1972), 603–637; errata, ibid. 11 (1974), 681. MR**0351086****12.**W. Kutta,*Beitrag zur näherungsweisen Integration totaler Differentialgleichungen*, Z. Math. Phys. (1901), 435--453.**13.**J. Douglas Lawson,*An order five Runge-Kutta process with extended region of stability*, SIAM J. Numer. Anal.**3**(1966), 593–597. MR**0216760****14.**G. Papageorgiou, Ch. Tsitouras, and S. N. Papakostas,*Runge-Kutta pairs for periodic initial value problems*, Rep. NA 93--1, Nat. Tech. Univ. Athens, Dept. Math., 1993.**15.**S. N. Papakostas, Ch. Tsitouras, and G. Papageorgiou,*A general family of explicit Runge-Kutta pairs of orders*, SIAM J. Numer. Anal.**33**(1996).**16.**L. F. Shampine,*Local extrapolation in the solution of ordinary differential equations*, Math. Comp.**27**(1973), 91–97. MR**0331803**, 10.1090/S0025-5718-1973-0331803-1**17.**E. Baylis Shanks,*Solutions of differential equations by evaluations of functions*, Math. Comp.**20**(1966), 21–38. MR**0187406**, 10.1090/S0025-5718-1966-0187406-1

Retrieve articles in *Mathematics of Computation of the American Mathematical Society*
with MSC (1991):
65L05

Retrieve articles in all journals with MSC (1991): 65L05

Additional Information

**S. N. Papakostas**

Affiliation:
Department of Mathematics, Zografou Campus, National Technical University of Athens, Athens 157 80, Greece

Email:
spapakos@theseas.ntua.gr

**G. Papageorgiou**

Affiliation:
Department of Mathematics, Zografou Campus, National Technical University of Athens, Athens 157 80, Greece

Email:
papag@nisyros.ntua.gr

DOI:
http://dx.doi.org/10.1090/S0025-5718-96-00718-1

Keywords:
Initial Value Problems,
Runge-Kutta,
pairs of embedded methods,
phase-lag

Received by editor(s):
September 7, 1993

Received by editor(s) in revised form:
September 5, 1994, and April 5, 1995

Article copyright:
© Copyright 1996
American Mathematical Society