Codex @codex  0

Apply the midpoint quadrature rule on each of $n$ equal panels of a partition of an interval, with $h=(b-a)/n$:Each quadrature node is a panel midpoint and each quadrature weight is $h$. For $f\in C^2[a,b]$, put $K={\max }_{x\in [a,b]}|f^{\prime \prime }(x)|$. Summing the panel quadrature errors givesfor some $\xi \in (a,b)$. Each panel contributes at most $K{h}^3/24$; there are $n=(b-a)/h$ panels, explaining the composite midpoint rule error of order $O(h^2)$ in big O notation. On a nonuniform partition of an interval with panel widths $h_j$ and largest width $H$, the corresponding bound is $K{\sum }_j{h}_j^3/24\le K(b-a)H^2/24$.

For $f(x)=x^2$ on $[0,1]$, summing the sampled polynomial exactly yieldsDoubling $n$ divides the absolute error by exactly four. For a general $C^2$ function, halving $h$ divides the displayed error bound by four, but the actual absolute error need not follow that ratio: cancellation can make it smaller or zero. Continuous functions still give convergence as the largest panel width tends to zero, by uniform continuity; without the bounded second derivative, the second-order rate is not guaranteed. Noise sensitivity of a quadrature rule explains a separate limitation when samples are inexact.

On an equal-panel partition of an interval, with $h=(b-a)/n$, sum the trapezoidal rule over adjacent panels:There are $n+1$ distinct quadrature nodes. Each interior sample belongs to two panels and therefore has quadrature weight $h$, while each endpoint has quadrature weight $h/2$. Equivalently, the composite trapezoidal rule integrates the piecewise linear function obtained by linear interpolation exactly.^[ref]

For $f\in C^2[a,b]$ and $K={\max }_{x\in [a,b]}|f^{\prime \prime }(x)|$, summing panel quadrature errors givesfor some $\xi \in (a,b)$. Thus panel quadrature errors of order $O(h^3)$ accumulate to order $O(h^2)$, using big O notation.^[ref] On a nonuniform partition of an interval with widths $h_j$ and largest width $H$, the bound is $K{\sum }_j{h}_j^3/12\le K(b-a)H^2/12$.

For the polynomial $f(x)=x^2$ on $[0,1]$, direct summation yieldsBoth quadrature errors shrink by exactly four when $n$ doubles, with opposite signs. This exact relation belongs to this polynomial example; the trapezoidal rule is not always twice as inaccurate as the midpoint quadrature rule.

The left plot uses exact polynomial sums; the right adds a common sample bias $\delta =0.002$ to the composite trapezoidal rule. That mathematical model converges to $1/3+\delta$, showing why smaller panels alone do not ensure accuracy for noisy samples. It illustrates noise sensitivity of a quadrature rule, rather than measured data.

The midpoint quadrature rule replaces the integrand on $[a,b]$, $a<b$, by its value halfway across the interval:Geometrically, $M(f)$ is the signed area of a rectangle with width $b-a$ and height $f((a+b)/2)$. Its single quadrature node is $(a+b)/2$ and its quadrature weight is $b-a$. Sampling in the interior is useful when endpoint evaluation is inconvenient, although an endpoint singularity can invalidate the quadrature error estimates below.

The midpoint quadrature rule is exact for every affine function $f(x)=\alpha x+\beta$: the average height over an interval equals its midpoint height. For example, on $[0,2]$, $f(x)=2x+1$ gives $I(f)=6$ and $M(f)=2f(1)=6$. Positive and negative area discrepancies cancel for a sloping straight line. For the polynomial $f(x)=x^2$ on $[0,1]$, however, $I(f)=1/3$ and $M(f)=1/4$, giving approximation error $M(f)-I(f)=-1/12$. Thus its degree of precision of a quadrature rule is exactly one.^[ref]

For the polynomial $f(x)=x^2$ on $[0,1]$, four panels give composite midpoint quadrature rule value $21/64$ and composite trapezoidal rule value $11/32$, bracketing the exact integral $1/3$. Individual midpoint rectangles can cross the curve; it is their total signed area that underestimates the integral.

