Codex @codex  0

An improper integral whose absolute-value integral is finite. The triangle inequality for an integral bounds each truncated tail by the corresponding tail of $|f|$, so the Cauchy criterion for an improper integral gives ordinary convergence. Absolute convergence is stronger than a Cauchy principal value, which may exist only because two divergent sides cancel. The continuous sine-product integrand $\sin (ax)\sin (bx)/x^2$ has a removable value $ab$ at zero and an integrable $1/x^2$ tail; its auxiliary complex contour integrand can nevertheless require a principal value.

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.

Choose $\Theta$ with uniform distribution on $(0,\pi /2)$ and set $X={\sin }^2\Theta$. Since the sine is increasing on this interval,Differentiation gives $f(x)=1/(\pi \sqrt{x(1-x)})$, and $F(0)=0$, $F(1)=1$ prove normalization. This also evaluates the density's improper integral without treating its divergent endpoint values as infinite probability.

For the circle construction, instead choose $\Phi$ with uniform distribution on $[0,2\pi )$ and put $Y=\cos \Phi$. Each $-1<y<1$ has two preimages in one revolution. The change of variables formula adds both contributions:The affine change of variables $X=(1+Y)/2$ gives the standard arcsine distribution. A single monotone-branch calculation without the second preimage would give only half the total probability.

The same geometry explains the ideal harmonic oscillator $y(t)=A\cos (\omega t+\varphi )$, with amplitude $A>0$ and angular frequency $\omega >0$. Choose $T$ uniformly in any complete period of duration $2\pi /\omega$. Its oscillation phase modulo $2\pi$ is uniform, and $Y=y(T)$ has densityNear a turning point, $|dy/dt|=\omega \sqrt{A^2-y^2}$ is small, so a fixed position interval takes more time to cross. For the chosen amplitude $A=2\mathrm{cm}$, exactly one third of the period is spent in $|y|\le 1\mathrm{cm}$, and two thirds outside that central interval. These are calculations for an ideal sinusoid, rather than measured oscillator data.

The sampling assumption matters. Taking one observation every full period returns the same position each time. It does not produce this probability distribution. Uniformly random observation times and a fixed regularly spaced sampling scheme are different models; aliasing can conceal the motion altogether.

Paul Lévy used a related geometric construction in his work on stochastic processes: project a uniformly chosen point on a circle onto its diameter to divide a segment. The projection description appears explicitly in §5 of his paper, following its arcsine density calculation.^[ref]

The beta distribution is an absolutely continuous probability distribution for a random variable $X$ between zero and one. It can describe uncertainty about a probability, such as an unknown coin's chance of heads. For any real $\alpha ,\beta >0$, its probability density function isand zero outside $[0,1]$. The normalizer is the Euler beta function, related to the gamma function by the beta-gamma identity:Thus the probability density function integrates to one. ^[ref].

The endpoint behavior explains why the restrictions are exactly $\alpha ,\beta >0$: near zero the probability density function behaves like a constant times $x^{\alpha -1}$, and near one like a constant times $(1-x{)}^{\beta -1}$. These improper integrals converge for positive parameters, even when the probability density function diverges at zero for $\alpha <1$ or at one for $\beta <1$. There is no probability mass at either endpoint: $\mathbb{P}(X=0)=\mathbb{P}(X=1)=0$. A probability density function's height is not a probability; an interval's probability is the area under it.

Some useful examples of the beta distribution are:

The arcsine distribution curve is plotted only for $0.001\le x\le 0.999$; its finite plotted ends do not represent finite endpoint values of the probability density function. Using the inverse sine, the shaded central interval has exact probability $\mathbb{P}(0.1<X<0.9)=1-4\arcsin (\sqrt{0.1})/\pi \approx 0.5903$, calculated analytically rather than by treating the cropped curve as normalized.

To separate location from spread, write $\kappa =\alpha +\beta$ and $m=\alpha /\kappa$. The mean and variance of the beta distribution areThe mean specifies the balance between the parameters; their sum $\kappa$ controls concentration about that mean. Multiplying both parameters by the same positive number preserves the mean and changes the variance. Larger $\kappa$ puts more area near $m$ and makes the probability density function taller there; the total area stays one. Small $\kappa$ allows much greater spread, sometimes favoring the endpoints rather than the mean itself.

