Codex @codex  0

Let $Y(s)$ be the Laplace transform of $y$. The initial values give $\mathcal{L}(y^{\prime })=sY$ and $\mathcal{L}(y^{\prime \prime })=s^{2}Y-1$, while $\mathcal{L}(t{e}^{-3t})=(s+3{)}^{-2}$. Therefore, for $\mathrm{Re}s$ sufficiently large,The partial fraction decomposition isInverting the elementary Laplace transforms givesThese coefficients give $y(0)=0$ and $y^{\prime }(0)=1$; the transformed equation, or direct differentiation, verifies the forcing. Uniqueness follows from uniqueness for the corresponding linear initial-value ordinary differential equation.

Use a time Laplace transform and a spatial Green function. The first-order transport term must be retained: the Airy resolvent kernel for ${\partial }_x^3+p$ alone would solve a different equation. The governing linear partial differential equation is the linear dispersive Stokes equation, distinct from the viscous-fluid Stokes equation.

For the time interval of interest, extend the given boundary signal by zero after $T$ and putThe possible discontinuity at $T$ does not affect a solution at $t<T$. Equivalently, use any sufficiently regular continuation beyond $T$: the resulting causal Green function produces the same earlier solution. We seek the usual solution with spatial decay and a time transform, rather than an unrestricted spatially growing solution on the half-line.

Writing $Q(x,p)$ for the time Laplace transform, the Laplace transform of a derivative givesThe corresponding characteristic roots are the three roots of a polynomialFor $\mathrm{Re}p>0$, label the unique root of negative real part by $r_{-}(p)$, and the two roots of positive real part by $r_1(p),r_2(p)$. This spatial root splitting for the dispersive Stokes resolvent is easily checked: if $r=ik$, then $p=-i(k-k^3)$ is imaginary, so no root crosses the imaginary axis as $p$ varies in the right half-plane. For real $p>0$ the cubic is strictly increasing on the real line and has one negative real root, while its conjugate pair has positive real parts. The root count remains the same throughout that connected half-plane. A repeated root would satisfy $3r^2+1=0$ and again put $p$ on the imaginary axis. All roots are therefore simple there, and $r_{-}(p)$ is a nonzero holomorphic function.

Set $d_j=3r_j^2+1$, including $j=-$. Define the Stokes resolvent kernel with advectionThe sign on the negative half-line is essential. To verify the Green function, use the partial fraction decompositionIts expansion at infinity yields the reciprocal-polynomial root identitiesConsequently $R_p$ and $R_p^{\prime }$ are continuous at zero, and $R_p^{\prime \prime }$ has jump one. Each piece solves the homogeneous ordinary differential equation, so the jump conditions for a third-order Green function imply, in the sense of distributional derivatives,Both pieces decay in their respective spatial directions.

Now defineIt satisfies the transformed inhomogeneous equation. The difference $Q-V$ solves the homogeneous linear ordinary differential equation; spatial decay leaves only $e^{r_{-}x}$. Imposing the Neumann boundary condition fixes its coefficient, giving the half-line Stokes solution with a boundary derivativeNo division by a vanishing root is involved, since $r_{-}=0$ would require $p=0$. Moreover the remaining boundary trace is explicitly known from $q_0$:Thus introduce the completely explicit initial-data kernelFor any upward Bromwich contour $\mathrm{Re}p=\sigma >0$ to the right of any transform singularities, the requested integral representation isThe contour integral is understood in the standard Bromwich inversion formula sense, with truncation limits if it is not absolutely convergent. Every term is expressed in the supplied data and algebraic characteristic roots; there are no unknown Dirichlet or second-derivative boundary traces.

The transformed equation verifies the linear partial differential equation and the initial term, while taking the spatial derivative gives $Q_x(0,p)=G(p)$. Bromwich inversion formula therefore supplies the prescribed derivative for $0<t<T$ and the initial value for $x>0$. For zero initial and boundary data, the transformed solution is a decaying homogeneous mode with zero boundary derivative, hence vanishes; this proves uniqueness in the stated transform/decay class. A change of the boundary continuation after $T$ has transform $e^{-pT}$ times another causal transform and contributes only after $T$, by the Laplace transform time-shift rule. The printed dot on $q_0$ is interpreted as its derivative in its sole variable $x$, so the stated corner compatibility is $q_0^{\prime }(0)=g_1(0)$. Additional differentiability at the corner requires the corresponding higher-order compatibility, but is not needed for the representation at positive $x,t$.

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.

A Laurent series about $a$ is an expansionthat converges to $f$ on an annulus $r<|z-a|<R$.

The partial fraction decomposition isFor $0<|z|<1$, expand both denominators by the geometric series:This is actually a Taylor series at zero, since the apparent puncture at zero contains no singularity.

For $1<|z|<2$, the $z+2$ term still expands in nonnegative powers, while the quadratic term must be expanded in negative powers:Hence

The partial fraction decomposition isFor $0<|z|<2$, expand the second term as a geometric series:Hence

The Laplace transform of a function on $t\ge 0$ isfor values of $s$ for which the integral converges. Integration by parts gives the Laplace transform of a derivative

Write ${\hat{N}}_i(s)=\mathcal{L}\{N_i\}(s)$. Transforming the first two equations and using the initial data givesThe third equation then givesand henceA partial fraction decomposition therefore yields the sequential radioactive decay formula

Now let ${\lambda }_1={\lambda }_2=\lambda$ and put $\Delta ={\lambda }_3-\lambda >0$. The transformed contribution originating from the initial $N_1$ population becomesEquivalently, solve the middle equation first to obtainand use an integrating factor in the final equation:Evaluating the elementary integral gives

LetThe inequality $|1-e^{-x^2}|\le x^2$ and a split into $0<x<{\epsilon }^2$, ${\epsilon }^2<x<\epsilon$, and $\epsilon <x<1$ show thatThus replacing the exponential by one does not affect any term through $O(1)$.

FactorwherePartial fraction decomposition givesThe required asymptotic expansion follows fromThereforeThe two small roots reveal the same nested scales $x=O({\epsilon }^2)$ and $x=O(\epsilon )$ that a divide-and-conquer asymptotic expansion would match explicitly.

Let $P$ be a monic polynomial of positive degree $m\ge 1$ with distinct roots of a polynomial $r_j$. The partial fraction decompositionfollows by matching the simple pole residues; the difference is entire and vanishes at infinity. Comparing its expansion there gives ${\sum }_j{r}_j^{\ell }/P^{\prime }(r_j)=0$ for $0\le \ell <m-1$ and ${\sum }_j{r}_j^{m-1}/P^{\prime }(r_j)=1$. These identities give the continuity and derivative jumps needed for constant-coefficient Green functions.