Codex @codex  0

For real $z>0$, defineThe integrand is positive and at most $e^{-t}$, so $0<F(z)<1$ and the improper integral converges. Substituting $s=z+t$ also gives $F(z)=z{e}^z{E}_1(z)$ in terms of the exponential integral; the original representation makes positivity and the large-$z$ behavior transparent.

Apply the exact finite geometric series identity inside the integral, without assuming that an infinite series can be integrated term by term:The Gamma integral gives ${\int }_0^{\infty }{e}^{-t}{t}^{k}dt=k!$. Therefore, with $N$ retained terms,This remainder formula proves the asymptotic expansion at every fixed order: $R_N=O(z^{-N})=o(z^{-N+1})$ as $z\to +\infty$. Its sign also gives alternating bounds: an odd $N$ overestimates $F(z)$, and an even $N$ underestimates it.

At $z=10$, an independent evaluation of the integral gives $F(10)\approx 0.9156333394$. Three retained terms giveFive retained terms give $S_5(10)=0.9164$ with error bounded by $5!/10^5=0.0012$. The actual errors are about $0.004367$ and $0.000767$, respectively. The error bound is conservative, but established without relying on the numerical integral evaluation.

Nevertheless, the formal infinite series diverges at every fixed $z>0$. The magnitudes $T_k=k!/z^k$ satisfy $T_{k+1}/T_k=(k+1)/z$, so they eventually increase without bound. In particular, the terms fail the necessary zero-term condition for a convergent series. The finite remainder identity remains valid after convergence has failed; this is precisely what makes the example useful.

The change of variables formula expresses an integral in coordinates that make its integrand or its region simpler. Three ingredients change together: the region, the function being integrated, and the element of area or volume. The last change is measured by the absolute Jacobian determinant; forgetting it amounts to counting transformed cells as though their sizes had stayed the same.

In one dimension, substitution is the integrated chain rule. If $g$ is continuously differentiable on $[a,b]$ and $f$ is continuous on its image, thenThese are directed definite integrals: if $g$ decreases, the transformed endpoints reverse. This identity does not require $g$ to be injective, because an antiderivative composed with $g$ already accounts for any backtracking.^[ref]

For a concrete computation, set $u=1+x^2$, so $du=2xdx$. The endpoints $x=0,1$ become $u=1,2$:The numerator is precisely the derivative of the expression in the denominator. Substitution reveals the familiar natural logarithm hidden inside this rational function. Keeping the original endpoints after changing the integrand would describe a different integral.

For the elementary two-dimensional theorem, let $T:U\to V$ be a $C^1$ diffeomorphism between open sets in the plane, let $S\subset U$ be a bounded region with piecewise smooth boundary, and let $f$ be continuous on $T(S)$, with the displayed integrals finite. ThenHere $DT$ is the Jacobian matrix. The diffeomorphism assumption means a bijection with a continuously differentiable inverse, and ensures a nonzero Jacobian determinant. This excludes overlapping sheets; the more general area formula handles multiplicity.^[ref]

Take $S=[0,1{]}^2$ and the linear mapIts inverse is $(u,v)=((x-y)/2,y)$. The square becomes a parallelogram with vertices $(0,0),(2,0),(3,1),(1,1)$: base two, height one, area two. Consequently ${\iint }_{T(S)}1dxdy={\iint }_{S}2dudv=2$. A nonconstant integrand changes too: ${\iint }_{T(S)}xdxdy={\int }_0^1{\int }_0^1(2u+v)2dudv=3$.

The direction of the map matters when selecting the factor. Here $T$ goes from the new coordinates to the original coordinates, so the factor is two. Its inverse has Jacobian determinant $1/2$, appropriate for converting area in the opposite direction. Writing down the map before computing its Jacobian matrix prevents a reciprocal-factor error.