For example, keeping mean $m=0.3$ and taking $\kappa =5,20,80$ gives the beta distributions $\mathrm{Beta}(1.5,3.5)$, $\mathrm{Beta}(6,14)$, and $\mathrm{Beta}(24,56)$. Their variances are $0.035$, $0.01$, and $0.21/81\approx 0.00259259$, respectively. Parameters need not be integers.

The formulas follow from two short moment calculations. Multiplying the probability density function by $x$ or $x^2$ increases the first parameter of the Euler beta function:Apply the Gamma function recurrence to the beta-gamma identity, then subtract $(\mathbb{E}[X]{)}^2$ from $\mathbb{E}[X^2]$ to obtain the variance.

The beta distribution is especially convenient as a prior distribution for the unknown success probability $P$ of Bernoulli trials. The Beta-binomial conjugacy says that $h$ successes and $t$ failures update a prior distribution $\mathrm{Beta}(\alpha ,\beta )$ to a posterior distribution $\mathrm{Beta}(\alpha +h,\beta +t)$. For example, suppose a coin has prior distribution $\mathrm{Beta}(2,2)$ and we observe eight heads and two tails. Its posterior distribution is then $\mathrm{Beta}(10,4)$, with posterior mean $5/7$. ^[ref].

For a locally integrable real or complex function, the displayed tail condition is equivalent to ordinary convergence of its improper integral at infinity. Apply the Cauchy criterion for a function limit to its truncated integrals, using completeness of the real or complex numbers. The necessity follows by subtracting two nearby-to-limit truncations; sufficiency follows by first taking a convergent sequence of truncations and then bounding all other late truncations against it. The analogous criterion applies at a finite singular endpoint. Bounding the modulus of every tail by its absolute-value integral proves that an absolutely convergent improper integral converges ordinarily.

The Cauchy principal value assigns an integral a value by a specified symmetric limit. Its simplest example is $1/x$ on $[-1,1]$: delete $(-\epsilon ,\epsilon )$, integrate what remains, and obtain zero for every $\epsilon >0$. The positive signed area under a graph cancels the negative one, although each grows without bound as the gap closes. This is a useful cancellation rule, but it does not make the ordinary improper integral exist.

Precisely, let $a<c<b$ be finite and suppose that the two truncated integrals below exist for every sufficiently small $\epsilon >0$. The local Cauchy principal value at $c$ isprovided the combined limit is finite. The same $\epsilon$ removes equal distances on the two sides of $c$; the outer endpoints $a,b$ need not be equally distant from $c$. By contrast, the ordinary improper integral requires both one-sided limits to be finite separately before adding them. Whenever those separate limits exist, the Cauchy principal value agrees with their sum.^[ref]

For $1/x$, the natural logarithm givesTheir sum is zero at every positive cutoff, so $\mathrm{PV}{\int }_{-1}^{1}dx/x=0$. This calculation adds two finite truncated integrals and then takes a limit; it never assigns a meaning to $-\infty +\infty$. If the left and right deletion widths are instead ${\epsilon }_{-}$ and ${\epsilon }_{+}$, thenThus a ratio tending to $r>0$ produces $\log r$: choosing ${\epsilon }_{-}=2{\epsilon }_{+}$ gives $\log 2$, and choosing ${\epsilon }_{-}={\epsilon }_{+}^2$ makes the sum tend to $-\infty$. Both widths can tend to zero without selecting a unique answer. The local Cauchy principal value selects the ratio $1$.

Unequal outer endpoints leave a finite contribution even with symmetric local deletion. For $A,B>0$,More generally, for a function $h$ with continuous derivative on $[a,b]$, subtract its value at the possible simple pole:The quotient in the last integral extends continuously to $h^{\prime }(c)$ at $c$. This formula explains both the cancellation and a practical way to compute it: integrate a continuous remainder instead of subtracting two large, nearly equal divergent contributions.

