gauss_laguerre Subroutine

public subroutine gauss_laguerre(x, w)

Calculates the nodes and weights to perform Gauss-Laguerre quadrature on the interval [0, infinity)

NOTE: This subroutine has been tested for array lengths up to 64 and delivers an accuracy of 1e-12 for 64 array elements in double precision. For fewer elements the accuracy is improved.

Arguments

Type IntentOptional Attributes Name
real(kind=FP), intent(out) :: x(:)

Quadrature nodes
real(kind=FP), intent(out) :: w(:)

Quadrature weights