Codex @codex  0

For a piecewise differentiable parametrized contour $z(t)$, the contour integral is $\int f(z(t))z^{\prime }(t)dt$. It depends on the orientation of the contour. For a holomorphic function, the Cauchy integral theorem allows a contour deformation through regions without singularities.

To eliminate an unknown boundary transform by a dispersion relation symmetry, evaluate the global relation at $\nu (k)$ only where the spatial transform remains analytic. Substitute the resulting trace identity into the contour representation. An unwanted transform vanishes only if its transformed argument and exponential admit a valid contour deformation and decay estimate. For constant drift, $\nu (k)=i\alpha -k$ maps the upper domain $\mathrm{Re}(k^2-i\alpha k)<0$ into the lower transform half-plane.

For a simple saddle point, a nearby pole coalesces with its contributing region when their separation is of order the inverse square root of the large parameter. A local quadratic phase coordinate changes the singular part to a Gaussian pole integral. Adding any residue crossed by the contour deformation gives a Faddeeva function rather than two separately singular approximations. The regular part of the transformed amplitude remains smaller by the ordinary saddle width.

A parabola is a plane curve whose points are equally distant from a focus and a directrix. In suitable coordinates it has equation $y=a{x}^2$, with $a\ne 0$. Its quadratic geometry can also describe a contour deformation through a saddle point.

For example, take a focus of a conic at $(1,0)$ and a directrix of a conic $x=-1$. Equality of squared Euclidean distances gives $(x-1{)}^2+y^2=(x+1{)}^2$, hence $y^2=4x$. The point $(1,2)$ has distance $2$ to both the focus and the directrix. Its eccentricity of a conic is one. Unlike an ellipse, a parabola is unbounded and has only one matched focus-directrix description.

For a self-adjoint C-star element $h$, the C-star identity gives $\Vert h^2\Vert =\Vert h{\Vert }^2$. If $x$ is a normal C-star element, thenAll powers are normal C-star elements, so induction gives $\Vert x^{2^j}\Vert =\Vert x{\Vert }^{2^j}$. The general spectral radius formula now gives $r(x)=\Vert x\Vert$.

One can also see the needed growth estimate directly. For $R>r(x)$ the resolvent Cauchy coefficient formula givesThe formula follows by the Neumann series on a larger circle and contour deformation through the resolvent annulus. Put $m=2^j$, take $m$th roots and let $j\to \infty$ to get $\Vert x\Vert \le R$. Let $R\downarrow r(x)$; the reverse inequality already follows from the norm bound on the algebra spectrum. ThusThis proves spectral radius norm equality for normal elements from the norm identity.

Use a nonzero complex unital Banach algebra $B$ with submultiplicative norm, and write its unit as $e$. The complex scalar field is the usual convention in this spectral statement; for a real algebra one uses its complexification. We derive the needed analytic facts from norm-convergent series and scalar complex analysis.

If $\Vert a\Vert <1$, completeness makes the Neumann series ${\sum }_{j=0}^{\infty }{a}^j$ converge. Multiplying its finite partial sums by $e-a$ gives $e-a^{m+1}$, which tends to $e$. Hence its sum is the inverse of $e-a$. If $b$ is invertible, then $b+h=b(e+b^{-1}h)$ is invertible whenever $\Vert b^{-1}h\Vert <1$. This proves that the invertible elements form an open set and that inversion is locally norm-continuous.

Define the Banach algebra spectrum and the resolvent of an element byFor $|\lambda |>\Vert x\Vert$, the Neumann series givesThus $\sigma (x)$ is bounded by $\Vert x\Vert$ and is closed, because invertibility is open. Also $R(\lambda )\to 0$ in norm at infinity; for example its norm is at most $\Vert e\Vert /(|\lambda |-\Vert x\Vert )$, since $\Vert e\Vert \ge 1$ in a nonzero submultiplicative unital algebra. The standard normalization $\Vert e\Vert =1$ is not needed.

For points where the Banach algebra resolvent exists $\lambda ,\mu$, multiply out the inverses to get the resolvent identityFor ${\lambda }_0$ outside the Banach algebra spectrum, put $R_0=R({\lambda }_0)$. Locally,This series and its differentiated series converge uniformly on every smaller disk. It therefore proves $R$ is a Banach-space-valued holomorphic function, establishing complex differentiability in the Banach-space norm and gives $R^{\prime }=-R^2$, without assuming a Banach-algebra-valued analytic theorem. In particular every continuous complex linear functional $\varphi$ makes $\varphi (R(\lambda ))$ a scalar holomorphic function.

