Numerical construction of Gaussian quadrature formulas for ∫₀¹(-𝐿𝑜𝑔𝑥)⋅𝑥^{𝛼}⋅𝑓(𝑥)⋅𝑑𝑥 and ∫₀^{∞}𝐸_{𝑚}(𝑥)⋅𝑓(𝑥)⋅𝑑𝑥

Abstract

Most nonclassical Gaussian quadrature rules are difficult to construct because of the loss of significant digits during the generation of the associated orthogonal polynomials. But, in some particular cases, it is possible to develop stable algorithms. This is true for at least two well-known integrals, namely \[ ∫ 0 1 − ( Log x ) ⋅ x α ⋅ f ( x ) ⋅ d x and ∫ 0 ∞ E m ( x ) ⋅ f ( x ) ⋅ d x . \int _0^1 { - ({\operatorname {Log}}\;x) \cdot {x^\alpha } \cdot f(x) \cdot dx\quad {\text {and}}\quad \int _0^\infty {{E_m}(x) \cdot f(x) \cdot } dx.} \] A new approach is presented, which makes use of known classical Gaussian quadratures and is remarkably well-conditioned since the generation of the orthogonal polynomials requires only the computation of discrete sums of positive quantities. Finally, some numerical results are given.