A high-order difference method for differential equations

Authors:
Robert E. Lynch and John R. Rice

Journal:
Math. Comp. **34** (1980), 333-372

MSC:
Primary 65L10

DOI:
https://doi.org/10.1090/S0025-5718-1980-0559190-8

MathSciNet review:
559190

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: This paper analyzes a high-accuracy approximation to the *m*th-order linear ordinary differential equation . At mesh points, *U* is the estimate of *u*; and *U* satisfies , where is a linear combination of values of *U* at stencil points (adjacent mesh points) and is a linear combination of values of *f* at *J* auxiliary points, which are between the first and last stencil points. The coefficients of , are obtained ``locally'' by solving a small linear system for each group of stencil points in order to make the approximation exact on a linear space *S* of dimension . For separated two-point boundary value problems, *U* is the solution of an *n*-by-*n* linear system with full bandwidth . For *S* a space of polynomials, existence and uniqueness are established, and the discretization error is the first divided differences of *U* tend to those of *u* at this rate. For a general set of auxiliary points one has ; but special auxiliary points, which depend upon *M* and the stencil points, allow larger *L*, up to . Comparison of operation counts for this method and five other common schemes shows that this method is among the most efficient for given convergence rate. A brief selection from extensive experiments is presented which supports the theoretical results and the practicality of the method.

**[G]**Garrett Birkhoff and Carl R. De Boor,*Piecewise polynomial interpolation and approximation*, Approximation of Functions (Proc. Sympos. General Motors Res. Lab., 1964 ), Elsevier Publ. Co., Amsterdam, 1965, pp. 164–190. MR**0189219****[G]**Garrett Birkhoff and Gian-Carlo Rota,*Ordinary differential equations*, Second edition, Blaisdell Publishing Co. Ginn and Co., Waltham, Mass.-Toronto, Ont.-London, 1969. MR**0236441****[R]**F. BOISVERT [1978],*The Effect on Accuracy of the Placement of Auxiliary Points in the HODIE Method for the Helmholtz Problem*, Dept. of Comput. Sci. Report CSD-TR 266, Purdue University, June.**[L]**Lothar Collatz,*The numerical treatment of differential equations. 3d ed*, Translated from a supplemented version of the 2d German edition by P. G. Williams. Die Grundlehren der mathematischen Wissenschaften, Bd. 60, Springer-Verlag, Berlin-Göttingen-Heidelberg, 1960. MR**0109436****[H]**H. B. Curry and I. J. Schoenberg,*On Pólya frequency functions. IV. The fundamental spline functions and their limits*, J. Analyse Math.**17**(1966), 71–107. MR**0218800**, https://doi.org/10.1007/BF02788653**[C]**Carl de Boor and Robert E. Lynch,*On splines and their minimum properties*, J. Math. Mech.**15**(1966), 953–969. MR**0203306****[C]**Carl de Boor and Blâir Swartz,*Collocation at Gaussian points*, SIAM J. Numer. Anal.**10**(1973), 582–606. MR**0373328**, https://doi.org/10.1137/0710052**[E]**J. DOEDEL [1976],*The Construction of Finite Difference Approximations to Ordinary Differential Equations*, Report, Dept. of Appl. Math., California Institute of Technology, Pasadena, California.**[E]**Eusebius J. Doedel,*The construction of finite difference approximations to ordinary differential equations*, SIAM J. Numer. Anal.**15**(1978), no. 3, 450–465. MR**0483481**, https://doi.org/10.1137/0715029**[S]**J. KARLIN &. W. J. STUDDEN [1966],*Tchebycheff Systems*:*With Applications in Analysis and Statistics*, Interscience, New York.**[R]**E. LYNCH & J. R. RICE [1975],*The HODIE Method*:*A Brief Introduction with Summary of Computational Properties*, Dept. of Comput. Sci. Report 170, Purdue University, Nov. 18.**[R]**Robert E. Lynch and John R. Rice,*High accuracy finite difference approximation to solutions of elliptic partial differential equations*, Proc. Nat. Acad. Sci. U.S.A.**75**(1978), no. 6, 2541–2544. MR**496774****[R]**E. LYNCH [1977a],*Accurate Finite Difference Approximation to Solutions of the Poisson Equation in Three Variables*, Dept. of Comput. Sci. Report CSD-TR 221, Purdue University, Feb. 15.**[R]**E. LYNCH [1977b],*Discretization Error Finite Difference Approximation to Solutions of the Poisson Equation in Three Variables*, Dept. of Comput. Sci. Report CSD-TR 230, Purdue University, April 19.**[R]**Robert E. Lynch and John R. Rice,*High accuracy finite difference approximation to solutions of elliptic partial differential equations*, Proc. Nat. Acad. Sci. U.S.A.**75**(1978), no. 6, 2541–2544. MR**496774****[R]**Robert E. Lynch and John R. Rice,*The Hodie method and its performance for solving elliptic partial differential equations*, Recent advances in numerical analysis (Proc. Sympos., Math. Res. Center, Univ. Wisconsin, Madison, Wis., 1978) Publ. Math. Res. Center Univ. Wisconsin, vol. 41, Academic Press, New York-London, 1978, pp. 143–175. MR**519061****[M]**M. R. Osborne,*Minimising truncation error in finite difference approximations to ordinary differential equations*, Math. Comp.**21**(1967), 133–145. MR**0223107**, https://doi.org/10.1090/S0025-5718-1967-0223107-X**[J]**James L. Phillips and Richard J. Hanson,*Gauss quadrature rules with 𝐵-spline weight functions*, Math. Comp.**28**(1974), no. 126, loose microfiche suppl., A1–C4. MR**0343551**, https://doi.org/10.2307/2005948**[H]**Carl de Boor (ed.),*Mathematical aspects of finite elements in partial differential equations*, Academic Press [A subsidiary of Harcourt Brace Jovanovich, Publishers], New York-London, 1974. Publication No. 33 of the Mathematics Research Center, The University of Wisconsin-Madison. MR**0349031****[R]**R. D. Russell and L. F. Shampine,*A collocation method for boundary value problems*, Numer. Math.**19**(1972), 1–28. MR**0305607**, https://doi.org/10.1007/BF01395926**[R]**R. D. Russell and J. M. Varah,*A comparison of global methods for linear two-point boundary value problems*, Math. Comput.**29**(1975), no. 132, 1007–1019. MR**0388788**, https://doi.org/10.1090/S0025-5718-1975-0388788-3**[R]**Robert D. Russell,*A comparison of collocation and finite differences for two-point boundary value problems*, SIAM J. Numer. Anal.**14**(1977), no. 1, 19–39. Papers on the numerical solution of two-point boundary-value problems (NSF-CBMS Regional Res. Conf., Texas Tech Univ., Lubbock, Tex., 1975). MR**0451745**, https://doi.org/10.1137/0714003**[B]**Carl de Boor (ed.),*Mathematical aspects of finite elements in partial differential equations*, Academic Press [A subsidiary of Harcourt Brace Jovanovich, Publishers], New York-London, 1974. Publication No. 33 of the Mathematics Research Center, The University of Wisconsin-Madison. MR**0349031**

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

Retrieve articles in all journals with MSC: 65L10

Additional Information

DOI:
https://doi.org/10.1090/S0025-5718-1980-0559190-8

Article copyright:
© Copyright 1980
American Mathematical Society