The Banach algebra spectrum is nonempty. Otherwise $R$ would be defined and analytic everywhere. Its local norm continuity and decay at infinity make it globally bounded. For every continuous linear functional $\varphi$, the scalar entire function $\varphi (R(\lambda ))$ is bounded and tends to zero at infinity, so it is identically zero: the scalar Cauchy estimate for its derivative on circles of radius $r$ is $M/r$, which tends to zero, and the limiting value then fixes its constant as zero. Such functionals separate points by the Hahn-Banach theorem. Explicitly a bounded real functional can be extended from the real span of any nonzero element using its norm, and complexifying it by $\varphi (v)=\ell (v)-i\ell (iv)$ gives a bounded complex linear functional nonzero on that element. Therefore $R$ itself would be zero, contradicting $(\lambda e-x)R=e$. We have proved the Banach algebra spectrum is a nonempty compact subset of the complex plane.

We next establish the spectral radius formula. The power norms are submultiplicative. If some $x^m=0$, their root norms are eventually zero. Otherwise, for any fixed $m$, write $n=qm+r$ with $0\le r<m$ to obtainTaking roots gives ${\mathrm{limsup}}_n\Vert x^n{\Vert }^{1/n}\le \Vert x^m{\Vert }^{1/m}$. Taking the infimum in $m$, and using that every root norm is at least that infimum, proves existence ofFor $|\lambda |>{\rho }_0$, the root test makes $\sum x^j/{\lambda }^{j+1}$ converge, even when $|\lambda |\le \Vert x\Vert$. The same telescoping multiplication proves it is a two-sided inverse. HenceFor the reverse bound fix $R>r_{\sigma }$. If $R_1>\max (R,\Vert x\Vert )$, integrate the uniformly convergent exterior series term by term on $|\lambda |=R_1$, obtainingThe same integral equals the one on $|\lambda |=R$: applying an arbitrary continuous linear functional reduces contour deformation through the Banach algebra spectrum-free annulus to the scalar Cauchy theorem, and separation of points restores equality in $B$. This proves the resolvent Cauchy coefficient formula rather than assuming vector-valued analytic calculus. Using the norm bound for the allowed vector integral givesTaking $n$th roots yields ${\rho }_0\le R$, and letting $R\downarrow r_{\sigma }$ finishes the proof:

For the first example take the commutative algebra of dual numbers $B=\mathbb{C}[\epsilon ]/({\epsilon }^2)$ with norm $\Vert a+b\epsilon \Vert =|a|+|b|$. Completeness follows from finite dimension, and the product estimate makes the norm submultiplicative. For $x=\epsilon$, $x\ne 0$ but $x^2=0$, so $\rho (x)=0$. Directly, $\lambda e-\epsilon$ has inverse ${\lambda }^{-1}e+{\lambda }^{-2}\epsilon$ for $\lambda \ne 0$, and its Banach algebra spectrum is $\{0\}$. For the second example take $B=\mathbb{C}$ with its absolute-value norm and $x=1$. Its Banach algebra spectrum is $\{1\}$, giving $\rho (x)=\Vert x\Vert =1$.

Take a Hankel contour $C_{\eta }$ around the positive real axis, going from infinity on its lower side toward a circle of radius $0<\eta <1$, clockwise around that circle, then out along the upper side. Use $0\le \arg z\le 2\pi$ and defineThis is an entire function of $t$: the small circle stays away from zero and the pole at $-1$, while the exponentially decaying rays converge uniformly for $t$ in compact subsets. Contour deformation shows independence of $\eta$. For $\Re t>-1$, let $\eta \to 0$. The upper ray contributes $F(t)$, the lower ray contributes $-e^{2\pi it}F(t)$, and the small-circle contribution vanishes. Thus the analytic continuation isThe only apparent singularities are at integers. At the nonnegative integers they are removable because the original integral is holomorphic throughout $\Re t>-1$. Hence singularities can occur only at $-1,-2,\dots$.

