The Lagrange interpolation polynomial is
The distinct nodes make all denominators nonzero. Each Lagrange basis polynomial has degree and satisfies , so the sum has degree at most and interpolates every specified value.
If and both interpolate the data and have degree at most , the polynomial has the distinct roots . A nonzero polynomial of degree at most has at most roots: repeated use of the factor theorem would otherwise make a product of linear factors divide it. Hence , proving the interpolating polynomial is unique.
Integrate the Lagrange interpolation polynomial and exchange the integral with the finite sum:
Therefore the interpolatory quadrature weights are
The resulting quadrature rule is exact for every polynomial of degree at most , since its interpolant is itself. For a general continuous , replacing by its interpolant supplies the approximation.
The first pivot is two. Gaussian elimination without partial pivoting uses and , giving
The second pivot is one; gives the upper triangular matrix
Both pivots used for division are nonzero, so no row exchange is needed.
The elimination multipliers are , , and . The unit lower triangular matrix in the LU factorization is therefore
Direct multiplication checks
The positive signs in undo the row subtractions used during Gaussian elimination.
Take the invertible matrix , whose determinant is . Suppose it had an LU factorization with any lower triangular and upper triangular :
The entry gives , so . The entry gives , so . But then the entry would be , contradicting . Thus this nonsingular matrix has no unpermuted LU factorization, even without imposing a unit diagonal on . A row permutation resolves the zero pivot, explaining the role of partial pivoting.

Articles by others on the same topic (0)

There are currently no matching articles.