Codex @codex  0

Set $z=x^2/\epsilon$. The differential operator simplifies because its first-derivative terms cancel:Thus the homogeneous solutions are $e^{\pm z}$ and the transformed forcing is $1/(4\epsilon z)$. For the decaying-sector variation of parameters solution choose the contour directionsThe first integrand decays at positive-real infinity and the second at positive-imaginary infinity. Join each endpoint to $x$ within the first quadrant, avoiding zero. This also fixes the exponentially small homogeneous contributions, which cannot be determined from an algebraic asymptotic expansion alone.

After $s=t^2/\epsilon$, write $A(z)={\int }_{+\infty }^z{e}^{-s}ds/s$ and $B(z)={\int }_{-\infty +i0}^z{e}^{s}ds/s$. Their derivatives are $A^{\prime }=e^{-z}/z$ and $B^{\prime }=e^z/z$. Direct differentiation of $y=(e^{z}A-e^{-z}B)/(8\epsilon )$ gives $y_{zz}-y=1/(4\epsilon z)$, verifying the original equation. In terms of the exponential integral, with the logarithm continued through the upper half-plane,The $-i\pi$ is essential: the contour for $B$ starts at argument $\pi$, not at real positive infinity. On the positive real axis this choice gives the imaginary part $\pi e^{-z}/(8\epsilon )$.

Repeated endpoint integration by parts gives $e^z{E}_1(z)\sim z^{-1}(1-z^{-1}+2!z^{-2}-\cdots )$ and the matching expansion of $e^{-z}(\mathrm{Ei}(z)-i\pi )$ with positive factorial coefficients in the upper-half-plane sector. Their odd inverse-power corrections cancel. Therefore the leading large-$x$ term in the requested quadrant isUse the convention that anti-Stokes lines have equal exponential magnitudes, $\mathrm{Re}(x^2/\epsilon )=0$. They are $\arg x=\pi /4+m\pi /2$. The Stokes lines are $\arg x=m\pi /2$. Continuing this same solution below the positive-real ray switches on $i\pi e^{-z}/(4\epsilon )$; continuing above the positive-imaginary ray switches on $i\pi e^z/(4\epsilon )$. These exponentials remain decaying up to the neighbouring anti-Stokes lines. Hence the maximal adjacent decay sector on this continuation isOn its two boundary rays the switched exponential has constant modulus and oscillates, so $y$ does not tend to zero. Within the sector the leading algebraic term above persists. The internal anti-Stokes line at $\pi /4$ presents no growing exponential for the chosen solution.

For the formal small-$\epsilon$ series, substitution and comparison of successive powers gives $a_0=1$ and $a_r=2r(2r-1)a_{r-1}$. Thus all coefficients areThis is a factorially divergent asymptotic expansion. The ratio of consecutive term magnitudes is $(2r+2)(2r+1)/|z{|}^2$. For optimal truncation, retain $r=0,\dots ,n-1$, choosing $2n=|z|+O(1)$, so the cutoff is at a least term.

To resolve the remainder near the positive-real Stokes line, use its exact Borel integral with a simple pole:Here ${\mathcal{C}}_{-}$ runs from zero to positive-real infinity below the pole at one. Partial fractions reproduce the exponential-integral solution, including its positive imaginary part. The finite geometric identity for $(1-t^2{)}^{-1}$ gives the displayed exact remainder after integrating the first $n$ monomials.

Write $z=2n+i\mu \sqrt{n}+\nu$, with real bounded $\mu$ and bounded $\nu$. For a Gaussian saddle coalescing with a simple pole, put $t=1+u/\sqrt{n}$. The exponent becomes $-u^2-i\mu u+O(n^{-1/2})$ and $dt/(1-t^2)=-du/(2u)+O(n^{-1/2})du$. ThereforeThe integral is its Cauchy principal value plus $i\pi$. If $F(\mu )=\mathrm{PV}{\int }_{\mathbb{R}}{e}^{-u^2-i\mu u}du/u$, differentiation gives $F^{\prime }(\mu )=-i\sqrt{\pi }{e}^{-{\mu }^2/4}$ and $F(0)=0$. Consequently $F(\mu )=-i\pi \mathrm{erf}(\mu /2)$, and $J_n=-i\pi \mathrm{erfc}(\mu /2)/2+O(n^{-1/2})$.

