Gaussian Quadrature

Gaussian quadrature is a numerical integration rule that chooses abscissas and weights to integrate polynomials of the highest possible polynomial degree exactly. The fundamental theorem of Gaussian quadrature states that the optimal abscissas of the -point Gaussian quadrature formulas are precisely the roots of the orthogonal polynomials of polynomial degree for the same interval and weighting function. For a nonnegative weighting function with finite moments of all orders and positive integral against the square of every nonzero polynomial, Gaussian quadrature fits all polynomials up to polynomial degree exactly. Slightly less optimal fits are obtained when nodes are constrained to include endpoints, as in Radau quadrature or Lobatto quadrature.

	interval	are roots of
1

To determine the weights corresponding to the Gaussian abscissas , compute a Lagrange interpolating polynomial for by letting

(1)

(where Chandrasekhar 1967 uses instead of ), so

(2)

Then fitting a Lagrange interpolating polynomial through the points gives

(3)

for arbitrary points . We are therefore looking for a set of points and weights such that for a weighting function ,

			(4)
			(5)

with weight

(6)

The weights are sometimes also called the Christoffel numbers (Chandrasekhar 1967). For the member of polynomial degree in a family of orthogonal polynomials, set above. Then

(7)

where is its leading coefficient. The weights are given by

			(8)
			(9)

(Hildebrand 1956, p. 322), where

(10)

Using the relationship

(11)

(Hildebrand 1956, p. 323) gives

(12)

(Note that Press et al. 1992 omit the factor .) In Gaussian quadrature, the weights are all positive. The error is given by

			(13)
			(14)

where (Hildebrand 1956, pp. 320-321).

Other curious identities are

(15)

and

			(16)
			(17)

(Hildebrand 1956, p. 323).

In the notation of Szegö (1975), let be an ordered set of points in , and let , ..., be a set of real numbers. If is an arbitrary function on the closed interval , write the Gaussian quadrature as

(18)

Here are the abscissas and are the Cotes numbers.

The nodes and weights can also be obtained from the leading block of the Jacobi matrix for the normalized orthogonal polynomials. Its eigenvalues are the nodes . If is a unit eigenvector for and , then

(19)

This is the method of Golub and Welsch (1969, pp. 221-223). Equivalently, the Christoffel function gives .

For example, the three-point Legendre-Gauss quadrature rule has

(20)

Its nodes are , 0, and , with corresponding weights , , and . The rule is exact through polynomial degree 5.

The rule's higher-degree exactness is illustrated above by the quartic polynomial and its quadratic Lagrange interpolating polynomial ,

			(21)
			(22)

Their values agree at the three nodes. The blue and orange shaded regions lie where is positive and negative, respectively. Their signed integrals cancel. Consequently,

(23)

The numbers beneath the nodes are the corresponding weights.