Subject: Calculus

# Gaussian Quadrature

Let us recall the more commonly used low-order Newton Cotes which is given by the formula

where x_i are the nodes and w_i are the weights. In deriving this formula, we did a two step process. First, we fixed the nodes { x_i } equidistantly and second, we find good weights { w_i } that minimizes the error by ensuring that the quadrature formula is interpolatory.

In Gauss quadrature, we would also want to optimize the location of the nodes so that we could minimize the error by finding good weights w_i and good nodes x_i. In this way, we can actually obtain a smaller error and at the same time giving us more degrees of freedom. However, if we could still recall from the previous section (Newton-Cotes Formula), we encounter a certain constraint as we go to higher order Newton-Cotes rule that would limit us to n nodes only. With the introduction of Gaussian quadrature, we get degrees of 2n – 1 which is an n-point quadrature formula that integrates all polynomials up to degree 2n – 1 exactly. This appears to be twice the degree that the Newton-Cotes formula.

Now, let us generalize the problem and consider weighted integrals (I_w(f)) of the type

Here, our w(x) is a positive integrable weight function. For simplicity, we will have the constraint at interval \left [ -1, 1\right ], which is no loss of generality. We approximate this weighted integral by the same type of interpolatory quadrature formula as before. Suppose the nodes -1 -1 \leq x_i \leq x_2 \leq \cdot \cdot \cdot \leq x_n \leq 1 and let p(x) be the unique interpolation polynomial in ℘_{n-1} for those nodes. Then we use

where

Now we introduce a theorem for us to be able to choose the right nodes,

### Theorem.

The formula Q_{n,w} (f) that is found in the previous equation is of the degree 2n – 1 if and only if

where

There is no set of nodes for which the degree is greater than 2n − 1. To back this up, let us provide a solid proof of the theorem.

### Proof.

Let f be an arbitrary polynomial in P_{2n−1}. We begin to show that

is a sufficient condition to conclude that Q_{n,w} (f)= I_w (f). Let q_{n−1} be the unique polynomial in ℘_{n-1} that interpolates f in the nodes {xj}. Then the nodes are roots of the polynomial f − q_{n−1} ∈ ℘_{2n-1} (among other roots) and we can write

For some polynomial π_{n-1} ∈ ℘_{n-1}$}. Thus, we can define

Since, Q_{n,w} is interpolatory, its degree d \geq n-1 and Q_{n,w} (q_{n-1})= I_w (q_{n-1} ). Moreover, by definition that w(x_i )=0 for all i. Therefore, if

holds, we will have

This proves that the equation is a sufficient condition. To show that it is also a necessary condition, let us suppose that there is a \tilde{p}(x) ∈ ℘_{n-1} for which the equation fails. Then, if f = ω_{n} \tilde{p}(x) ∈ ℘_{2n-1}, we get

but

This shows that the equation is necessary to be able to conclude that Q_{n,w} has degree 2n − 1, and the first part of the theorem is proved.

**BACK TO**: Numerical integration