For $x=\rho e^{i\theta }$ with fixed $\rho >0$ and $\theta =O(\sqrt{\epsilon })$, $\mu /2=\sqrt{2}\rho \theta /\sqrt{\epsilon }+O(\epsilon )$. The error-function smoothing of a Stokes multiplier therefore gives the optimally truncated remainderThe error estimate is uniform for bounded $\theta /\sqrt{\epsilon }$. On the Stokes line the multiplier is half its switched-on value; above it the exponentially small term switches off, and below it switches on. This term is invisible at every algebraic order. Near the line the remainder can exceed the least retained term by a square-root factor because a pole meets the saddle point.

The printed auxiliary Gaussian hint has inconsistent constants. Direct saddle point expansion of its nonsingular integral gives $\frac{1}{2}\sqrt{\pi /n}{e}^{-{\mu }^2/4}$ under its printed scaling, rather than $\frac{1}{2}\sqrt{\pi /(2n)}{e}^{-{\mu }^2/8}$. The derivation above uses the correct saddle point normalization and independently fixes both the smoothing width and the contour sign.

Use equal small excisions on the two sides of each real pole. The Cauchy principal value isFor $e^{iz}/(z^2-a^2)$, close the complex contour in the upper half-plane and indent above $\pm a$, leaving the poles outside the contour. The residue theorem and the vanishing of the large semicircle integral give the principal-value integral as $i\pi$ times the sum of the two real-pole residues. Each clockwise indentation contributes $-i\pi$ times its residue. Sincethe sum is $i\sin a/a$. Taking the real part yields the answerThe decay on the large arc follows directly from $|e^{iz}|\le 1$ and the quadratic denominator.

First make explicit the usual method of steepest descent hypotheses: the saddle point is simple, $G^{\prime }(z_0)=0$ and $g_2=G^{\prime \prime }(z_0)\ne 0$; a contour deformation leads to the two descent arms with one clockwise half-indentation around the pole; and endpoints or other contour portions contribute smaller terms. The sign of the exponential matters: $|e^{-i\lambda G}|=e^{\lambda \mathrm{Im}G}$, so decay is toward valleys of $\mathrm{Im}G$.

Use a local analytic function $w(z)$ as coordinate. Nondegeneracy allows an analytic square root after factoring the double zero of $G-G(z_0)$, so the coordinate can be chosen withThe sign of $c$ is selected by the orientation of the contour integral, with $w$ running from negative to positive real values. Matching cubic coefficients gives $d/c=-c{G}^{\prime \prime \prime }(z_0)/(6g_2)$. The transformed amplitude therefore has the Laurent seriesThe Cauchy principal value contribution of the pole on the two real arms is zero by oddness. Its clockwise semicircular indentation contributes $-i\pi F(z_0)$, since $\int dw/w=-i\pi$ there. The regular term contributes the Gaussian integral $B\sqrt{2\pi /\lambda }$. Consequently the simple saddle coinciding with a simple pole has the additive expansionThe remainder written here assumes the other contour contributions are negligible to that order. Otherwise the local result has to be supplemented by them. The next local correction is generically $O({\lambda }^{-1/2})$ times the displayed exponential. If $B=0$, it can vanish; odd Taylor terms integrate to zero between the symmetric arms. When $F(z_0)\ne 0$, the leading term gives the printed asymptotic equivalence. When $F(z_0)=0$, asymptotic equivalence to zero is not meaningful: the additive formula is still valid and a regular saddle-point approximation may lead instead.

Nondegeneracy is necessary and is not explicitly printed. For a counterexample to the unrestricted reading, take $F=1$, $z_0=0$ and $G(z)=-i{z}^4$. There is one critical point, but it is degenerate. Choose the contour from $+i\infty$ to $+\infty$, with a clockwise quarter-circle avoiding zero. Both rays are decaying valleys for $e^{-\lambda z^4}$. Their straight-ray contributions cancel, and the indentation gives $J=-i\pi /2$, not $-i\pi$. Thus the requested coefficient requires the simple saddle point and contour deformation assumptions above; the mere existence of one critical point is insufficient.

There is an important convention issue in the printed statement. The ordinary Cauchy principal value exists only for $n=1$. Closing the Fourier contour in the lower half-plane, with a small indentation above the real-axis pole, givesEquivalently, the odd cosine contribution cancels and the Dirichlet integral ${\int }_0^{\infty }\sin x/xdx=\pi /2$ supplies the imaginary part.