In fact all of those points are simple poles. At $t=-n-1$ the two rays cancel and the clockwise circle gives $H(-n-1)=-2\pi i{c}_n$, whereSince the denominator has derivative $-2\pi i$ at each integer, the residue at $-n-1$ is $c_n$. The parameter has been called $t$ consistently; the final use of $z$ for this parameter in the printed question is harmless notation.

Set $\omega (k)=k^2-i\alpha k$ and $\nu (k)=i\alpha -k$, so $\omega (\nu )=\omega (k)$. For boundary traces $h_j(s)={\partial }_x^{j}q(0,s)$, defineThe Half-range Fourier transform of the equation is obtained by two integrations by parts:Solving this scalar equation gives the drift-diffusion half-line global relationIn particular the reflected relation is valid on $L$, since $\mathrm{Im}\nu (k)\le 0$ there.

In the region between the real axis and $L$, $\mathrm{Re}\omega \ge 0$. Thus a boundary contribution $e^{-\omega t}{G}_j={\int }_0^t{e}^{-\omega (t-s)}{h}_j(s)ds$ admits an upward contour deformation from the real axis to $L$. Large-contour estimates follow from the spatial exponential and the time integral; at the upper time endpoint one can integrate in time once or use a small contour displacement. The initial-data transform is retained on the real axis. Hence, if $I_C[F]=(2\pi {)}^{-1}{\int }_C{e}^{ikx-\omega (k)t}F(k)dk$, ordinary inversion first givesFor the Dirichlet boundary condition, take $h_0=g_0$. Subtract the null integral of the reflected time-dependent transform, using the inversion from part (a) with $c=-1$. The unknown derivative transform cancels because $i\nu +\alpha =-ik$. This derives the Dirichlet spectral formula for half-line constant drift:Only initial and prescribed boundary data occur in this expression.

For the Neumann boundary condition, take $h_1=g_1$. The reflected global relation reads $e^{\omega t}\hat{q}(\nu ,t)=Q(\nu )-G_1+ik{G}_0$. Solve for $G_0$ and substitute it in the preliminary representation. The leftover term is a null integral with multiplier $(ik+\alpha )/(ik)$, analytic above $L$; its pole at zero is below the contour. Thus the Neumann spectral formula for half-line constant drift isThere is no residue from the origin, since this elimination never moves the contour across it. The given derivative is the inward-coordinate derivative $q_x$, not the outward derivative $-q_x$. The initial compatibility for this case is $q_0^{\prime }(0)=g_1(0)$.

Define the negative-exponential Half-range Fourier transformswhich are holomorphic for $\mathrm{Im}k<0$. Write $g_j(s)={\partial }_x^{j}q(0,s)$ and define the finite-time spectral boundary transformsIntegrate the local relation over $0<x<\infty$, $0<s<t$. The spatial boundary term at infinity vanishes by decay, and the term at $x=0$ is subtracted. Consequently the half-line linear dispersive Stokes global relation isFourier inversion on the real axis gives an initial-data integral plus an integral of $B$. To put the boundary term on a complex contour, setFor $k=u+iv$, $\mathrm{Re}\omega =v(3u^2-v^2-1)$, so $D_{+}$ is the connected upper domain $v>0$, $v^2>3u^2-1$. Orient $\partial D_{+}$ with this domain on the left: from upper-left infinity down to $-1/\sqrt{3}$, along the real segment to $1/\sqrt{3}$, and then to upper-right infinity.

The finite-time spectral boundary transforms are entire functions. In the upper regions outside $D_{+}$, where $\mathrm{Re}\omega \ge 0$, their time-integrated factors satisfy $|e^{-\omega t}{e}^{\omega s}|\le 1$ for $0\le s\le t$. Together with $e^{ikx}$, this permits contour deformation of the two remaining real-axis pieces to the curved boundary of $D_{+}$. The resulting representation isThe contour orientation accounts for the plus sign in this convention. This formula still contains the two unknown derivative traces through $G_1,G_2$.

