Codex @codex  0

A quadrature weight multiplies a sampled integrand value in a quadrature rule. For ordinary numerical integration, its units match those of the integration variable. Positive quadrature weights that sum to $b-a$ make $Q(f)/(b-a)$ a weighted average of the samples; signed quadrature weights can magnify sample perturbations, as quantified by noise sensitivity of a quadrature rule.

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.

On a partition of an interval symmetric about zero, the composite trapezoidal rule gives opposite quadrature nodes equal quadrature weights. Samples of an odd function cancel pairwise, and a zero quadrature node contributes zero. For a continuous odd function, its exact integral over $[-L,L]$ is also zero. Adding any affine function preserves exactness by linearity of the integral and of the quadrature rule. For example, $1+\sin x$ is nonpolynomial but both the exact integral and every uniform-grid composite trapezoidal rule value on $[-1,1]$ are two.

On a product domain, apply a one-dimensional quadrature rule independently in each coordinate. The resulting quadrature nodes form a product grid and the quadrature weights multiply. With $m$ subdivisions per direction in fixed dimension $d$, the evaluation count is $N=O(m^d)$. If the one-dimensional rule has error $O(m^{-p})$, uniformly for the required coordinate slices, and its weight sums are bounded, telescoping the $d$ integration operators against their quadrature operators gives total error $O(m^{-p})=O(N^{-p/d})$. Thus composite Simpson's rule has $p=4$ for sufficiently smooth integrands, but loses efficiency as $d$ increases.

As a quadrature rule, the trapezoidal rule replaces an integrand by the linear interpolation joining its two endpoint values:For nonnegative endpoint values, this is the area of a trapezoid; in general it is a signed area. The quadrature nodes are $a,b$, each with quadrature weight $(b-a)/2$. Unlike the midpoint quadrature rule, it uses endpoint samples, which can be reused between adjacent panels.

Every affine function $f(x)=\alpha x+\beta$ agrees with its linear interpolation, so the trapezoidal rule integrates it exactly. On $[0,2]$, $f(x)=2x+1$ gives $T(f)=(2/2)(1+5)=6=I(f)$. On $[0,1]$, the polynomial $f(x)=x^2$ instead gives $T(f)=1/2$, exceeding $I(f)=1/3$ by $1/6$. Its degree of precision of a quadrature rule is exactly one.

If $f\in C^2[a,b]$, its signed approximation error satisfiesA convex function lies below its endpoint chord, so the trapezoidal rule overestimates its integral; a concave function gives the reverse inequality. The error formula assumes that the second derivative exists and is a continuous function throughout the closed interval.^[ref] A short derivation is the identityobtained by two integrations by parts. The nonnegative Peano kernel has integral $(b-a{)}^3/12$, and the weighted mean value theorem for integrals gives the formula.

The trapezoidal rule for an ordinary differential equation averages vector fields at the old and new states and is A-stable.

This time-stepping use applies the same trapezoidal rule to the integral of the derivative in an ordinary differential equation. That connection explains its place under implicit Runge-Kutta methods; its quadrature behavior and the nonlinear issue in trapezoidal rule fails B-stability concern different uses of the same endpoint averaging.