Zusammenfassung: | Given a weight function ω(x) on (α, β), and a system of polynomials {p<sub>k</sub>(x)}<sup>∞</sup><sub>k = 0</sub>, with degree p<sub>k</sub>(x) = k, we consider the problem of constructing Gaussian quadrature rules <tex-math>$\int^\beta_\alpha f(x)\omega(x)dx \doteq = \sum^n_{r = 1} \lambda^{(n)}_r f(\xi^{(n)}_r}$</tex-math> from "modified moments" ν<sub>k</sub> = ∫<sup>β</sup><sub>α</sub> p<sub>k</sub>(x)ω(x)dx. Classical procedures take p<sub>k</sub>(x) = x<sup>k</sup>, but suffer from progressive ill-conditioning as n increases. A more recent procedure, due to Sack and Donovan, takes for {p<sub>k</sub>(x)} a system of (classical) orthogonal polynomials. The problem is then remarkably well-conditioned, at least for finite intervals [ α, β ]. 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.
|