Use the half-line transform $\hat{u}(k,t)={\int }_0^{\infty }{e}^{-ikx}u(x,t)dx$ and let $g_j(t)={\partial }_x^{j}u(0,t)$. Integrating the heat equation twice by parts givesHere ${\hat{u}}_0(k)=(1+ik{)}^{-2}$ andwith removable singularities at the apparent denominator zeros. Let $D^{+}=\{\mathrm{Im}k>0,\mathrm{Re}{k}^2<0\}$, and orient its boundary from infinity on the $3\pi /4$ ray to zero and then out on the $\pi /4$ ray. Transform inversion and contour deformation giveTo eliminate the unknown boundary derivative, evaluate the transform identity at $-k$: ${\tilde{g}}_1={\hat{u}}_0(-k)+ik{\tilde{g}}_0-e^{k^{2}t}\hat{u}(-k,t)$. The resulting extra contour integral of $e^{ikx}\hat{u}(-k,t)$ vanishes by analyticity and decay in the upper half-plane. This proves that the boxed expression involves only the prescribed data. It is intended for $x,t>0$, with boundary and initial values obtained as limits.

The Malgrange–Ehrenpreis theorem states that every nonzero constant-coefficient linear differential operator on Euclidean space has a distributional fundamental solution. With $D_j={\partial }_{x_j}$ and $P(D)={\sum }_{|\alpha |\le N}{a}_{\alpha }{D}^{\alpha }$, the conclusion is an $E\in {\mathcal{D}}^{\prime }({\mathbb{R}}^d)$ such that $P(D)E={\delta }_0$. If $N=0$, simply take $E=a_0^{-1}{\delta }_0$; hence assume $N\ge 1$.

We construct a Hörmander staircase in frequency space. The highest-degree homogeneous part $P_N$ is not identically zero on real vectors, since a polynomial vanishing on all real vectors has all coefficients zero. After an orthogonal transformation of coordinates, we can arrange that $P_N(e_1)\ne 0$. Thus, for each real ${\xi }^{\prime }\in {\mathbb{R}}^{d-1}$, the polynomialhas degree $N$ in $z$ and the same nonzero leading coefficient $b=i^N{P}_N(e_1)$, independent of ${\xi }^{\prime }$. This constant leading coefficient of a polynomial is what makes a uniform staircase possible.

For a fixed ${\xi }^{\prime }$, factor $p(z,{\xi }^{\prime })=b{\prod }_{\nu =1}^N(z-z_{\nu })$, counting repeated roots of a polynomial. Consider the $N+1$ heights $h_j=3j$, $j=0,\dots ,N$. A given root has imaginary part within distance less than one of at most one of these heights. By the pigeonhole principle, some height avoids every root, and thenThis is finite-height polynomial root avoidance. Root labels need not be chosen continuously or even measurably. Instead define the closed setsContinuity in $s$ extends the inequality from rational $q$ to every real $s$. These Borel sets partition ${\mathbb{R}}^{d-1}$. The Hörmander staircase assigns the horizontal contour $s+i{h}_j$ over $A_j$; its height is bounded by $3N$, and the denominator has the uniform lower bound $|b|$. For $d=1$, the transverse space is a single point and only one horizontal contour is needed.

Use the Fourier transform convention $\hat{\phi }(\zeta )=\int e^{-ix\cdot \zeta }\phi (x)dx$, with inverse factor $(2\pi {)}^{-d}$. Define the candidate fundamental solution of a linear differential operator byFor a test function supported in a fixed compact set $K$, Fourier decay in a bounded complex strip gives, for any integer $M$,Indeed, this is the real Fourier transform of $e^{-h_j{x}_1}\phi (x)$ with reversed frequency, and repeated integration by parts with $(1-\Delta {)}^M$ proves the estimate. Taking $2M>d$ and using the denominator bound proves absolute convergence and a continuity estimate on ${\mathcal{D}}_K$. Thus $E$ is a distribution; no unsupported interpretation of a divergent inverse Fourier transform is being used.

For its distributional derivatives, the formal transpose of a differential operator is $P(-D)$. Sinceapplying $P(D)$ cancels the denominator. For each fixed ${\xi }^{\prime }$, the numerator is an entire function of $z$. The Cauchy integral theorem, applied to a rectangle between the lines $\mathrm{Im}z=0$ and $\mathrm{Im}z=h_j$, givesThe two vertical edges tend to zero by the same bounded-strip decay. The resulting contour deformation is performed separately for each transverse frequency, so discontinuities of the staircase height introduce no additional boundary terms. Absolute convergence permits integration over the partition $A_j$. The Fourier inversion theorem then yieldsThis proves the Malgrange–Ehrenpreis theorem. Undoing the orthogonal change of coordinates gives the fundamental solution of a linear differential operator for the original operator; the Dirac delta distribution is unchanged by that change of coordinates.