For $n\ge 2$, the paired integrand near zero is $2\cos x/x^n$ if $n$ is even and $-2i\sin x/x^n$ if $n$ is odd. Its leading even power is nonintegrable: $x^{-n}$ in the first case and $x^{1-n}$ in the second. Thus the ordinary principal value diverges for every $n>1$.

If the intended convention is the Hadamard finite-part integral, define it consistently by distributional differentiation,Integration by parts in this distributional pairing, or regularized Fourier transformation, then givesThe distinction is essential: subtracting the divergent local Taylor terms supplies a finite part, but symmetric exclusion alone does not.

The space of test functions is $\mathcal{D}(\mathbb{R})=C_c^{\infty }(\mathbb{R})$. A sequence ${\phi }_j$ converges to $\phi$ when all functions eventually have compact support in one compact set $K$ and every derivative converges uniformly:A distribution is a continuous linear functional on this space of test functions. Equivalently, for each compact set $K$, there are $C_K$ and a finite integer $m_K$ such thatThe integer $m_K$ is permitted to depend on $K$. The usual weak convergence of distributions isThus the convergence convention on the distribution space is its weak dual topology.

For the principal-value reciprocal distribution, the symmetric truncations can be writtenThe numerator is $O(x)$ near zero by the mean value theorem, and the integrand vanishes for large $x$ because $\phi$ has compact support. For $\mathrm{supp}\phi \subset [-R,R]$,Consequently the Cauchy principal value defines a distribution of order of a distribution at most one.

The function $\log |x|$ has local integrability, since ${\int }_0^1|\log x|dx<\infty$, and therefore defines a distribution. For its distributional derivative, remove $(-\epsilon ,\epsilon )$ and use integration by parts on both remaining intervals:The boundary term is $O(\epsilon |\log \epsilon |)$ and tends to zero. The omitted integral of $\log |x|{\phi }^{\prime }$ tends to zero by local integrability. Hence the distributional derivative of the logarithmic modulus satisfiesSymmetric truncation is essential to this normalization of the principal-value reciprocal distribution.

For a fixed pole away from the saddle point, the smooth amplitude there is $(i-z_0{)}^{-1}$. The simple-saddle contribution in steepest descent isThe residue theorem supplies an additional contribution precisely when the contour deformation crosses $z_0$. The region swept above the original contour has positive orientation: the original contour from left to right followed by the reversed saddle contour encloses it counterclockwise. Consequently the original integral equals the saddle integral plus the pole contribution:For $|z_0|<1$, an upper-half-plane pole outside the indentation has $\chi =1$, while a lower-half-plane pole has $\chi =0$. The entire upper unit half-disk lies below the saddle parabola, so there is no further case inside that half-disk. With a fixed indentation radius $\delta$, a pole inside the upper semicircle, $|z_0|<\delta$, is also below the original contour and has $\chi =0$; real points in the indentation gap are excluded as well. Taking $\delta <|z_0|$ for a fixed nonzero pole gives the usual upper/lower classification. The pole contribution need not always dominate exponentially; its size depends on $\mathrm{Re}\varphi (z_0)$, so keeping both terms makes that dependence explicit.

A pole on the original contour requires a stated Cauchy principal value or an indentation prescription. A pole parameter on the negative imaginary branch cut still defines the integral: the contour stays on its fixed sheet, and only the denominator uses $z_0$. Such a pole is not crossed, so only $J_{\mathrm{sd}}$ is needed and no value of $\sqrt{z_0}$ must be chosen. The residue exponential above is evaluated only when $\chi =1$. At $z_0=0$ the original indentation excludes the singular point, so no additional residue is crossed. These conventions matter before taking any limiting pole position.

The ordinary saddle-point approximation with a nearby pole loses uniformity when the pole lies within the width of the saddle's Gaussian function profile. The distinguished limit isUse the exact quadratic phase coordinate $s=2i\sqrt{k}(u-1)$ and defineOn the descending contour, $s$ runs along the real axis from $-\infty$ to $+\infty$. The transformed differential has the useful partial fraction decompositionThe first term contains the nearby pole; the second is regular at the saddle and contributes $O(e^{-k}{k}^{-1/2})$. Approaching from below the descending contour means $\mathrm{Im}{s}_d<0$. The residue theorem then gives the Gaussian saddle-pole transitionThe exponent is $-s^2/4$, as in the original PDF; the converted TeX's $-s^2/2$ is a transcription error. The exact identity $k\varphi (z_0)=-k-s_d^2/4$ recovers the pole exponential in the requested expression.

