Quadrature Rules of Gauss-Laguerre Type

**QUADRATURE_RULES_LAGUERRE**
is a dataset directory which
contains examples of quadrature rules of Gauss-Laguerre type.

Gauss-Laguerre quadrature rules are designed to approximate integrals on semi-infinite intervals. By default, this interval is [0,+oo) but in some cases, the interval may be specified with a different initial point, which we will suggest by [A,+oo).

In the following discussion, we will represent the left endpoint by
**A**, while understanding that it is usually taken to be 0.

Standard Gauss-Laguerre quadrature assumes that the integrand we are considering has a form like:

Integral ( A <= x < +oo ) exp(-x) f(x) dxwhere the factor

A *standard Gauss-Laguerre quadrature rule* is a set of **n**
positive weights **w** and abscissas **x** so that

Integral ( A <= x < +oo ) exp(-x) f(x) dxmay be approximated by

Sum ( 1 <= I <= N ) w(i) * f(x(i))

It is often convenient to consider approximating integrals in which
the weighting factor **exp(-x)** is implicit. In that case, we
are looking at approximating

Integral ( A <= x < +oo ) f(x) dxand it is easy to modify a standard Gauss-Laguerre quadrature rule to handle this case directly.

A *modified Gauss-Laguerre quadrature rule* is a set of **n**
positive weights **w** and abscissas **x** so that

Integral ( A <= x < +oo ) f(x) dxmay be approximated by

Sum ( 1 <= I <= N ) w(i) * f(x(i))

When using a Gauss-Laguerre quadrature rule, it's important to know whether the rule has been developed for the standard or modified cases. Basically, the only change is that the weights of the modified rule have been multiplied by an exponential factor evaluated at the corresponding abscissa.

For this directory, a quadrature rule is stored as three files, containing the weights, the points, and a file containing two points defining the endpoints of the region. The first endpoint is either 0 or A; the second endpoint is really infinity, and we will simply place a "large" value here, trusting that the user will already know that this is merely a placeholder for infinity.

We consider a modified Gauss-Laguerre quadrature rule of order 4.

Here is the text of the "W" file storing the weights of such a rule:

```
0.60315
0.35741
0.03888
0.00053
```

(We have reduced the number of digits printed in this presentation)
Here is the text of the "X" file storing the abscissas of such a rule:

```
0.32254
1.74576
4.53662
9.39507
```

Here is the text of the "R" file storing the lower and upper limits of the region:

```
0.0
1.0E+30
```

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

INT_EXACTNESS_LAGUERRE, a C++ program which tests the polynomial exactness of Gauss-Laguerre quadrature rules.

LAGUERRE_POLYNOMIAL, a C++ library which evaluates the Laguerre polynomial, the generalized Laguerre polynomial, and the Laguerre function.

LAGUERRE_RULE, a C++ program which can compute and print a Gauss-Laguerre quadrature rule.

QUADRATURE_RULES_GEN_LAGUERRE, a dataset directory which contains quadrature rules for integration on a semi-infinite interval, using a generalized Gauss-Laguerre rule.

TEST_INT_LAGUERRE, a C++ library which defines test integrands for integration over [A,+oo).

Standard Gauss-Laguerre Rules, Order 1 through 10, 16, 32, 64:

- lag_o01_x.txt, the abscissas.
- lag_o01_w.txt, the weights.
- lag_o01_r.txt, the range.
- lag_o02_x.txt, the abscissas.
- lag_o02_w.txt, the weights.
- lag_o02_r.txt, the range.
- lag_o03_x.txt, the abscissas.
- lag_o03_w.txt, the weights.
- lag_o03_r.txt, the range.
- lag_o04_x.txt, the abscissas.
- lag_o04_w.txt, the weights.
- lag_o04_r.txt, the range.
- lag_o05_x.txt, the abscissas.
- lag_o05_w.txt, the weights.
- lag_o05_r.txt, the range.
- lag_o06_x.txt, the abscissas.
- lag_o06_w.txt, the weights.
- lag_o06_r.txt, the range.
- lag_o07_x.txt, the abscissas.
- lag_o07_w.txt, the weights.
- lag_o07_r.txt, the range.
- lag_o08_x.txt, the abscissas.
- lag_o08_w.txt, the weights.
- lag_o08_r.txt, the range.
- lag_o09_x.txt, the abscissas.
- lag_o09_w.txt, the weights.
- lag_o09_r.txt, the range.
- lag_o10_x.txt, the abscissas.
- lag_o10_w.txt, the weights.
- lag_o10_r.txt, the range.
- lag_o16_x.txt, the abscissas.
- lag_o16_w.txt, the weights.
- lag_o16_r.txt, the range.
- lag_o32_x.txt, the abscissas.
- lag_o32_w.txt, the weights.
- lag_o32_r.txt, the range.
- lag_o64_x.txt, the abscissas.
- lag_o64_w.txt, the weights.
- lag_o64_r.txt, the range.

Modified Gauss-Laguerre Rules, Order 1, 2, 4, 8, 16, 32, 64:

- lag_o01_modified_x.txt, the abscissas.
- lag_o01_modified_w.txt, the weights.
- lag_o01_modified_r.txt, the range.
- lag_o02_modified_x.txt, the abscissas.
- lag_o02_modified_w.txt, the weights.
- lag_o02_modified_r.txt, the range.
- lag_o04_modified_x.txt, the abscissas.
- lag_o04_modified_w.txt, the weights.
- lag_o04_modified_r.txt, the range.
- lag_o08_modified_x.txt, the abscissas.
- lag_o08_modified_w.txt, the weights.
- lag_o08_modified_r.txt, the range.
- lag_o16_modified_x.txt, the abscissas.
- lag_o16_modified_w.txt, the weights.
- lag_o16_modified_r.txt, the range.
- lag_o32_modified_x.txt, the abscissas.
- lag_o32_modified_w.txt, the weights.
- lag_o32_modified_r.txt, the range.
- lag_o64_modified_x.txt, the abscissas.
- lag_o64_modified_w.txt, the weights.
- lag_o64_modified_r.txt, the range.

You can go up one level to the DATASETS page.