pce_legendre, a MATLAB code which sets up the system matrix for a polynomial chaos expansion, using Legendre polynomials with a linear factor, applied to a 2D partial differential equation (PDE) with a stochastic diffusion coefficient.

We wish to analyze a stochastic PDE of the form:

        -div A(X,Y) grad U(X,Y) = F(X)

We let X be a space of finite element functions generated by piecewise linear functions associated with a particular triangular dissection of the unit square.

We let Y be the space of polynomials over R^N with total degree at most P.

We seek solutions U in XxY using a polynomial chaos expansion approach.


[ nxy, bxy, fxy ] = pce_legendre_linear_assemble ( n, p )
where the input is: and the output is


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


pce_legendre is available in a MATLAB version.

Related Data and Programs:


cnoise, a MATLAB code which generates samples of noise obeying a 1/f^alpha power law, by Miroslav Stoyanov.

legendre_polynomial, a MATLAB code which evaluates the Legendre polynomial and associated functions.

legendre_rule, a MATLAB code which computes a Gauss-Legendre quadrature rule.

pce_burgers, a MATLAB code which defines and solves a version of the time-dependent viscous Burgers equation, with uncertain viscosity, using a polynomial chaos expansion, by Gianluca Iaccarino.

pce_ode_hermite, a MATLAB code which sets up a simple scalar ordinary differential equation (ODE) for exponential decay with an uncertain decay rate, using a polynomial chaos expansion in terms of Hermite polynomials.

polpak, a MATLAB code which evaluates a variety of mathematical functions.

stochastic_diffusion, a MATLAB code which implements several versions of a stochastic diffusivity coefficient.

subset, a MATLAB library which enumerates combinations, partitions, subsets, index sets, and other combinatorial objects.


  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. Roger Ghanem, Pol Spanos
    Stochastic Finite Elements: A Spectral Approach,
    Revised Edition,
    Dover, 2003,
    ISBN: 0486428184,
    LC: TA347.F5.G56.
  4. Dongbin Xiu,
    Numerical Methods for Stochastic Computations: A Spectral Method Approach,
    Princeton, 2010,
    ISBN13: 978-0-691-14212-8,
    LC: QA274.23.X58.
  5. 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 07 October 2022.