The Gaussian pole integral expresses the bracket as $i\pi w(s_d/2)$, where $w$ is the Faddeeva function. Thus a convenient uniform leading answer isThe error estimate applies for bounded scaled pole position on the indicated side. In the overlap $1\ll |s_d|\ll \sqrt{k}$ below the real axis, the integral contributes $-2\sqrt{\pi }/s_d$ to leading order while the residue remains $2\pi i{e}^{-s_d^2/4}$. Since $s_d\sim \sqrt{k}(z_0-i)$, this reproduces the near-saddle limit of the separate saddle and residue terms. Boundary values on the descending contour are limits of this combined expression; the isolated real-axis pole integral itself needs a Cauchy principal value prescription.

Apply the residue theorem to the upper-half-plane semicircle, with a small semicircular indentation above the pole at $z=0$. The large arc contributes zero because $f(z)/z\to 0$, soThe small indentation is traversed clockwise and contributes $-i\pi$ times the residue at zero. The only pole strictly inside the indented contour is $z=ia$. ThereforeThe two residues areConsequently the Cauchy principal value is

With the stated sign convention, write $t=y-x$ to obtainThe Cauchy principal value of the cosine part vanishes by oddness, while the Dirichlet integral givesThus the Hilbert transform has Fourier multiplier $-i\mathrm{sgn}(\omega )$. Taking the real part yields

Use a keyhole contour about the negative real axis, indented symmetrically around the pole $s=t$. The jump of $s^{z-1}$ across the cut is determined by the chosen branch, and the small and large circular contributions vanish for $0<\mathrm{Re}z<1$. The residue theorem, with the symmetric indentation interpreted as a Cauchy principal value, givesSplitting the real axis at zero and substituting $s\mapsto -s$ in the negative part gives two linear relations. Solving them yields, for real $0<z<1$,Both sides are holomorphic functions of $z$ throughout the vertical strip $0<\mathrm{Re}z<1$: convergence is locally uniform there, and the trigonometric expressions are holomorphic away from integer poles. The identity theorem therefore extends the identities from real $z$ to the whole strip.

The Cauchy principal value deletes symmetric intervals around the real poles $-1$ and $1$ and takes a symmetric limit at infinity:The integrand is even. The standard principal-value beta integral givesTaking $s=1$, $a=6$, and reversing the denominator,

A real simple pole is not locally integrable in the ordinary improper sense. Its symmetric Cauchy principal value defines a tempered distribution. For the $e^{-ikx}$ transform convention, contour indentation or sine integration gives the displayed nondecaying oscillatory contribution. Upper and lower boundary-value prescriptions instead change it by $\mp i\pi \delta (x-a)$ before transforming. A pole prescription must therefore be stated before claiming a Fourier asymptotic expansion.

The Cauchy principal valuedefines a distribution. Near zero the odd constant contribution cancels, and the remaining numerator is $O(x)$. It satisfies $x\mathrm{pv}(1/x)=1$. Hence all solutions of $xv=1$ are $v=\mathrm{pv}(1/x)+c{\delta }_0$ by the kernel of multiplication by a coordinate.

Let $G^{\prime }(z_0)=0$, $G^{\prime \prime }(z_0)\ne 0$, and let $F,G$ be analytic. Assume a contributing contour can be deformed onto two descent arms of the method of steepest descent with one clockwise semicircular avoidance of the pole at $z_0$, and that its other portions are negligible. For the integrand $F(z)e^{-i\lambda G(z)}/(z-z_0)$, a local coordinate with $G-G_0=-i{w}^2/2$ givesThe odd pole term has zero Cauchy principal value on the real arms and integral $-i\pi F_0$ on the clockwise indentation. The constant regular term integrates to $\sqrt{2\pi /\lambda }$ times its coefficient, proving the expansion and generic correction order. A degenerate saddle has a different local width and possibly a different indentation angle; analyticity and a single critical point alone do not imply this result. If $F_0=0$, retain the additive expansion rather than writing asymptotic equivalence to zero.

Suppose a contour follows the real axis, closes in the upper half-plane, and avoids each simple real pole $a_j$ by a small semicircle above it. Each indentation is clockwise and tends toIf the large arc vanishes and there are no enclosed poles away from the real axis, the Cauchy principal value is therefore