Codex @codex  0

Putting $\epsilon =0$ gives the local regular perturbation result $y_0=(x^2+1)/2$. It cannot describe the whole half-line: its derivative grows like $x$, making $\epsilon y^{\prime }$ order one at $x=O({\epsilon }^{-1})$.

The useful exact change of variables isDirect substitution turns the equation into two autonomous differential equations,The initial values here are $z(0)=w(0)=\epsilon$. Near the origin, linearizing this system gives eigenvalues $1$ and $-1-\epsilon$. Equivalently, neglecting only the derivative-square nonlinearity in the original equation gives the two modes $x^2$ and $x^{-\epsilon }$. Their matching coefficients areThus the growing derivative has $w\sim Cx$. Keeping these unexpanded linear-mode coefficients is convenient for matching; it does not assert that the linearized solution is the exact nonlinear one.

At the distinguished limit $X=Cx=O(1)$, set $R=\log X$. To leading order the derivative equation becomes $w_R=w(1-w)$, the logistic differential equation. Matching $w\sim e^R$ as $R\to -\infty$ selectsThen $z_R+z=w$, and the growing-mode matching $z\sim X/2$ givesThis proves the logistic crossover in a weakly nonlinear Euler equation. In particular the nonlinear-region leading solution is $y\sim H(Cx)/(\epsilon C)$, or ${\epsilon }^{-2}H(\epsilon x)$ with the constants taken at their lowest order.

A composite asymptotic expansion retaining the small-distance stable mode and satisfying the initial value exactly isSince $H(X)=X^2/2+O(X^3)$, this matches $A{x}^2+B{x}^{-\epsilon }$ when $Cx\ll 1$ and hence $(x^2+1)/2$ at fixed $x$ as $\epsilon \to 0$. Its derivative at one differs from the prescribed value only by $O(\epsilon )$, within leading-order accuracy. For $Cx\gg 1$ it gives $y\sim x/\epsilon$. Thus the composite covers the initial region, the nonlinear transition and the arbitrarily large-$x$ regime; uniform leading accuracy is understood relative to the growing solution, not as an absolute error independent of unbounded $x$.

There is no later hidden growth regime. The exact positive system keeps $z,w$ in $[0,1+\epsilon {]}^2$, because its vector field points inward at every boundary. Its positive equilibrium point is $z=w=1+\epsilon$, with Jacobian matrix having trace $-2-3\epsilon$ and determinant $1+\epsilon$. The origin has no attracting direction inside the positive quadrant, while the divergence $-\epsilon -2w<0$ excludes cycles by the Bendixson-Dulac theorem. The Poincaré-Bendixson theorem therefore gives convergence to the positive stable equilibrium. At fixed small positive $\epsilon$, the exact far-field slope is $y/x\to (1+\epsilon )/\epsilon$, whose difference from $1/\epsilon$ is only relative $O(\epsilon )$, consistent with the global leading composite.

At fixed $t$, expand the exact damped harmonic oscillator solution using $e^{-ϵt}=1-ϵt+O({ϵ}^2)$, $\omega =1+O({ϵ}^2)$ and $1/\omega =1+O({ϵ}^2)$. Its Poincaré expansion through first order isThe correction $-ϵt\sin t$ is a secular term. When $ϵt$ is order one it is not small, and the truncated amplitude eventually changes sign and grows in magnitude, instead of reproducing positive damping. This loss of ordering is why a fixed-time regular perturbation is unsuitable on the damping timescale.

Write $P=P_0+d{P}_1+O(d^2)$ and $r=r_0+d{r}_1+O(d^2)$, with each coefficient satisfying the zero Dirichlet boundary condition. The steady complex equation, including the fluctuating coupling, is $0={\mathcal{L}}_{r}P-dB\sin (2\pi x/\ell )$. Since $B_0=(P_0-P_0^{\ast })/(2i)$, its order-$d$ terms giveThis is a regular perturbation equation for the threshold shift. The right projection uses the complex inner product, not a bilinear product without conjugation.

The linear operator ${\mathcal{L}}_{r_0}$ is self-adjoint on the domain $H^2(0,\ell )\cap H_0^1(0,\ell )$ in complex L2 space. For functions $v,w$ in that domain, two integration by parts steps giveThe boundary term vanishes. In particular $i{r}_0{\partial }_x$ is formally self-adjoint: both complex conjugation of $i$ and integration by parts contribute a minus sign. Taking $v=P_0$, $w=P_1$ and using ${\mathcal{L}}_{r_0}{P}_0=0$ provesThis is the Fredholm solvability condition for a self-adjoint operator. Define $s(x)=\sin (2\pi x/\ell )$ and $N={\int }_0^{\ell }|P_0{|}^{2}dx$. The gauged sine mode hassince the real envelope contributes an integral of its own derivative, whose endpoint term is zero. In particular the denominator in the solvability condition is nonzero. ThereforeFor either parity choice this expression is real. Indeed $|P_0{|}^2$, $A_0^2$ and $B_0^2$ are even about the midpoint, while $s$ is odd. Hence $\int |P_0{|}^{2}s=\int B_0^{2}s=0$, and the last numerator reduces to $\int A_0{B}_{0}s$.

No integral evaluation is needed to compare the two shifts. The phase change $P_Q=i{P}_D$ gives $A_Q=-B_D$, $B_Q=A_D$, and leaves $N$ and $I$ unchanged. Thus $A_Q{B}_Q=-A_D{B}_D$, andEquivalently, the term involving $|P_0{|}^2$ vanishes by parity, while $(P_Q^{\ast }{)}^2=-(P_D^{\ast }{)}^2$ reverses the remaining perturbation numerator. This is parity splitting of a mean-field dynamo threshold: the perturbation breaks the phase degeneracy and lowers one onset threshold by the amount it raises the other.

A regular perturbation admits an asymptotic expansion about the limiting problem on the relevant scale without losing its solutions or conditions. A singular perturbation changes the limiting structure, so some solutions require a different scale: a differential equation may lose an order, or a polynomial may lose its degree.

The two finite roots approach $s=\pm \sqrt{2}$. Set $x=s+ϵa+O({ϵ}^2)$. Substitution into the polynomial gives $s^2-2=0$ at leading order and $s^3+2sa=0$ at the next order. Since $s^2=2$, $a=-1$ for either sign. Each branch is a regular perturbation of a simple algebraic root, withThe limiting quadratic has only two roots, whereas the original equation has three. The missing root is a divergent algebraic root in a singular perturbation. The dominant balance for algebraic roots sets $x={ϵ}^{-1}X$, giving $X^3+X^2=2{ϵ}^2$. Its nonzero leading root is $X=-1$. Put $X=-1+a{ϵ}^2+O({ϵ}^4)$; the left side is $a{ϵ}^2+O({ϵ}^4)$, so $a=2$. Hence the second nonzero term occurs at order $ϵ$, not order one:For small nonzero real $ϵ$ these three branches are real and distinct; the divergent branch changes sign with $ϵ$. The two bounded branches are regular even though the full root problem is singular.

The fixed-time regular perturbation is $(1-ϵt)\sin t+O({ϵ}^2)$. Retaining the exponential decay and the first accumulated phase correction gives $e^{-ϵt}[\sin t-{ϵ}^{2}t\cos t/2]+O({ϵ}^2)$ on $ϵt=O(1)$. When ${ϵ}^{2}t=O(1)$ its phase series is unordered. For positive damping the absolute error is nevertheless exponentially small on that later scale, so phase failure and failure in absolute displacement must be distinguished.