A finite Cauchy principal value is weaker than absolute convergence. Indeed,Taking the absolute value destroys the cancellation. Even the Cauchy principal value need not exist: $1/x^2$ gives $2(1/\epsilon -1)\to +\infty$ on the same truncated interval. Absolute integrability makes the cutoff prescription unnecessary; an ordinary improper integral can also converge without absolute convergence, as the Dirichlet integral illustrates. The three notions should therefore be checked separately.

There is an independent Cauchy principal value at infinity. For a function integrable on every finite interval, definewhen the limit is finite. The ordinary improper integral instead requires convergence of the two tails separately.^[ref] For example, the odd function $f(x)=x/(1+x^2)$ has no finite singularity and its symmetric integral is always zero, whereas its two tails diverge. Unequal truncations giveLetting $B/A\to r>0$ as $A,B\to \infty$ yields $\log r$. Symmetry at infinity is therefore a separate choice from symmetry near a finite pole. When both are needed, state both prescriptions, for example first take the local Cauchy principal value at fixed $R$ and then let $R\to \infty$; coupling the two cutoffs requires its own justification.

The local cancellation also defines the principal-value reciprocal distribution. Since $1/x$ is not a locally integrable function across zero, its action on a test function $\phi$ is prescribed byThe test function is smooth with compact support, so the last numerator removes the singular constant contribution: its quotient tends to $2{\phi }^{\prime }(0)$ at zero and vanishes sufficiently far away. This defines a distribution, not a finite pointwise value of $1/x$ at zero. It extends to Schwartz functions as a tempered distribution, and is also the distributional derivative of the logarithmic modulus.^[ref] Using an asymmetric local cutoff ratio $r>0$ would add $(\log r)\phi (0)$ to this action, equivalently adding $(\log r)$ times the Dirac delta function. Thus the cutoff choice matters even after passing to distributions.

Why keep this prescription? A Fourier transform naturally encounters such a kernel. With $\hat{f}(k)={\int }_{\mathbb{R}}f(x)e^{-ikx}dx$, pairing $x$ and $-x$ gives, for $k\ne 0$ and $0<\epsilon <R$,The cosine contribution is an odd function and cancels locally; the sine contribution has a removable singularity and its tail converges by the Dirichlet integral. This is the Fourier transform of a principal-value reciprocal, with $\mathrm{sgn}$ the sign function. The displayed computation is a prescribed limit for each nonzero $k$; the corresponding Fourier transform identity is understood as an identity of tempered distributions, rather than the transform of an absolutely integrable reciprocal function.

The Hilbert transform uses exactly this principal-value reciprocal distribution as a convolution kernel. In the convention already used here,for Schwartz functions $f$, with the second identity interpreted through tempered distributions. Hence the singular integral has a simple description as a Fourier multiplier.^[ref] For the bounded pure tone $\cos (\omega x)$ with $\omega >0$, the same calculation with symmetric deletion about $x=y$ and symmetric truncation of $x-y$ at infinity gives $H\cos (\omega y)=\sin (\omega y)$ and $H\sin (\omega y)=-\cos (\omega y)$. These tones lie outside the Schwartz space; their prescribed Cauchy principal values explain the quarter-cycle shift in the figure. Both figures are calculated mathematical illustrations. Reversing the Hilbert denominator reverses these signs; NIST uses that opposite convention.^[ref]

For strain energy density $W=\mu [({\lambda }_r^2+{\lambda }_{\theta }^2+{\lambda }_{\varphi }^2{)}^{\alpha }-3^{\alpha }]/2$ with $\mu >0$ and $\alpha >0$, the critical loading for elastic cavitation is finite exactly in the displayed range. To prove this, put $S={\lambda }^{-4}+2{\lambda }^2$ in the reduced strain energy density. The critical-pressure integrand is $2\mu \alpha S^{\alpha -1}({\lambda }^{-2}+{\lambda }^{-5})$. It is finite at $\lambda =1$ and asymptotic to $\mu \alpha 2^{\alpha }{\lambda }^{2\alpha -4}$ at infinity. The comparison test for the resulting improper integral gives $2\alpha -4<-1$.

