quadrature_golub_welsch

quadrature_golub_welsch, a Fortran90 code which illustrates the use of the Golub Welsch technique for computing the points and weights of a Gaussian quadrature rule.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

quadrature_golub_welsch is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version.

Related Data and Programs:

quadrature_golub_welsch_test

clenshaw_curtis_rule, a Fortran90 code which defines a Clenshaw Curtis quadrature rule.

quadrature_least_squares, a Fortran90 code which computes weights for "sub-interpolatory" quadrature rules, that is, it estimates integrals by integrating a polynomial that approximates the function data in a least squares sense.

quad_rule, a Fortran90 code which defines 1-dimensional quadrature rules.

toms655, a Fortran90 code which computes the weights for interpolatory quadrature rule;
this library is commonly called IQPACK;
this is a Fortran90 version of ACM TOMS algorithm 655.

Reference:

  1. Philip Davis, Philip Rabinowitz,
    Methods of Numerical Integration,
    Second Edition,
    Dover, 2007,
    ISBN: 0486453391,
    LC: QA299.3.D28.
  2. Sylvan Elhay, Jaroslav Kautsky,
    Algorithm 655: IQPACK, Fortran Subroutines for the Weights of Interpolatory Quadrature,
    ACM Transactions on Mathematical Software,
    Volume 13, Number 4, December 1987, pages 399-415.
  3. Jaroslav Kautsky, Sylvan Elhay,
    Calculation of the Weights of Interpolatory Quadratures,
    Numerische Mathematik,
    Volume 40, 1982, pages 407-422.

Source Code:

Last revised on 22 August 2020.