Codex @codex  0

Polynomial interpolation makes Simpson's rule exact for every polynomial of polynomial degree at most two. To understand the extra polynomial degree, centre the coordinate at the midpoint: $y=x-m$, with endpoints $-h,h$. For any cubic polynomial $p$, let $p_2$ be its interpolating quadratic polynomial. The polynomial $p-p_2$ vanishes at all three quadrature nodes, sowhere $A$ is the coefficient of $x^3$ in $p$. As a function of $y$, this is an odd function, whose integral over the symmetric interval $[-h,h]$ is zero. Hence integrating $p_2$ also integrates $p$ exactly. This is the symmetry mechanism in the nodal-polynomial criterion for quadrature exactness.

In the worked example, the polynomial difference factors asIts definite integrals on $[-1,1/2]$ and $[1/2,2]$ are respectively $81/64$ and $-81/64$. Their cancellation concerns signed integrals, rather than pointwise accuracy of polynomial interpolation.

The degree of precision of a quadrature rule stops at three: on the same interval, the monomial $x^4$ gives $S(x^4)=69/8$, whereas ${\int }_{-1}^2{x}^{4}dx=33/5$. Their difference, exact integral minus approximation, is $-81/40$.

A cubic polynomial in one variable has polynomial degree three. Over a splitting field, its factorization of polynomials contains three polynomial factors of polynomial degree one when multiplicity of a root is counted.

Move the constant term to the left and factor the cubic polynomial:Setting each factor to zero gives all three polynomial roots:

Write the remaining quadratic polynomial as $x^2+bx+c$. Comparing coefficients in $(x-2)(x^2+bx+c)$ with the given cubic polynomial gives $b-2=-1$ and $-2c=-2$, hence $b=c=1$. The coefficient of $x$ is then $c-2b=-1$, as required. Thus

Let $p(x)=x^3-5x-1$. Its values at useful points areSince a polynomial is a continuous function, the intermediate value theorem gives one polynomial root in each of the disjoint intervals $(-3,-2)$, $(-1,0)$ and $(2,3)$. A nonzero cubic polynomial has at most three distinct polynomial roots. Hence it hasThese three distinct polynomial roots already exhaust its degree of a polynomial, so each has multiplicity of a root one.

Simpson's rule approximates a definite integral by integrating the quadratic polynomial obtained by polynomial interpolation at two endpoints and their midpoint. For an interval $[a,b]$ with $a<b$, write $m=(a+b)/2$ and $h=(b-a)/2$. The three quadrature nodes are $a,m,b$, and the quadrature rule isThus its quadrature weights are $h/3,4h/3,h/3$. Their sum is $2h=b-a$, as required for exactness on constant functions. This is the three-node Newton-Cotes closed quadrature; its degree of precision of a quadrature rule is three, even though the interpolating polynomial has polynomial degree at most two.^[ref]

For an original worked example, take the cubic polynomial $f(x)=x^3-2x+4$ on the interval $[-1,2]$. Here $m=1/2$ and $h=3/2$, so the function evaluations giveSimpson's rule givesAn antiderivative supplies an independent exact definite integral:The quadratic polynomial through the same quadrature nodes is $p_2(x)=\frac{3}{2}{x}^2-\frac{1}{2}x+3$. It differs from $f$ between the quadrature nodes, but its definite integral agrees exactly. The picture shows why: the signed polynomial interpolation error has two lobes with opposite integrals.

Three function evaluations therefore capture more information about an integral than about the whole function. When an antiderivative is unavailable, composite Simpson's rule repeats this inexpensive local construction over a partition of an interval. For a sufficiently smooth function, halving the node spacing reduces its truncation error by approximately a factor of sixteen.

Take the stipulated polynomial roots $-2,1,5$. Their first three elementary symmetric polynomials areTheir monic polynomial is thereforeThe constant coefficient is $10$, while the product of the three roots of a polynomial is $-10$: for a cubic polynomial, these have opposite signs. Expanding the factors directly verifies the same polynomial identity.

For the nonmonic polynomial $-3x^3+12x^2+21x-30$, the same ratios give sum $-12/(-3)=4$, pairwise-product sum $21/(-3)=-7$, and product $-(-30)/(-3)=-10$. It is $-3$ times the monic polynomial above, so its roots of a polynomial are the same. Checking only the sum and product would leave the middle elementary symmetric polynomial unchecked; all three are needed to recover a monic polynomial of polynomial degree three.