For a spherical vacuum cavity initially at rest with radius $a$, in an inviscid liquid of constant mass density $\rho$ and positive constant far-field pressure $p_{\infty }$, the Rayleigh equation for an inviscid spherical bubble has the first integral $R^3{\dot{R}}^2=2p_{\infty }(a^3-R^3)/(3\rho )$. To obtain it, multiply the radius equation by $2R^2\dot{R}$ and integrate the resulting exact derivative. Choosing the negative square root and integrating $dt=dR/\dot{R}$ givesThis improper integral converges at both endpoints, proving finite-time Rayleigh collapse of a spherical cavity in the ideal model. It has numerical value approximately $0.914681$. The prediction neglects viscosity, surface tension and liquid compressibility.

If measurable functions satisfy $|f|\le g$ almost everywhere and $g$ is nonnegative and Lebesgue integrable, then $f$ is Lebesgue integrable and $\int |f|\le \int g$. This follows from monotonicity of the Lebesgue integral. Conversely, $0\le g\le |f|$ with $\int g=\infty$ rules out Lebesgue integrability of $f$. For nonnegative improper integrals, the same comparison applies after taking increasing finite-interval limits.

Complex analysis studies functions of a complex number, especially holomorphic functions. Its distinctive feature is that a single complex derivative controls variation in every direction of the complex plane. The map $f(z)=z^2$ gives a first picture: straight coordinate lines become curved, yet small perpendicular arrows remain perpendicular wherever $f^{\prime }(z)\ne 0$.

The bottom panels show displacements rather than absolute positions. Multiplication by $f^{\prime }(1+i)=2+2i$ predicts the image of each small arrow; the exact image differs by $h^2$. The origin is exceptional because $f^{\prime }(0)=0$. This distinction between local conformal maps and a critical point is central: complex differentiability imposes strong geometric structure, without making every holomorphic function globally one-to-one.

The remarkable next step is that holomorphic functions are determined inside a circle by their values on its boundary. Precisely, let $R>0$, let $D=\{z:|z-a|<R\}$, and suppose $f$ is a holomorphic function on an open set containing the closed disc $\overline{D}$. Orient its boundary counterclockwise. The Cauchy integral formula states that, for every $w\in D$,The contour integral combines boundary values with a kernel depending on the requested interior point. It is a reconstruction formula, not merely a trick for evaluating an integral. The hypothesis about the closed disc ensures that all boundary values and the contour integral are well defined. ^[ref]

At the centre, substituting $z=a+R{e}^{i\theta }$ cancels the denominator against $dz=iR{e}^{i\theta }d\theta$:Thus the centre value is the angular average of the boundary values. For the polynomial $f(z)=z^2$, $a=1+i$ and $R=1/2$, expand the boundary data asBoth oscillating terms have zero average, leaving $f(a)=2i$. Away from the centre the weights change: on the unit circle,With $f(z)=z^2$ and $w=1/2$, expanding the denominator as a geometric series leaves only its second-power term after integration, giving $f(1/2)=1/4$. This also explains how the Cauchy integral formula recovers an entire Taylor series from boundary data.

When a holomorphic function has isolated singularities, a different boundary calculation detects their residues. Suppose $\Gamma$ is a counterclockwise simple closed contour and $f$ is holomorphic on a neighbourhood of it and its interior except at finitely many interior isolated singularities $a_j$. The residue theorem givesA residue is the coefficient of $(z-a_j{)}^{-1}$ in the Laurent series. A short proof idea is to deform the contour to small circles around the isolated singularities using the Cauchy integral theorem. Integrating the Laurent series on each circle kills every power except $-1$. ^[ref]

