QUADRATURE_RULES_GEN_LAGUERRE
Quadrature Rules
Generalized Gauss-Laguerre Type


QUADRATURE_RULES_GEN_LAGUERRE is a dataset directory which contains examples of quadrature rules of Generalized Gauss-Laguerre type.

Generalized Gauss-Legendre quadrature rules are a generalization of the basic Gauss-Legendre quadrature rule. The generalization involves including a factor of xalpha in the weight factor. (The parameter alpha is a real value which may take on any value greater than -1.)

Generalized 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 generalized Gauss-Laguerre quadrature assumes that the integrand we are considering has a form like:

        Integral ( A <= x < oo ) xalpha * exp(-x) * f(x) dx
      
where the factor xalpha * exp(-x) is regarded as a weight factor.

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

        Integral ( A <= x < oo ) xalpha * exp(-x) * f(x) dx
      
may 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 xalpha * exp(-x) is implicit. In that case, we are looking at approximating

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

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

        Integral ( A <= x < oo ) f(x) dx
      
may be approximated by
        Sum ( 1 <= I <= N ) w(i) * f(x(i))
      

When using a generalized 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 the inverse of the weight function 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.

Example:

We consider a standard generalized Gauss-Laguerre quadrature rule of order 4, with ALPHA = 0.5.

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


        0.4530087825780014    
        0.3816169905147275    
        0.5079463161099657E-01
        0.8065912142432840E-03
      

Here is the text of the "X" file storing the abscissas of such a rule:


        0.5235260767382691    
        2.156648763269094    
        5.137387546176711    
        10.18243761381592    
      

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


        0.0
        1.0E+30
      

The modified rule of order 4 and ALPHA = 0.5 has the same values of X and R, but the modified weights are now:


        1.056822251392305    
        2.245735732700531    
        3.815800963697504    
        6.681963149567727  
      

Licensing:

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

Related Data and Programs:

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

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

Sample Files:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = -0.5:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 0.0:

Standard Generalized Gauss-Laguerre Rule, Order 1, ALPHA = 0.5:

Standard Generalized Gauss-Laguerre Rule, Order 2, ALPHA = 0.5:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 0.5:

Standard Generalized Gauss-Laguerre Rule, Order 8, ALPHA = 0.5:

Standard Generalized Gauss-Laguerre Rule, Order 16, ALPHA = 0.5:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 1.0:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 2.0:

Standard Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 3.0:

Modified Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 0.0:

Modified Generalized Gauss-Laguerre Rule, Order 1, ALPHA = 0.5:

Modified Generalized Gauss-Laguerre Rule, Order 2, ALPHA = 0.5:

Modified Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 0.5:

Modified Generalized Gauss-Laguerre Rule, Order 8, ALPHA = 0.5:

Modified Generalized Gauss-Laguerre Rule, Order 16, ALPHA = 0.5:

Modified Generalized Gauss-Laguerre Rule, Order 4, ALPHA = 1.0:

You can go up one level to the DATASETS page.


Last revised on 04 September 2007.