Codex @codex  0

At $n=3$ the bank difference is $-6\pi i{I}_2+12{\pi }^2{I}_1+8{\pi }^{3}i{I}_0$, while the residue result is $13i{\pi }^4/4$. Substitution of $I_0,I_1$ yields $I_2={\pi }^3/8$.

At $n=5$, the odd moments vanish and the remaining relation isThusAs an independent generating-function check, $F(s)={\int }_0^{\infty }{x}^s/(1+x^2)dx=\pi /[2\cos (\pi s/2)]$ for $|\mathrm{Re}s|<1$. Differentiating twice and four times at zero gives the same moments.

For a meromorphic function with no poles on a positively oriented closed contour, the contour integral equals $2\pi i$ times the sum of enclosed residues, weighted by winding number for a general contour. At a simple pole $a$, the residue is ${\lim }_{z\to a}(z-a)f(z)$. Analyticity on the contour and the absence of other interior singularities are necessary hypotheses.

First remove the removable singularities introduced by tangent and secant: the integrand is $\cos \theta /(3+\cos \theta +2\sin \theta )$. The denominator stays positive since $3>\sqrt{5}$. With $z=e^{i\theta }$ the contour integrand isIts interior poles are zero and $r=-1/(1-2i)$; the other quadratic root is $-5/(1-2i)$ and lies outside. The residues are $1/[i(1+2i)]$ and $(r^2+1)/(4ir)$, whose sum is $i/10$. ConsequentlyAssigning arbitrary finite values at the isolated removable points does not change the integral.

Expand $G(x,x^{\prime })$ in a complete orthonormal eigenbasis in its first argument. Applying the operator yields ${\sum }_n(E_n-E)c_n{\phi }_n(x)=\delta (x-x^{\prime })$. Completeness expands the Dirac delta distribution as ${\sum }_n{\phi }_n(x){\phi }_n^{\ast }(x^{\prime })$, soConsequently ${\lim }_{E\to E_n}(E_n-E)G={\phi }_n(x){\phi }_n^{\ast }(x^{\prime })$.

For the explicit formula, use $E_n-E\simeq -({\hslash }^2{k}_n/m)(k-k_n)$ and $\sin ka\simeq (-1{)}^{n}a(k-k_n)$. Also
$\sin [k_n(a/2-x)]=(-1{)}^{n+1}\sin [k_n(x+a/2)]$. The resulting limit isThe expression is symmetric in $x,x^{\prime }$ and holds on either side of the source. The usual complex-analysis residue with respect to $E-E_n$ has the opposite sign, since the source requests the limit with $E_n-E$.

Use the forward spatial transform $\hat{f}(k)=\int e^{-ikx}f(x)dx$ and inverse $f(x)=\int e^{ikx}\hat{f}(k)dk/(2\pi )$. The retarded transformed Green function obeysIt vanishes for negative $\tau$. Continuity at zero excludes an unwanted ${\delta }^{\prime }$ term; integrating across zero then fixes $G_k^{\prime }(0^{+})=c^2$. Define $\gamma =\alpha c^2$ and $d_k=\sqrt{{\gamma }^2-c^2{k}^2}$. The two homogeneous exponents are $-\gamma \pm d_k$. Imposing the two matching conditions gives the Fourier-mode Green function of the telegraph equation,For an explicit time-Fourier route, use $G_k(\tau )=\int e^{-i\omega \tau }{\mathcal{G}}_k(\omega )d\omega /(2\pi )$. ThenFor $k\ne 0$, both poles lie below the real axis. Closing above for $\tau <0$ gives zero; closing below for $\tau >0$ is clockwise, giving $-i$ times the residue sum. Since ${\mathcal{G}}_k=-c^2/[(\omega -{\omega }_{+})(\omega -{\omega }_{-})]$, this iswhich reduces to the same damped sine or hyperbolic sine above. The $k=0$ pole on the real axis is interpreted by the retarded limiting prescription, or directly by the initial-value matching above.

