laguerre_product


laguerre_product, an Octave code which computes weighted integrals of products of Laguerre polynomials.

LAGUERRE_LINEAR_PRODUCT takes as input the maximum degree P and the (nonnegative integer) exponent E and computes

Tij = Integral ( 0 <= X < +oo ) X^E * L(i)(X) * L(j)(X) * exp(-X) dx
for I and J between 0 and P.

LAGUERRE_EXPONENTIAL_PRODUCT takes as input the maximum degree P and the coefficient B and computes

Tij = Integral ( 0 <= X < +oo ) exp(B*X) * L(i)(X) * L(j)(X) * exp(-X) dx
for I and J between 0 and P.

When the polynomial chaos expansion is used to study stochastic differential equations, it is a common task to have to form and compute integrals of the sort considered here.

Note that, because of the orthonormality of the Laguerre polynomials, LAGUERRE_LINEAR_PRODUCT will return the identity matrix when E=0, and LAGUERRE_EXPONENTIAL_PRODUCT will return the identity matrix when B=0.

What is more interesting is that, because of the recursion relationship

i*L(i+1)(X) = (2*i-1 - X) * L(i)(X) -(i-1)*L(i-1)(X)
LAGUERRE_LINEAR_PRODUCT will return a symmetric tridiagonal matrix when E=1.

Licensing:

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

Languages:

laguerre_product is available in a MATLAB version and an Octave version.

Related Data and Programs:

laguerre_product_test

hermite_polynomial, an Octave code which evaluates the physicist's hermite polynomial, the probabilist's hermite polynomial, the hermite function, and related functions.

laguerre_rule, an Octave code which computes a gauss-laguerre quadrature rule.

legendre_product, an Octave code which computes weighted integrals of products of legendre polynomials.

polpak, an Octave code which evaluates a variety of mathematical functions.

quad_rule, an Octave code which defines quadrature rules on a variety of intervals with different weight functions.

Reference:

  1. Milton Abramowitz, Irene Stegun,
    Handbook of Mathematical Functions,
    National Bureau of Standards, 1964,
    ISBN: 0-486-61272-4,
    LC: QA47.A34.
  2. Philip Davis, Philip Rabinowitz,
    Methods of Numerical Integration,
    Second Edition,
    Dover, 2007,
    ISBN: 0486453391,
    LC: QA299.3.D28.
  3. Daniel Zwillinger, editor,
    CRC Standard Mathematical Tables and Formulae,
    30th Edition,
    CRC Press, 1996,
    ISBN: 0-8493-2479-3,
    LC: QA47.M315.

Source Code:


Last modified on 09 February 2019.