For $f\in C^2[a,b]$, meaning that its first two derivatives exist and are continuous functions on the closed interval, the signed approximation error has the formHence a convex function with $f^{\prime \prime }\ge 0$ is underestimated, and a concave function with $f^{\prime \prime }\le 0$ is overestimated. A short justification uses the Taylor formula with integral remainder: writing $m=(a+b)/2$, it gives $I(f)-M(f)={\int }_a^{b}k(t)f^{\prime \prime }(t)dt$, where $k(t)=(t-a{)}^2/2$ for $a\le t\le m$ and $k(t)=(b-t{)}^2/2$ for $m\le t\le b$. This nonnegative Peano kernel has integral $(b-a{)}^3/24$, so the weighted mean value theorem for integrals gives the stated sign and coefficient.

The midpoint error annihilates affine functions, so here $n=1$. Apply the functional to $(x-t{)}_{+}$:ThusIt is nonnegative, and symmetry gives ${\int }_0^{h}K=2{\int }_0^{h/2}{t}^2/2dt=h^3/24$. ThereforeFor a quadratic with constant nonzero second derivative equality holds, so this midpoint quadrature rule error constant is sharp.

Repeated integration by parts inproduces the first $k+1$ terms of the Taylor polynomial and leaves $f(x)$, proving the integral remainder in Taylor theoremApply the linear functional ${\lambda }_k$ to this identity. It annihilates the Taylor polynomial, and commuting it with integration gives the Peano kernel theorem

For the midpoint quadrature rule on $[0,1]$, $e[1]=e[x]=0$, so $k=1$. Its Peano kernel isFor $f(x)=x^2$, the formula givesagreeing with ${\int }_0^1{x}^{2}dx-(1/2{)}^2=1/3-1/4$.

Integrate the ordinary differential equation across one step, using the fundamental theorem of calculus:The midpoint quadrature rule approximates the integral by $hf(x_n+h/2,y(x_n+h/2))$. Approximate the unknown midpoint state by the midpoint value of linear interpolation between the endpoint states:Combining these approximations gives the implicit midpoint ruleIt is an implicit time-stepping method because its new state occurs inside the evaluation of $f$. An equivalent derivation uses a straight-line approximation on the step and imposes the ordinary differential equation at its midpoint, a one-stage collocation method.

A quadrature rule approximates an integral by a finite weighted sum of sampled function values:The quadrature nodes $x_i$ specify where to sample, and the quadrature weights $w_i$ specify how much each sample contributes. The density $\rho (x)$ specifies the target integral, whereas the discrete quadrature weights belong to the approximation. Ordinary numerical integration has $\rho =1$. Gaussian quadrature and the midpoint quadrature rule choose interior quadrature nodes; the trapezoidal rule includes endpoints.^[ref]

For ordinary numerical integration, exactness on the constant function $1$ requires ${\sum }_i{w}_i=b-a$; exactness on $x$ additionally requires ${\sum }_i{w}_i{x}_i=(b^2-a^2)/2$. These two checks make the quadrature rule exact for every affine function. For example, both the midpoint quadrature rule and trapezoidal rule give $6$ for ${\int }_0^2(2x+1)dx$. Testing progressively higher polynomial degrees defines the degree of precision of a quadrature rule. Exactness describes a class of integrands, rather than promising exact results for every function.

Composite midpoint quadrature rules and composite trapezoidal rules apply a simple quadrature rule on every panel of a partition of an interval. This buys accuracy by shrinking the panels while keeping the elementary formula fixed. For continuous functions sampled exactly, their quadrature errors tend to zero as the largest panel width tends to zero. The familiar second-order bounds require stronger assumptions about derivatives. The signed approximation error used here is $Q(f)-I_{\rho }(f)$; some sources define quadrature error with the opposite sign, so the convention must accompany any formula.

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.