For a one-dimensional operator $L=P(D)$ of degree $N\ge 1$, with leading coefficient $a_N$, we may choose the single staircase contour below every pole. Take $\gamma >\max \{\mathrm{Re}\lambda :P(\lambda )=0\}$ and integrate on $\mathrm{Im}\zeta =-\gamma$. The corresponding formula is the Bromwich contour version of the construction above, after putting $\lambda =i\zeta$. Its poles all lie above the frequency contour. For $x<0$, close that contour downwards; the exponential function $e^{ix\zeta }$ decays and there are no enclosed poles. For $x>0$, close upwards, where the exponential function again decays, and apply the residue theorem. The bound $1/P(i\zeta )=O(|\zeta {|}^{-N})$ justifies the large arcs, including $N=1$ by the Jordan lemma away from their endpoints. Thus away from $x=0$ the retarded fundamental solution is $H(x)u(x)$, whereThe sum is over distinct roots of a polynomial, with the residue including the whole multiplicity. For simple characteristic roots of a constant-coefficient differential equation, it reduces to $u(x)={\sum }_{\lambda }{e}^{\lambda x}/P^{\prime }(\lambda )$. Repeated roots give exponential polynomial solutions of a constant-coefficient differential equation through differentiation of $e^{zx}$.

To establish the equality also at the origin, rather than leave a possible point-supported term undecided, verify the distributional jump formula for a Heaviside product. The residue expression is an entire function of $x$, andIts initial derivatives areThese identities follow by integrating over a large circle: for $j\le N-2$ the integrand is $O(|z{|}^{-2})$, whereas for $j=N-1$ its coefficient of $z^{-1}$ is $a_N^{-1}$. The distributional jump formula for a Heaviside product therefore gives $L(Hu)={\delta }_0$. Moreover, since $e^{-\gamma x}Hu$ is a tempered distribution, its Fourier transform satisfiesThere are no real zeros of this polynomial, so its reciprocal is the transform. This is exactly the shifted-contour construction, proving equality there as a distribution as well. HenceThe sign of $a_N$ matters: for the operator in the preceding solution, $a_N=-1$ and the last initial derivative is $-1$.

Use the Fourier transform pair $\hat{q}(k)=\int q(x)e^{ikx}dx$, $q(x)=(2\pi {)}^{-1}\int \hat{q}(k)e^{-ikx}dk$. The point force transforms to $F$. With the dispersion relation $D$ defined above, the sheet equation becomesThe point-force radiation from a fluid-loaded sheet is therefore represented exactly byThe causal contour and outgoing acoustic square-root branch are fixed first with $\mathrm{Im}\omega <0$, then continued to the desired real frequency. This prescription fixes how poles and the branch points are passed.

For the acoustic far field $x=r\cos \theta$, $y=r\sin \theta$, take $0<\theta <\pi$ bounded away from grazing and $k_{0}r\gg 1$. The method of steepest descent saddle point is $k_s=k_0\cos \theta$, with ${\gamma }_s=i{k}_0\sin \theta$. The supplied saddle point rule, including its $\sin \theta$ factor, gives, provided the contour deformation crosses no poles,A convenient simplification, free of division by $T$, isEquivalently, when $T\ne 0$,The expression printed in the PDF is missing sound-speed factors for general dimensional $c_0$. It agrees with this result if $c_0=1$ in fully normalized units; when $c_0$ is retained as an arbitrary sound speed, the numerator needs $c_0^{-3/2}$ and the structural term needs $T/c_0^2$ in the last form. These factors arise respectively from cylindrical spreading, the pressure-density relation, and $k_s=\omega \cos \theta /c_0$.

A direct countercheck is the transparent-sheet limit $m=T=0$. The sheet jump condition then gives ${\hat{p}}_{+}(k,0)=F/2$, so the saddle point rule requiresThe printed expression, interpreted continuously after multiplying out its structural factor, instead gives $F(2c_0{)}^{-1}\sqrt{\omega /(2\pi r)}\sin \theta$ times the same phase. It differs by a factor $c_0^{3/2}$; for example it is eight times too large when $c_0=4$. This limit also verifies the normalization of the corrected density field independently of the sheet's elastic-sheet tension.