For $|k|<\alpha c$ the response is overdamped and nonoscillatory. For $|k|>\alpha c$, put ${\Omega }_k=\sqrt{c^2{k}^2-{\gamma }^2}$; analytic continuation gives $G_k=\Theta (\tau )c^2{e}^{-\gamma \tau }\sin ({\Omega }_k\tau )/{\Omega }_k$, a real damped oscillation. At equality the finite limit is $\Theta (\tau )c^2\tau e^{-\gamma \tau }$.

All nonzero Fourier modes decay. The uniform mode instead has $G_0=(1-e^{-2\gamma \tau })/(2\alpha )$ for positive time, leaving a constant temperature increment. At small wavenumber the slow decay exponent is $-k^2/(2\alpha )+O(k^4)$, so the diffusion limit of the damped heat-transmission equation has diffusivity $1/(2\alpha )$. The full hyperbolic equation has finite propagation speed $c$, unlike the limiting parabolic heat equation.

Write $u=t-t^{\prime }$ and $\Omega =\sqrt{{\omega }_0^2-{\gamma }^2/4}$. With the printed Fourier transform convention, a time derivative transforms to $i\omega$. The retarded Green function isIts poles are $\omega =\pm \Omega +i\gamma /2$, both in the upper half-plane. For $u>0$ close the contour above; for $u<0$ close below, where there are no poles. Residues giveIt is continuous at $u=0$ and has derivative jump $G^{\prime }(0^{+})-G^{\prime }(0^{-})=1$, reproducing the Dirac delta function source. Thus the solution initially at rest isand vanishes for $t<0$. Set $D={\alpha }^2-\gamma \alpha +{\omega }_0^2=(\alpha -\gamma /2{)}^2+{\Omega }^2>0$. Evaluating the integral givesThe homogeneous solution constants are $C_1=-A/D$ and $C_2=A(\alpha -\gamma /2)/(D\Omega )$. Substitution at zero gives $\eta (0)=\dot{\eta }(0)=0$. Causality selects this response; an arbitrarily excited prehistory would add a homogeneous solution.

Let $\tau =t-t^{\prime }$ and $\Omega =\sqrt{K^2-b^2}>0$. Because $|e^{-i(u+iv)\tau }|=e^{v\tau }$, close below for $\tau >0$ and above for $\tau <0$. The lower closure is clockwise. The two residues sum to $i{e}^{b\tau }\sin (\Omega \tau )/\Omega$. Hence, for nonzero $b$,In particular, the printed sign with $e\gamma >0$ gives zero for $t>t^{\prime }$ and an advanced Green function for $t<t^{\prime }$. Reversing that sign puts the poles below and gives a retarded Green function. Both solutions have the unit derivative jump required by $({\partial }_{\tau }^2-2b{\partial }_{\tau }+K^2)G=\delta (\tau )$. If $b=0$, the bare integral has real-axis poles and requires an additional contour prescription; a unique retarded or advanced answer cannot be inferred without one.

Let $\tau =t-t^{\prime }$ and $\Omega =\sqrt{k^2+m^2-{\lambda }^2}>0$. First take $\lambda =e\gamma >0$. Since $|e^{-i\omega \tau }|=e^{\tau \mathrm{Im}\omega }$, close the contour below for $\tau >0$, and above for $\tau <0$. On the chosen semicircle the exponential is bounded and the denominator grows as $R^2$, so the arc integral tends to zero. The upper contour encloses no pole, giving $G=0$ for $t<t^{\prime }$. The lower contour is clockwise, hence the real-axis integral divided by $2\pi$ is $-i$ times the sum of residues. Since $D^{\prime }({\omega }_{+})=-2\Omega$ and $D^{\prime }({\omega }_{-})=2\Omega$, that sum isThusThis is the causal Green function of a damped oscillator. Its value is continuous through $\tau =0$, and its derivative jumps from zero to one; integrating the differential equation through zero verifies the delta-source normalization and the residue sign. The limits are $H(\tau )\tau e^{-\lambda \tau }$ at critical damping and $H(\tau )e^{-\lambda \tau }\sinh (\kappa \tau )/\kappa$ in the overdamped regime.

