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:

${\displaystyle \int _{0}^{+\infty }e^{-x}f(x)\,dx.}$

In this case

${\displaystyle \int _{0}^{+\infty }e^{-x}f(x)\,dx\approx \sum _{i=1}^{n}w_{i}f(x_{i})}$

where x_i is the i-th root of Laguerre polynomial L_n(x) and the weight w_i is given by

${\displaystyle w_{i}={\frac {x_{i}}{\left(n+1\right)^{2}\left[L_{n+1}\left(x_{i}\right)\right]^{2}}}.}$

The following Python code with the SymPy library will allow for calculation of the values of ${\displaystyle x_{i}}$ and ${\displaystyle w_{i}}$ to 20 digits of precision:

from sympy import *

def lag_weights_roots(n):
    x = Symbol("x")
    roots = Poly(laguerre(n, x)).all_roots()
    x_i = [rt.evalf(20) for rt in roots]
    w_i = [(rt / ((n + 1) * laguerre(n + 1, rt)) ** 2).evalf(20) for rt in roots]
    return x_i, w_i

print(lag_weights_roots(5))