Gauss–Laguerre quadrature

In numerical analysis Gauss–Laguerre quadrature (named after Carl Friedrich Gauss and Edmond Laguerre) is an extension of the Gaussian quadrature method for approximating the value of integrals of the following kind:

In this case

where xi is the i-th root of Laguerre polynomial Ln(x) and the weight wi is given by [1]

For more general functions

To integrate the function we apply the following transformation

where . For the last integral one then uses Gauss-Laguerre quadrature. Note, that while this approach works from an analytical perspective, it is not always numerically stable.

Generalized Gauss–Laguerre quadrature

More generally, one can also consider integrands that have a known power-law singularity at x=0, for some real number , leading to integrals of the form:

In this case, the weights are given[2] in terms of the generalized Laguerre polynomials:

where are the roots of .

This allows one to efficiently evaluate such integrals for polynomial or smooth f(x) even when α is not an integer.[3]

gollark: Recursive *string rewrites*?
gollark: The code is kind of horrifying.
gollark: Which is implemented in Rust.
gollark: osmarkscalculator™.
gollark: I just wrote these definitions, which are very definitive:``` ("Fac[0]", "1"), ("Fac[n]", "Fac[n-1]*n"), ("Fib[0]", "0"), ("Fib[1]", "1"), ("Fib[n]", "Fib[n-1] + Fib[n-2]")```

References

  1. Equation 25.4.45 in Abramowitz, M.; Stegun, I. A. Handbook of Mathematical Functions. Dover. ISBN 978-0-486-61272-0. 10th reprint with corrections.
  2. Weisstein, Eric W., "Laguerre-Gauss Quadrature" From MathWorld--A Wolfram Web Resource, Accessed March 9, 2020
  3. Rabinowitz, P.; Weiss, G. (1959). "Tables of Abscissas and Weights for Numerical Evaluation of Integrals of the form ". Mathematical Tables and Other Aids to Computation. 13: 285–294. doi:10.1090/S0025-5718-1959-0107992-3.

Further reading

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.