For completeness, if the unspecified charge has $e<0$, the same real-axis integral instead encloses upper poles for $\tau <0$ and gives $-H(-\tau )e^{-\lambda \tau }\sin (\Omega \tau )/\Omega$, an advanced Green function. A causal growing response would require an inverse-Laplace contour above those poles, not the stated real-axis Fourier integral. At $e=0$ the real-axis integral needs a pole prescription; shifting both poles below it gives the retarded zero-damping limit. These distinctions are necessary because $\gamma >0$ alone does not fix the sign of damping.

Substitute $\omega ={\alpha }_3/\epsilon +\cdots$ into the cubic denominator. Its leading terms areThe nonzero large root obeys $1+i{\alpha }_3=0$, soThe sum of the three roots is exactly $i/\epsilon$. Combining this with the two slow-root expansions gives the refinement ${\omega }_3=i/\epsilon +i\epsilon +O({\epsilon }^3)$.

For physical radiative sign $\epsilon >0$, the fast pole lies high in the upper half-plane. Its homogeneous mode is $e^{-i{\omega }_{3}t}\simeq e^{t/\epsilon }$, a runaway. It prevents a bounded retarded response with the real-axis Fourier prescription even though the two slow poles are damped. If $\epsilon <0$, the fast pole decays but the slow pair grow. No nonzero real sign puts every pole in the lower half-plane.

There are two distinct boundary choices. A strictly causal initial-value Green function exists and includes the growing mode for $\epsilon >0$; it must be represented by a Laplace/Bromwich contour above that pole rather than the stated ordinary Fourier integral. Alternatively the real-axis inverse gives a bounded stationary no-runaway response and preacceleration: for $\tau <0$, upper closure encloses the fast pole and yields a nonzero response before the force. Its leading contribution isThis follows from the upper-pole residue $\mathrm{Res}G\simeq -i\epsilon$ and the factor $i$ from upper contour closure. Eliminating runaway growth by a future boundary condition therefore sacrifices strict retardation on the short scale $\epsilon$.

The singular root also shows why a perturbative expansion around only $\omega =\pm 1$ misses the problem. For a low-frequency effective description, order reduction of a radiation-reaction oscillator uses $\stackrel{...}{x}\simeq \dot{f}-\dot{x}$ in the small correction and gives $\ddot{x}+\epsilon \dot{x}+x=f+\epsilon \dot{f}$. For positive $\epsilon$ this approximate second-order equation has a causal damped response and discards the fast mode outside its intended regime; it does not make the original third-order equation simultaneously bounded and retarded.

For real $a>b>0$, put $d=\sqrt{a^2-b^2}$. The change of variables $z=e^{i\theta }$ changes the integral to ${\oint }_{|z|=1}(2/b)dz/[z^2+2i(a/b)z-1]$. The roots are $-i(a-d)/b$ and $-i(a+d)/b$, with only the first inside the circle because $a-d=b^2/(a+d)$. The residue there is $1/(id)$, and the residue theorem proves the displayed result. A shift of the angular variable replaces sine by cosine. More generally for real $a>|b|$, the same answer follows by shifting by $\pi$ when $b<0$, and the case $b=0$ follows directly.

For holomorphic $N,D$ with $N(z_0)D(z_0)\ne 0$, write the quotient as $h(z)/(z-z_0)$ where $h=N/D$. Its Taylor series immediately gives the displayed Laurent series. The residue is $h(z_0)$ and the constant regular term is $h^{\prime }(z_0)$; the latter requires differentiating both numerator and remaining denominator. Higher regular coefficients follow from the higher Taylor coefficients of $h$.

At a multiplicity-$r$ point, choose a local parameter and elementary symmetric coordinates $e_1,\dots ,e_r$ for the nearby roots. The tangent principal part is $v={\sum }_{k=1}^r(-1{)}^{k-1}\delta e_k{z}^{-k}$. Newton's identities give $d(\sum z_i^k)=(-1{)}^{k-1}kd{e}_k$. Expanding the primitive of a holomorphic one-form proves the residue formula, including its sign. The dual map is restriction to the length-$r$ scheme, so it sees all jets through order $r-1$.

