An extension of Olver's method for the numerical solution of linear recurrence relations

Author:
J. R. Cash

Journal:
Math. Comp. **32** (1978), 497-510

MSC:
Primary 65Q05

DOI:
https://doi.org/10.1090/S0025-5718-1978-0483578-8

MathSciNet review:
0483578

Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: An algorithm is developed for computing the solution of a class of linear recurrence relations of order greater than two when unstable error propagation prevents the required solution being found by direct forward recurrence. By abandoning an appropriate number of initial conditions the original problem may be replaced by an inexact but well-conditioned boundary value problem, and in certain circumstances the solution of this new problem is a good approximation to the required solution of the original problem. The required solution of this reposed problem is generated using an algorithm based on Gaussian elimination, and a technique developed by Olver is extended to estimate automatically the truncation error of the proposed algorithm.

**[1]**J. R. Cash,*High order methods for the numerical integration of ordinary differential equations*, Numer. Math.**30**(1978), no. 4, 385–409. MR**502523**, https://doi.org/10.1007/BF01398507**[2]**Walter Gautschi,*Computational aspects of three-term recurrence relations*, SIAM Rev.**9**(1967), 24–82. MR**0213062**, https://doi.org/10.1137/1009002**[3]**J. D. Lambert,*Computational methods in ordinary differential equations*, John Wiley & Sons, London-New York-Sydney, 1973. Introductory Mathematics for Scientists and Engineers. MR**0423815****[4]**D. W. LOZIER,*A Stable Algorithm for Computing any Solution of an Arbitrary Linear Difference Equation*, Ph.D. Thesis, University of Maryland.**[5]**J. C. P. MILLER,*British Association for the Advancement of Science*:*Bessel Functions, Part*II.*Mathematical Tables*, Vol. 10, Cambridge Univ. Press, Cambridge, 1952.**[6]**Kenneth S. Miller,*Linear difference equations*, W. A. Benjamin, Inc., New York-Amsterdam, 1968. MR**0227644****[7]**J. Oliver,*Relative error propagation in the recursive solution of linear recurrence relations*, Numer. Math.**9**(1966/1967), 323–340. MR**0213064**, https://doi.org/10.1007/BF02162423**[8]**J. Oliver,*The numerical solution of linear recurrence relations*, Numer. Math.**11**(1968), 349–360. MR**0226893**, https://doi.org/10.1007/BF02166688**[9]**J. Oliver,*An extension of Olver’s error estimation technique for linear recurrence relations*, Numer. Math.**12**(1968), 459–467. MR**0239782**, https://doi.org/10.1007/BF02161370**[10]**F. W. J. Olver,*Numerical solution of second-order linear difference equations*, J. Res. Nat. Bur. Standards Sect. B**71B**(1967), 111–129. MR**0221789****[11]**F. W. J. Olver,*Bounds for the solutions of second-order linear difference equations*, J. Res. Nat. Bur. Standards Sect. B**71B**(1967), 161–166. MR**0229407****[12]**J. H. Wilkinson,*The algebraic eigenvalue problem*, Clarendon Press, Oxford, 1965. MR**0184422****[13]**R. V. M. ZAHAR, "A mathematical analysis of Miller's algorithm,"*Numer. Math.*, v. 27, 1977, pp. 427-447.

Retrieve articles in *Mathematics of Computation*
with MSC:
65Q05

Retrieve articles in all journals with MSC: 65Q05

Additional Information

DOI:
https://doi.org/10.1090/S0025-5718-1978-0483578-8

Article copyright:
© Copyright 1978
American Mathematical Society