On the construction of Gaussian quadrature rules from modified moments.
HTML articles powered by AMS MathViewer
- by Walter Gautschi PDF
- Math. Comp. 24 (1970), 245-260 Request permission
Abstract:
Given a weight function $\omega (x)$ on $(\alpha ,\beta )$, and a system of polynomials $\left \{ {{p_k}(x)} \right \}_{k = 0}^\infty$, with degree ${p_k}(x) = k$, we consider the problem of constructing Gaussian quadrature rules $\int _\alpha ^\beta {f(x)\omega (x)dx = \sum \nolimits _{r = 1}^n {{\lambda _r}^{(n)}f({\xi _r}^{(n)})} }$ from "modified moments" ${v_k} = \int _\alpha ^\beta {{p_k}(x)\omega (x)dx}$. Classical procedures take ${p_k}(x) = {x^k}$, but suffer from progressive ill-conditioning as $n$ increases. A more recent procedure, due to Sack and Donovan, takes for $\{ {p_k}(x)\}$ a system of (classical) orthogonal polynomials. The problem is then remarkably well-conditioned, at least for finite intervals $[\alpha ,\beta ]$. In support of this observation, we obtain upper bounds for the respective asymptotic condition number. In special cases, these bounds grow like a fixed power of $n$. We also derive an algorithm for solving the problem considered, which generalizes one due to Golub and Welsch. Finally, some numerical examples are presented.References
- A. Erdélyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Tables of integral transforms. Vol. I, McGraw-Hill Book Co., Inc., New York-Toronto-London, 1954. Based, in part, on notes left by Harry Bateman. MR 0061695
- Géza Freud, Orthogonale Polynome, Lehrbücher und Monographien aus dem Gebiete der Exakten Wissenschaften, Mathematische Reihe, Band 33, Birkhäuser Verlag, Basel-Stuttgart, 1969 (German). MR 0481888, DOI 10.1007/978-3-0348-7169-3 D. Galant, "Gauss quadrature rules for the evaluation of $2{\pi ^{ - 1/2}}\smallint _0^\infty \exp ( - {x^2})f(x)dx$," Math. Comp., v. 23, 1969, p. 674.
- Walter Gautschi, On inverses of Vandermonde and confluent Vandermonde matrices. II, Numer. Math. 5 (1963), 425–430. MR 164437, DOI 10.1007/BF01385906 W. Gautschi, "Algorithm 236—Bessel functions of the first kind," Comm. ACM, v. 7, 1964, pp. 479–480.
- Walter Gautschi, Construction of Gauss-Christoffel quadrature formulas, Math. Comp. 22 (1968), 251–270. MR 228171, DOI 10.1090/S0025-5718-1968-0228171-0 W. Gautschi, "Algorithm 331—Gaussian quadrature formulas," Comm. ACM, v. 11, 1968, pp. 432–436.
- Gene H. Golub and John H. Welsch, Calculation of Gauss quadrature rules, Math. Comp. 23 (1969), 221-230; addendum, ibid. 23 (1969), no. 106, loose microfiche suppl, A1–A10. MR 0245201, DOI 10.1090/S0025-5718-69-99647-1
- I. P. Mysovskih, On the construction of cubature formulas with the smallest number of nodes, Dokl. Akad. Nauk SSSR 178 (1968), 1252–1254 (Russian). MR 0224284
- G. I. Natanson, A two-sided estimate for the Lebesgue function of the Lagrange interpolation process with Jacobi nodes, Izv. Vysš. Učebn. Zaved. Matematika 1967 (1967), no. 11 (66), 67–74 (Russian). MR 0221158
- I. P. Natanson, Konstruktivnaya teoriya funkciĭ, Gosudarstvennoe Izdatel′stvo Tehniko-Teoretičeskoĭ Literatury, Moscow-Leningrad, 1949 (Russian). MR 0034464 R. A. Sack & A. F. Donovan, An Algorithm for Gaussian Quadrature Given Generalized Moments, Department of Mathematics, University of Salford, Salford, England, 1969.
- N. M. Steen, G. D. Byrne, and E. M. Gelbard, Gaussian quadratures for the integrals $_{0}^{\infty }\,\textrm {exp}(-x^{2})f(x)dx$ and $_{0}^{b}\,\textrm {exp}(-x^{2})f(x)dx$, Math. Comp. 23 (1969), 661–671. MR 247744, DOI 10.1090/S0025-5718-1969-0247744-3 G. Szegö, Orthogonal Polynomials, 2nd rev. ed., Amer. Math. Soc. Colloq. Publ., vol. 23, Amer. Math. Soc., Providence, R. I., 1959. MR 21 #5029.
- Ion Zamfirescu, An extension of Gauss’ method for the calculation of improper integrals, Acad. R. P. Romîne. Stud. Cerc. Mat. 14 (1963), 615–631 (Romanian, with French and Russian summaries). MR 184434
Additional Information
- © Copyright 1970 American Mathematical Society
- Journal: Math. Comp. 24 (1970), 245-260
- MSC: Primary 65.55
- DOI: https://doi.org/10.1090/S0025-5718-1970-0285117-6
- MathSciNet review: 0285117