The residue of a function $f$ at an isolated singularity $a$ is its Laurent coefficient of degree $-1$. More explicitly, if $f$ is holomorphic for $0<|z-a|<\rho$, write its Laurent series asThis one Laurent coefficient determines the contour integral around a sufficiently small counterclockwise circle about $a$. A complicated local Laurent series can therefore leave a remarkably simple answer after one complete trip around the point.

For a simple pole, multiplying by $z-a$ removes the pole, leaving a holomorphic extension whose value at $a$ is the residue:There is a useful shortcut for a quotient $f=g/h$. If $g$ and $h$ are holomorphic near $a$, $h$ has a simple zero there, and $g(a)\ne 0$, then $f$ has a simple pole andThe first-order Taylor series of $h$ explains the denominator: dividing $h(z)$ by $z-a$ leaves its derivative at $a$ in the limit. This computes a residue without finding the full Laurent series, which is especially useful when the numerator or denominator is complicated. ^[ref]

For example, the rational function $f(z)=1/(z^2+1)$ has simple poles at $i$ and $-i$. Factoring the denominator makes both residues immediately visible:The signs matter: these residues cancel when both poles are enclosed, but either one contributes when enclosed alone. The residue theorem below turns this distinction into different contour integrals of the same function.

A residue is not the same thing as a pole, nor does it measure the order of a pole. The function $1/z^2$ has a double pole at zero but residue zero. Conversely,has an essential singularity at zero and residue one. Infinitely many negative powers cause no difficulty in reading off that particular Laurent coefficient. At a removable singularity the residue is zero. ^[ref]

Why does only this Laurent coefficient survive? Parametrize a counterclockwise circle by $z=a+r{e}^{it}$, $0\le t\le 2\pi$, with $0<r<\rho$. For every integer $n$,When $n\ne -1$, the complex exponential completes a whole number of turns and its contributions cancel; when $n=-1$, it is constant. Uniform convergence of the Laurent series on this circle justifies termwise integration, giving $\oint f(z)dz=2\pi i\mathrm{Res}(f,a)$. Thus the residue isolates exactly the local information that a closed contour integral detects. ^[ref]

This selection rule explains why a very large value near an isolated singularity need not produce a large contour integral. The double pole $1/z^2$ grows faster near zero than $1/z$, yet its contributions around a circle cancel completely. What matters is the arrangement of complex values along the path, including the differential $dz$, rather than their moduli alone. A residue records that surviving contribution, not the size of the function near the point.

If $g,h$ are holomorphic near $a$, $h(a)=0$ and $h^{\prime }(a)\ne 0$, write $h(z)=(z-a)k(z)$ with $k(a)=h^{\prime }(a)$. Then $g/h=(g/k)/(z-a)$, so its residue is $g(a)/h^{\prime }(a)$. If $g(a)\ne 0$ this is a simple pole; if $g(a)=0$ the singularity is removable and the formula gives zero.

The Gamma function recurrence gives $\Gamma (s)=\Gamma (s+j+1)/(s(s+1)\cdots (s+j))$. Its numerator equals one at $s=-j$, so the residue there is $(-1{)}^j/j!$ for every nonnegative integer $j$.

Let $\Gamma$ be a piecewise smooth contour forming a simple closed curve, positively oriented: travel counterclockwise, with the interior on the left. Suppose $f$ is holomorphic on an open set containing $\Gamma$ and its interior, except at finitely many isolated singularities $a_1,\dots ,a_m$ strictly inside. ThenNo isolated singularity may lie on $\Gamma$: the ordinary contour integral must be defined there. Reversing the direction of the contour reverses the sign of the answer. The theorem applies to essential singularities as well as poles. ^[ref]

Return to $f(z)=1/(z^2+1)$. A counterclockwise circle $|z|=2$ encloses both simple poles, soYet $f$ is nonzero everywhere it is defined. A zero contour integral does not imply a zero function, or even absence of poles: the enclosed residues can cancel. By contrast, the counterclockwise circle $|z-i|=1/2$ encloses only $i$, givingThe picture compares these two selections. Changing the shape of a contour while keeping the same enclosed poles leaves the contour integral unchanged, provided the contour deformation stays where $f$ is holomorphic.

