To evaluate the class of integrals $\int^1_{-1}e^{-\alpha x}f(x) dx$, where $\alpha \in \R^+$ and the function f(x) is known only approximately in a tabular form, we wish to use a Gaussian quadrature formula. Nodes and weights have to be computed using the family of monic orthogonal polynomials, with respect to the weight function $e^{-\alpha x}$, obtained through the three-term recurrence relation $P_{k+1}(x) = (x+B_{k+1})P_k(x)-C_{k+1}P_{k-1}(x)$. To guarantee a good precision, we must evaluate carefully the values for the coefficients $B_{k+1}$ and $C_{k+1}$. Such evaluations are made completely formally through a Mathematica program to obtain great precision. A comparison between various methods, starting from moments and modified moments, is shown. Numerical results are also presented.
Computing the coefficients of a recurrence formula for numerical integration by moments and modified moments
REDIVO ZAGLIA, MICHELA
1993
Abstract
To evaluate the class of integrals $\int^1_{-1}e^{-\alpha x}f(x) dx$, where $\alpha \in \R^+$ and the function f(x) is known only approximately in a tabular form, we wish to use a Gaussian quadrature formula. Nodes and weights have to be computed using the family of monic orthogonal polynomials, with respect to the weight function $e^{-\alpha x}$, obtained through the three-term recurrence relation $P_{k+1}(x) = (x+B_{k+1})P_k(x)-C_{k+1}P_{k-1}(x)$. To guarantee a good precision, we must evaluate carefully the values for the coefficients $B_{k+1}$ and $C_{k+1}$. Such evaluations are made completely formally through a Mathematica program to obtain great precision. A comparison between various methods, starting from moments and modified moments, is shown. Numerical results are also presented.File | Dimensione | Formato | |
---|---|---|---|
10.1016-0377-0427(93)90152-2.pdf
accesso aperto
Tipologia:
Published (publisher's version)
Licenza:
Accesso libero
Dimensione
650.15 kB
Formato
Adobe PDF
|
650.15 kB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.