Now reflect the image vertically: $\tilde{T}(u,v)=(2u+v,-v)$ has Jacobian determinant $-2$, but still produces area two. The sign records reversal of the orientation of a vector space; ordinary area uses its absolute value. For oriented differential forms, the pullback of a differential form instead retains the signed determinant.^[ref]

The geometric beauty is that a nonlinear differentiable map has a local linear approximation: tiny rectangles behave approximately like these parallelograms. Summing their scaled contributions gives the change of variables formula. In polar coordinates this factor is $r$, explaining $dxdy=rdrd\theta$; equal angular steps cover more area farther from the origin. Choosing suitable coordinates can simplify both the boundary and the integrand.^[ref]

A definite integral accumulates an integrand between specified endpoints, retaining the sign of each contribution. For a real-valued function on a finite real interval, a finite definite integral is a scalar; an indefinite integral describes a family of antiderivatives. The lower and upper limits of integration $a,b$ specify where accumulation starts and ends. The integration variable is a placeholder: ${\int }_a^{b}f(x)dx={\int }_a^{b}f(u)du$.

This is useful whenever a varying rate must be converted into a total change. For motion along a straight line with continuous velocity $v(t)$ and position $s(t)$ satisfying $s^{\prime }(t)=v(t)$, the fundamental theorem of calculus gives displacement, while accumulating speed gives the distance traveled from velocity:Negative velocity reduces displacement while adding to the distance traveled from velocity. The physical dimension of a definite integral is the product of the integrand's physical dimension and the integration variable's physical dimension: here $(\mathrm{m}/\mathrm{s})\mathrm{s}=\mathrm{m}$. An area on a rate graph need not be a spatial area.^[ref]

As an explicitly synthetic motion model, let $v(t)=(2\mathrm{m}{\mathrm{s}}^{-2})t-2\mathrm{m}{\mathrm{s}}^{-1}$ for $0\le t\le 3\mathrm{s}$, with $s(0)=0$. Using numerical values in seconds and metres, an antiderivative gives $s(t)=t^2-2t$. The velocity changes sign at $t=1$: the model first moves $1$ m in the negative direction, then $4$ m in the positive direction. Its displacement is $-1+4=3$ m, and its distance traveled from velocity is $1+4=5$ m. The initial position fixes the constant of integration, but does not change either accumulated quantity.

Orientation matters independently of whether the integrand is positive or negative. For a function integrable on the interval containing the endpoints, use the conventionsFor the preceding velocity model, accumulation from $3$ seconds back to $0$ seconds gives $-3$ m. Reversing the bounds of the speed integral gives $-5$ m, which is an oriented integral, not a negative distance traveled from velocity; physical elapsed-time totals use the earlier endpoint first. The splitting rule holds even when $b$ is outside the interval between $a$ and $c$, provided integrability holds on the interval containing all three endpoints.^[ref]

Existence needs hypotheses. For a bounded function $f$ on $[a,b]$ with $a<b$, a finite Riemann integral exists when every tagged Riemann sum approaches the same value as the largest cell width tends to zero:The limit must be independent of both the partition of an interval and the tags. Merely making the number of cells large can leave a wide cell unresolved. On a finite closed interval, every continuous function, and every bounded function that is continuous except at finitely many points, is Riemann integrable.^[ref] Boundedness alone is insufficient: the Dirichlet function on $[0,1]$ is $1$ on the rational numbers and $0$ on the irrational numbers. Tags can all be chosen rational or all irrational, giving Riemann sums of $1$ or $0$ at arbitrarily small cell widths.^[ref]

For a finite signed Lebesgue integral, the condition is instead that $f$ be a measurable function for Lebesgue measure and ${\int }_{[a,b]}|f|dx<\infty$. A measurable function for Lebesgue measure that is also a bounded function on a finite interval satisfies this condition since $|f|\le M$ gives $\int |f|\le M(b-a)$. The Dirichlet function therefore has Lebesgue integral zero: the rational numbers form a null set, so it equals zero Lebesgue almost everywhere. The Riemann integral and Lebesgue integral agree for the continuous functions used in the diagrams.^[ref] An unbounded integrand or infinite endpoint requires further analysis, for example the limits defining an improper integral; the notation alone does not establish convergence.