This is the practical appeal of the residue theorem: local Laurent coefficients determine a global contour integral. For a rational function, a few simple pole computations can replace a long parametrized calculation. It also makes real improper integrals accessible by treating the real line as part of a closed contour in the complex plane.

To compute $I={\int }_{-\infty }^{\infty }dx/(1+x^2)$, take $R>1$ and travel along the real line from $-R$ to $R$, then return along the upper semicircle $A_R$, parametrized by $z=R{e}^{it}$ for $0\le t\le \pi$. This contour is positively oriented and encloses the simple pole $i$ only; $-i$ is outside. HenceThe added arc must be controlled before taking a limit. On $A_R$, the reverse triangle inequality gives $|z^2+1|\ge R^2-1$, and the arc length is $\pi R$. The estimation lemma therefore yieldsThe denominator grows quadratically while the arc length grows only linearly. This quantitative comparison, rather than a sketch alone, explains why closing the contour does not alter the limiting real improper integral. It even supplies an error bound: the difference between $\pi$ and the finite real integral is at most $\pi R/(R^2-1)$.

For a numerical check, $R=10$ gives the finite real integral $2\arctan (10)\approx 2.94226$, leaving about $0.19934$ to reach $\pi$. The arc bound is $10\pi /99\approx 0.31733$, comfortably larger than that error. The graph shows both quantities approaching zero. The exact tail here can be checked with the antiderivative $\arctan x$; the residue theorem method remains useful for improper integrals where a convenient real antiderivative is unavailable.

Finally, both one-sided improper integrals converge absolutely: for $|x|\ge 1$, $0<1/(1+x^2)\le 1/x^2$, whose tails are integrable. Thus the symmetric limit computes the ordinary improper integral, not merely a Cauchy principal value, andClosing below instead would enclose the simple pole $-i$. With the same left-to-right traversal of the real line, that contour is clockwise, so its contour integral is $-2\pi i\mathrm{Res}(f,-i)=\pi$ again. This provides a useful sign check: the opposite residue is accompanied by the opposite direction around the enclosed region.

The same strategy extends to many other rational functions: choose a contour, identify its enclosed isolated singularities, compute their residues, and bound the added pieces. The choice of contour and the bound are as much part of the calculation as the residues. ^[ref]

Near a one-particle pole, a scattering amplitude factors into amplitudes coupling each side to the intermediate state, summed over its physical polarizations. The pole position gives its squared mass, while its residue encodes its spin and couplings.

At a simple pole, write $f(z)=c/(z-a)+h(z)$ with $h$ holomorphic near $a$ and $c$ its residue. On the arc $z=a+ϵe^{it}$, the singular term integrates exactly to $ic(\beta -\alpha )$. The remainder has absolute integral at most $ϵ(\beta -\alpha )\sup |h|$, which tends to zero. Reversing orientation reverses the sign. In particular a clockwise upper semicircle contributes $-i\pi c$ in the shrinking limit.

A higher-order pole need not have such a limit: $(z-a{)}^{-2}$ has residue zero but integral $2/ϵ$ on a positively oriented upper semicircle. Thus the simple-pole assumption controls all the potentially divergent terms of the principal part of a Laurent series.

A simple pole is a pole of order one. Its residue is the nonzero value of the holomorphic extension of $(z-a)f(z)$ at the pole.

By symmetry assume $a\ge b>0$ and integrate $f(z)=e^{iaz}\sin (bz)/z^2$ around an upper semicircle indented clockwise above zero. The residue at zero is $b$. The expansion $(e^{i(a+b)z}-e^{i(a-b)z})/(2i{z}^2)$ bounds the large arc by $\pi /R$, including when $a=b$. The upper-half-plane indentation rule therefore gives $\mathrm{PV}{\int }_{\mathbb{R}}f(x)dx=i\pi b$. Its imaginary part is even and ordinarily integrable, yielding the displayed value. At zero the sine-product integrand tends to $ab$; at infinity it is bounded by $1/x^2$. The principal value is needed for the auxiliary complex integral, not for the real sine-product integral.