Here is how a complex pole determines a real improper integral. Set $g(z)=1/(1+z^2)$ and use the contour consisting of the segment from $-R$ to $R$, followed by the upper semicircle $A_R$, with $R>1$. Only the simple pole $i$ lies inside, andThe arc must actually disappear. On $A_R$, the reverse triangle inequality gives $|1+z^2|\ge R^2-1$, and its length is $\pi R$, soThe real integrand is bounded by $x^{-2}$ for $|x|\ge 1$, so both tails converge separately. Taking the limit therefore proves the ordinary improper integralFor example, at $R=10$ the real segment contributes $2\arctan (10)\approx 2.942255$, and the arc contributes $0.199337$, below the bound $10\pi /99\approx 0.317333$. The usefulness of the residue theorem is that the same strategy works for many rational functions whose real antiderivatives are cumbersome: find the enclosed poles, compute their residues, and control the connecting arc. ^[ref]

Historically, Cauchy's 1825 memoir on integrals with imaginary limits develops complex contour integration and its use in transforming and evaluating definite integrals. The opening sections define complex path integrals, discuss changing the path, and examine what happens near singular points. These are recognizable ancestors of the modern link between holomorphic functions, contours and residues. ^[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.

Suppose $h$ is continuously differentiable on $[0,L]$. For $0<x\le L$, the improper integral $g(x)={\int }_0^{x}h(t)(x-t{)}^{-1/2}dt$ exists, but substituting directly into the ordinary Leibniz integral rule produces an undefined endpoint term. First use integration by parts to obtainThe new endpoint term at $t=x$ vanishes, and the differentiated kernel is locally integrable. Differentiating gives the displayed formula. One rigorous justification rewrites the second term as ${\int }_0^x[{\int }_0^s{h}^{\prime }(t)(s-t{)}^{-1/2}dt]ds$ by Fubini's theorem; the inner integral is continuous for $s>0$, so the fundamental theorem of calculus applies. Moreover $g(x)=2h(0)\sqrt{x}+O(x^{3/2})$: if $h(0)\ne 0$, the right difference quotient at zero diverges rather than defining a finite derivative.

Suppose $q(t)=A/(t-c)+r(t)$, where $r$ is continuous, the pole $c$ is fixed, and differentiable endpoints obey $a(x)<c<b(x)$ locally. ThenApply the chain rule and the fundamental theorem of calculus to getThe same boundary expression as in the Leibniz integral rule thus applies to this explicitly specified Cauchy principal value. If $A\ne 0$, the two ordinary one-sided improper integrals diverge, so their sum is not an ordinary integral. A finite boundary derivative alone does not make that divergent integral meaningful.

If $F(t)={\int }_a^{t}f(x)dx$ is bounded and $g$ is a nonnegative differentiable decreasing function tending to zero, then ${\int }_a^{\infty }f(x)g(x)dx$ converges as an improper integral. Integration by parts bounds the tail ${\int }_A^{B}fg$ by a constant multiple of $g(A)$, because $F$ is bounded and ${\int }_A^B|g^{\prime }|=g(A)-g(B)$. In particular ${\int }_1^{\infty }{x}^{-\alpha }\sin xdx$ converges for every $\alpha >0$, even though absolute integrability of that tail requires $\alpha >1$.

For nonevolving comoving number density $n_0$, identical spectra $L_{\nu }\propto {\nu }^{\alpha }$, and no attenuation, the cosmological background intensity from comoving emissivity in an Einstein-de Sitter universe is proportional to ${\int }_0^{\infty }(1+z{)}^{\alpha -5/2}dz$. The improper integral is finite exactly for $\alpha <3/2$, with value $1/(3/2-\alpha )$. Equality gives logarithmic divergence; larger $\alpha$ gives a power-law divergence. A formation epoch or a spectral break can regularize this idealized model.

A pair of oppositely directed contours on a branch cut produces the difference of the two boundary values integrated over their actual endpoint interval. Replacing finite endpoints by zero and infinity requires a limiting argument and convergence of the combined improper integral. Individual side integrals may diverge even when their difference converges, so cancellation must precede separation of endpoint limits.

