Abstract
Two-parameter families of predictor-corrector methods based upon a combination of Adams- and Nyström formulae have been developed. The combinations use correctors of order one higher than that of the predictors. The methods are chosen to give optimal stability properties with respect to a requirement on the form and size of the regions of absolute stability. The optimal methods are listed and their regions of absolute stability are presented. The efficiency of the methods is compared to that of the corresponding Adams methods through numerical results from a variable order, variable stepsize program package.
Similar content being viewed by others
References
R. L. Crane and R. M. Klopfenstein,A predictor-corrector algorithm with an increased range of absolute stability, J.A.C.M. 12 (1965), 227–241.
D. M. Himmelblau,Applied Nonlinear Programming, McGraw-Hill, New York (1972).
R. Hook and T. A. Jeeves,“Direct search” solution of numerical and statistical problems, J.A.C.M. 8 (1961), 212–229.
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.
R. W. Klopfenstein and R. S. Millman,Numerical stability of one evaluation predictor-corrector algorithm for numerical solution of ordinary differential equations, Math. Comp. 22 (1968), 557–564.
F. T. Krogh,Predictor-corrector methods of high order with improved stability characteristics, J.A.C.M. 13 (1966), 374–385.
F. T. Krogh,On testing a subroutine for the numerical integration of ordinary differential equations, J.A.C.M. 20 (1973), 545–562.
J. D. Lambert,Computational methods in ordinary differential equations, John Wiley and Sons, London (1973).
J. J. H. Miller,On the location of zeros of certain classes of polynomials with applications to numerical analysis, J. Inst. Math. Appl. 8 (1971), 397–406.
K. Scheen,Fifth and sixth order PECE algorithms with improved stability properties, SIAM J. Numer Anal. 8,2 (1971), 244–248.
L. F. Shampine and M. K. Gordon,Computer solution of ordinary differential equations. The initial value problem, W. H. Freeman and Company, San Francisco (1975).
L. F. Shampine,Stiffness and non-stiff differential equation solvers, inNumerische Behandlung von Differentialgleichungen (L. Collatz, ed.). Int. Series Numer. Math. 27 (1975), 287–301.
P. G. Thomsen, and Z. Zlatev,A two-parameter family of multistep methods, Institute for Numerical Analysis, Technical University of Denmark, Rep. No. NI-75-03 (1975).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Thomsen, P.G., Zlatev, Z. Two-parameter families of predictor-corrector methods for the solution of ordinary differential equations. BIT 19, 503–517 (1979). https://doi.org/10.1007/BF01931267
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01931267