To decide about poles, track the roots of $\Delta (k,\omega )$ on the chosen square-root sheet and deform the original causal contour to the steepest descent contour. A root contributes a residue exactly when it lies in the region swept out by that deformation; its sign is fixed by the contour orientation. Branch cuts must be retained throughout this comparison. Which roots are crossed can depend on observation angle, producing a change of the modal contribution when a pole meets the deformation boundary. A saddle point approaching a pole or a grazing endpoint requires an approximation uniform in that limit, rather than the isolated saddle point formula above.

The crossed poles are the free fluid-sheet modes of the preceding solution. Real subsonic roots represent evanescent acoustic surface waves carrying energy along the sheet, with normal decay; complex continuations represent leaky or radiating modes. Their residues must be added to the saddle point sound when the causal contour selects them. The specification “no poles contribute” is therefore a substantive condition on the contour, not permission to ignore zeros of the dispersion relation.

The dispersion symmetry elimination of a boundary trace uses the symmetry of the dispersion relationFor $k\in \overline{D^{+}}$, $\mathrm{Im}\nu (k)\le 0$, so the global relation is valid at $\nu (k)$. Since $i\nu +\alpha =-ik$, it givesSubstitution into the contour integral representation produces an unwanted integral ${\int }_{\partial D^{+}}{e}^{ikx}\hat{u}(\nu (k),t)dk$. Its integrand is analytic in $D^{+}$ and decays on closing the contour upwards; Jordan lemma makes this integral zero for $x>0$. Hence the Fokas method eliminates the unknown normal derivative:All quantities here are determined by the prescribed initial and Dirichlet boundary data up to time $t$.

For verification and for numerical evaluation it is useful to evaluate the spectral contour integrals, giving a half-line drift reflection kernel. Putand define the half-line drift boundary kernelFubini's theorem, the Gaussian Fourier transform and contour deformation give the equivalent causal formulaFor the reflected initial term, ${\hat{u}}_0(i\alpha -k)={\int }_0^{\infty }{e}^{\alpha y+iky}{u}_0(y)dy$ on the contour; the evaluated Gaussian supplies the necessary large-$y$ decay. If $e^{\alpha y}{u}_0(y)$ is not integrable, truncate the initial conditions first, evaluate, and pass to the limit using Gaussian bounds. No extra exponential-decay assumption on the original data is needed for this kernel formula.

The boundary kernel follows particularly simply fromThis calculation independently checks both the sign and the coefficient of the boundary forcing.

Apply Fourier inversion to the global relation on the real axis, extending $u$ by zero to negative $x$. At an interior point $x>0$ this givesThe finite-time spectral boundary transforms are entire functions of $k$. Define the upper decay domainOrient $\partial D^{+}$ from its left infinite end through $i\alpha$ to its right infinite end, so that $D^{+}$ lies to the left. Its finite point is $i\alpha$, not zero. On the contour, $\mathrm{Re}\omega =0$; the factor $e^{ikx}$ decays because $x>0$.

For the boundary term, write $e^{-\omega t}{G}_j={\int }_0^t{e}^{-\omega (t-s)}{g}_j(s)ds$. Between the real axis and $\partial D^{+}$, $\mathrm{Re}\omega \ge 0$, so this factor has no exponential growth. Contour deformation and Jordan lemma, with the usual cutoffs for oscillatory integrals, therefore giveThis is the requested complex-plane representation with the unknown finite-time spectral boundary transform still present. The decay domain must agree with the sign of the drift in the dispersion relation.

Use the principal complex logarithm in the upper half-plane. Deform the interval upward into a half-strip: the vertical side at zero is traversed upward, and the side at one downward. The top side vanishes exponentially; a small indentation at zero contributes $O(r|\log r|)$ and vanishes. Thus the contour deformation gives the exact representation, for $x>0$,Since $\log (is)=\log s+i\pi /2$, scaling $u=xs$ and using the integral defining Euler's constant givesFor the second endpoint, $\log (1+is)=is+s^2/2+O(s^3)$. Watson's lemma, or integration of these powers against $e^{-xs}$ with an exponentially small tail, yieldsIn particular the required first terms areThe two different endpoint scales come from the logarithmic singularity at zero and the simple zero of the amplitude at one.

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.