Multiple quadrature with central differences on one line
HTML articles powered by AMS MathViewer
- by Herbert E. Salzer PDF
- Math. Comp. 16 (1962), 244-248 Request permission
Abstract:
The coefficients $A_{2m}^n$ in the n-fold quadrature formulas for the stepwise integration of (1) ${y^{(n)}} = f(x,y,y’, \cdots , {y^{(n - 1)}})$, at intervals of h, namely, for n even, (2) ${\delta ^n}{y_0} = {h^n}\sum \nolimits _{m = 1}^{10} {(1 + A_{2m}^n{\delta ^{2m}}){f_0} + \cdots }$, for n odd, (3) $\mu {\delta ^n}{y_0} = {h^n}\sum \nolimits _{m = 1}^{10} {(1 + A_{2m}^n{\delta ^{2m}}){f_0} + \cdots }$, are tabulated exactly for n = 1(1)6, m = 1(1)10. They were calculated from the well-known symbolic formulas (4) ${{\delta }^{n}}y = {{({\delta }/{D})}^{n}}f$, (5) ${{({\delta }/{D})}^{n}} = {{({{{\delta h}/{2 \sinh }}^{-1}}({\delta }/{2}))}^{n}}$ and (6) $\mu = {{(1 + {{{\delta }^{2}}}/{4})}^{{1}/{2}}} = 1 + \frac {{{\delta }^{2}}}{8} - \frac {{{\delta }^{4}}}{128} + \frac {{{\delta }^{6}}}{1024} - \frac {5{{\delta }^{8}}}{32768} + \cdots$. For calculating ${y^{(r)}}$, replace n by n - r in (2) and (3). Use of (2) and (3) avoids the solution of (1) by simultaneous lower-order systems for $n > 1$, as well as mid-interval tabular arguments, requires only even-order differences, on a single line, and provides great accuracy due to rapid decrease of $A_{2m}^n$ as m increases. However, the integration may be slowed down by the need to estimate and refine iteratively the later values of $y,y’, \cdots , {y^{(n - 1)}}$ required in ${\delta ^{2m}}{f_0}$. Reference to earlier collected formulas of Legendre, Oppolzer, Thiele, Lindow, Salzer, Milne and Buckingham, reveals that Thiele and Buckingham come closest to (2), (3), as their works contain schemes that involve just tabular arguments throughout. For n odd, they give formulas that are based upon the series in ${\delta ^{2m}}$ for $({1}/{\mu }){{({\delta }/{D})}^{n}}$ instead of $\mu {{({\delta }/{D})}^{n}}$ as in the present arrangement.References
-
A. M. Legendre, Traité des Fonctions Elliptiques, v. 2, Paris, 1826, Chapter 3, p. 41-60 (For errors, see MTAC, v. 5, 1951, p. 27).
T. R. Oppolzer, Lehrbuch zur Bahnbestimmung der Cometen und Planeten, v. 2, W. Engelmann, Leipzig, 1880, p. 35, 53-54, 545, 596.
T. N. Thiele, Interpolationsrechnung, B. G. Teubner, Leipzig, 1909, p. 95-97. (Some misprints are noted in Math. Comp., v. 15, 1961, p. 321.)
M. Lindow, Numerische Infinitesimalrechnung, F. Dümmler, Berlin and Bonn, 1928, p. 170-171.
- Herbert E. Salzer, Coefficients for mid-interval numerical integration with central differences, Philos. Mag. (7) 36 (1945), 216–218. MR 13924 W. E. Milne, Numerical Calculus, Princeton, 1949, p. 196-197.
- Hebert E. Salzer, Coefficients for repeated integration with central differences, J. Math. Physics 28 (1949), 54–61. MR 0029284 R. A. Buckingham, Numerical Methods, Pitman Publishing Corp., New York and London, 1957, p. 150-154. (For errors, see Math. Comp., v. 15, 1961, p. 319.).
Additional Information
- © Copyright 1962 American Mathematical Society
- Journal: Math. Comp. 16 (1962), 244-248
- MSC: Primary 65.55
- DOI: https://doi.org/10.1090/S0025-5718-1962-0145655-9
- MathSciNet review: 0145655