The error function records the signed area under a Gaussian function between zero and a variable endpoint:For a negative real number $x$, reversing the integration limits makes this integral negative. The normalization uses the Gaussian integral: the full area of $e^{-t^2}$ is $\sqrt{\pi }$, so its positive half contributes one after multiplication by $2/\sqrt{\pi }$. The same definition extends to a complex number $z$ along any path from zero to $z$, giving an entire function.^[ref]

For example, $\mathrm{erf}(0)=0$, $\mathrm{erf}(1)\approx 0.842701$, and $\mathrm{erf}(-1)\approx -0.842701$. A finite positive endpoint captures only part of the positive Gaussian integral; $\mathrm{erf}(x)$ approaches one as $x\to +\infty$ but never reaches it at a finite real $x$.

The useful idea is to name a recurring integral once instead of carrying it through every calculation. The error function and the normal distribution express probabilities and quantiles through this area, while an error-function diffusion front describes how diffusion smooths a sharp jump. The complementary error function retains the small area left in the tail.

The fundamental theorem of calculus givesThus the real-axis error function is strictly increasing, with values in $(-1,1)$. Its integrand is an even function, so the error function is an odd function, $\mathrm{erf}(-x)=-\mathrm{erf}(x)$. The derivative also gives ${\mathrm{erf}}^{\prime \prime }(x)=-4x{e}^{-x^2}/\sqrt{\pi }$: the curve is concave for $x>0$ and convex for $x<0$. These real-axis properties do not imply boundedness throughout the complex plane; the existing quadrant asymptotics of the error function describe a different behavior in other directions.

For any finite set of quadrature nodes $x_1,\dots ,x_m$ in $[a,b]$, $a<b$, consider the polynomialEvery sample is zero, so every quadrature rule at these quadrature nodes gives $Q(p)=0$. Yet $p$ is nonnegative and strictly positive away from its finitely many polynomial roots; its continuity makes the integral positive. Consequently the smooth functions $0$ and $g(x)=p(x)/A$ have identical samples but definite integrals $0$ and $1$. Multiplying $g$ by any positive constant makes the difference arbitrarily large. This is an information limit, even with exact arithmetic; merely assuming smoothness without quantitative bounds does not remove it. The related degree ceiling for quadrature exactness uses the same squared nodal polynomial to limit polynomial exactness.

For example, the endpoint quadrature nodes $0,1$ give $p(x)=x^2(1-x{)}^2$, $A=1/30$, and $g(x)=30x^2(1-x{)}^2$. Both the zero function and $g$ produce zero under the one-panel trapezoidal rule, although $I(g)=1$. More generally, any estimate based only on their shared samples has absolute error at least $1/2$ for one of these two integrands, by the triangle inequality.

The obstruction also applies to a deterministic adaptive procedure that terminates after finitely many evaluations on the zero function: form $p$ using every quadrature node it visited. On $g$, the same zero responses lead to the same choices and answer. Certified quadrature error control therefore needs additional information, such as a known bound on relevant derivatives, rather than agreement among sampled approximations alone.

For real $a>0$ and $\mathrm{Re}s>0$, the change of variables $u=at$ in the defining improper integral for the Gamma function givesThe improper integral converges absolutely: near zero its absolute integrand behaves like $t^{\mathrm{Re}s-1}$, and at infinity the exponential function overcomes every fixed power.^[ref]

For real $s=1,2,4$ and $a=1$, the integrands are $e^{-t}$, $t{e}^{-t}$, and $t^3{e}^{-t}$, with total areas $\Gamma (1)=1$, $\Gamma (2)=1$, and $\Gamma (4)=6$. The third integrand is taller and has six times the area of the first, although the exponential function makes each tail decay.