Consider the vector field $\mathbf{F}=(p(x),q(x),r(z))$ in the open unit cube, where $p$ and $r$ extend to continuously differentiable functions on $[0,1]$, and $q$ is continuous on $(0,1]$ with ${\int }_0^1|q(x)|dx<\infty$. Even if $q$ is unbounded as $x\downarrow 0$, the normal component on an $x$-face has the finite trace $\pm p(x)$. The two $y$-faces have finite improper integral fluxes $\pm {\int }_0^{1}q(x)dx$ and therefore cancel. The remaining face contributions giveTo justify this use of the divergence theorem, first apply it to $[\epsilon ,1]\times [0,1{]}^2$, on which the field is regular. Its divergence is $p^{\prime }(x)+r^{\prime }(z)$, independent of $q$. Its outward flux is $p(1)-p(\epsilon )+(1-\epsilon )[r(1)-r(0)]$. Taking $\epsilon \downarrow 0$ yields the displayed answer, and absolute convergence makes the two $y$-face limits finite separately. For example, $q(x)=1/[x(1+(\ln x{)}^2)]$ diverges at zero butThis argument defines the normal traces and the face integrals explicitly; it does not assert that the whole vector field has a finite value on the singular face.

The fundamental strip of a Mellin transform is the interior of the set of real parts $\sigma$ for whichViewed in the complex plane, it is a vertical strip $\alpha <\mathrm{Re}s<\beta$; either bound may be infinite, and the interior may be empty. Both endpoints of the improper integral must be checked. Convergence on a boundary line, whether absolute convergence or conditional convergence, requires a separate test. The Mellin transform is holomorphic throughout a nonempty fundamental strip.^[ref]

A practical sufficient test uses Big O notation: if $f(x)=O(x^p)$ near zero and $f(x)=O(x^{-q})$ near infinity, with real $p,q$, then comparison with elementary improper integrals givesThus $-p<\sigma <q$ guarantees absolute convergence. These upper bounds on $|f|$ need not locate the maximal Mellin fundamental strip; faster actual decay can widen it. The strip property itself follows by bounding the weight at any intermediate $\sigma$ by the sum of two weights already known to be integrable. On a closed narrower strip, neighboring integrable weights also dominate the natural logarithm factors introduced by differentiation under the integral sign, making the transform a holomorphic function.

An exact two-endpoint example isHere both bounds are sharp: at $\sigma =-2$ the zero-end absolute convergence test contains ${\int }_0^{1}dx/x$, and at $\sigma =3$ the infinity-end test contains ${\int }_1^{\infty }dx/x$. The displayed rational function continues meromorphically outside the Mellin fundamental strip and has simple poles at $-2$ and $3$. Its other finite continued values outside the strip are values of a meromorphic continuation, not of a convergent defining integral. More general endpoint information leads to meromorphic continuation of a Mellin transform from an asymptotic expansion.

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$,

The Gaussian integral gives the exact area under a Gaussian function. For a positive real number $a$,This is a finite improper integral over the whole real line. For example, $I(1)=\sqrt{\pi }\approx 1.772454$ and $I(4)=\sqrt{\pi }/2\approx 0.886227$. Increasing $a$ from one to four narrows the curve and halves its area, while its height at zero stays one. The change of variables $u=\sqrt{a}x$ explains the factor $a^{-1/2}$.

The useful consequence is an exact normalization. For $a>0$ and any real center $\mu$,This is a probability density function for a normal distribution with expected value $\mu$ and variance $1/(2a)$. A narrower probability density function must be taller to keep total probability one. In the usual notation $a=1/(2{\sigma }^2)$, the prefactor is $1/(\sigma \sqrt{2\pi })$, where $\sigma$ is the standard deviation.

The same area calculation makes the one-dimensional heat kernelintegrate to one in $x$, with $D$ the diffusion coefficient. A spreading Gaussian function can therefore preserve total mass while its height falls. This links the Gaussian integral with both the normal distribution and the diffusion equation.

Finite intervals lead to the error function, rather than this single constant:The complementary tail is described by the complementary error function. These definitions and their real-axis limits are given in the NIST reference. ^[ref]

The polar-coordinate proof of the Gaussian integral explains why $\pi$ appears in an integral of a function defined on a line: squaring its area creates a rotationally symmetric integral on a plane. Substitution $t=x^2$ also gives $\Gamma (1/2)=2{\int }_0^{\infty }{e}^{-x^2}dx=\sqrt{\pi }$, connecting the result with the gamma function.