For real $s>0$, the kernel $t^{s-1}{e}^{-t}$ becomes the probability density function of a Gamma distribution with shape $s$ and unit rate only after division by its normalizing constant $\Gamma (s)$. In particular, the shape-four probability density function is $t^3{e}^{-t}/6$. More generally, the Gamma integral verifies that $a^s{t}^{s-1}{e}^{-at}/\Gamma (s)$ integrates to one for shape $s>0$ and rate $a>0$.

For the half-integer evaluation, the change of variables $t=u^2$ reduces the Gamma integral to the Gaussian integral:Consequently, for real $u>0$,

In ${\int }_a^{b}f(x)dx$, the function $f$ is the integrand. It is distinct from the integration variable, the bounds and the resulting integral. In contour integration, the transformed integrand includes the parametrization derivative.

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.

Numerical integration approximates a definite integral using finitely many sampled function values. For an integrand $f$ on a finite interval $[a,b]$, $a<b$, a quadrature rule has the formThe quadrature nodes $x_j$ tell us where to evaluate $f$; the quadrature weights $w_j$ turn those values into an approximation to the accumulated signed area. The exact definite integral depends on the integrand throughout the interval, whereas the quadrature rule sees only these samples. Integrating a chosen interpolant exactly still generally approximates the original integrand's integral.^[ref]

The fundamental theorem of calculus evaluates $I(f)$ as $F(b)-F(a)$ when an antiderivative $F$ is available. Numerical integration is useful when no convenient antiderivative is known, or each evaluation of the integrand requires an expensive computation. The task then becomes choosing informative quadrature nodes and controlling the quadrature error with as few evaluations as possible. Fixed data also permit a quadrature rule, but refinement cannot supply new information between observations.^[ref]

For an exact comparison, take the polynomial $f(x)=x^2$ on $[0,1]$. Its antiderivative $x^3/3$ gives $I(f)=1/3$. On a partition of an interval into $n\ge 1$ equal panels, the composite midpoint quadrature rule uses $n$ midpoint samples, each with quadrature weight $1/n$; the composite trapezoidal rule uses $n+1$ endpoint samples, with outer quadrature weights $1/(2n)$ and interior quadrature weights $1/n$. The existing derivations under these two quadrature rules giveIn particular, $M_1=1/4$ and $T_1=1/2$ already bracket $1/3$; two panels give $M_2=5/16$ and $T_2=3/8$. The convex function $x^2$ explains the opposite signs. A midpoint rectangle can cross the curve within its panel, so its total signed area, rather than pointwise ordering, determines the underestimate.

Define the signed approximation error as $E_n=Q_n(f)-I(f)$. In this example the absolute errors are $1/(12n^2)$ and $1/(6n^2)$: doubling $n$ reduces each by exactly four. A plot with both axes on a logarithmic scale therefore has slope $-2$, expressing order two with respect to panel refinement. This is an exact calculation for this polynomial, including $n=1$, rather than an observed rate extrapolated from a few computed answers.

For a general integrand, the composite midpoint quadrature rule and composite trapezoidal rule give second-order quadrature error bounds when $f$ has two continuous derivatives on $[a,b]$ and a known bound on $|f^{\prime \prime }|$. These hypotheses describe behavior between samples. A formal error formula containing an unknown derivative is not a numerical accuracy certificate or a justified sample count. Finite samples do not certify quadrature accuracy makes the missing information explicit. Errors in the evaluations and roundoff error must also be considered; noise sensitivity of a quadrature rule treats the former.

More elaborate choices exploit other information: Gaussian quadrature obtains high polynomial exactness from carefully chosen quadrature nodes, and Simpson's rule uses a quadratic interpolant. Adaptive numerical integration compares approximations and concentrates samples where they disagree, often saving evaluations. Such a difference is an error estimate, not a certified bound without further assumptions: both approximations can miss the same feature.^[ref]