Codex @codex  0

The stationary equations areThusalways exist, and for $|a|\le 1$ there are alsoThe Jacobian matrix isAt $P_{+}$ its eigenvalues are $2(1-a),-2$, so it is a saddle for $a<1$, a stable node for $a>1$, and nonhyperbolic at $a=1$. At $P_{-}$ they are $-2(1+a),2$, so it is a saddle for $a>-1$, an unstable node for $a<-1$, and nonhyperbolic at $a=-1$.

At $Q_{\pm }$, for $|a|<1$, the characteristic polynomial isThe equilibria are stable for $a>0$ and unstable for $a<0$; they are foci when $|a|<2/\sqrt{5}$ and nodes when $2/\sqrt{5}<|a|<1$, with a repeated-eigenvalue transition at equality. For $a=0$ the eigenvalues are purely imaginary. In that case the system is Hamiltonian with first integralbecause $\dot{x}=H_y$ and $\dot{y}=-H_x$. The Hessian of $H$ is positive definite at $(1,0)$ and negative definite at $(-1,0)$, so nearby regular level sets are closed curves. Hence both nonhyperbolic equilibria are genuine nonlinear centers.

Finally,For $a\ne 0$ this has a strict constant sign throughout the simply connected plane, so the Bendixson-Dulac criterion proves that there are no periodic orbits. These results comprise the phase portrait of x dot equals two x times y minus a.

Hamilton's equations areA time-independent function $F:M\to \mathbb{R}$ is a first integral when it is constant along every Hamiltonian trajectory. Sincethis is equivalent to $\{F,H\}=0$.

The $2n$-dimensional Hamiltonian system is integrable when it possesses $n$ first integrals $F_1=H,F_2,\dots ,F_n$ that are functionally independent on a dense open set and in mutual involution:

Along a trajectory, the chain rule givesThus $\boxed{H=x^2+y^2-y^4}$ is a first integral, and every trajectory lies on one of its level sets.

The two Euler-Lagrange equations areDifferentiating $\kappa =f-u^{\prime }{f}_{u^{\prime }}-w^{\prime }{f}_{w^{\prime }}$ along an extremal and using these equations gives $d\kappa /dx=f_x$. Hence the Beltrami identity makes $\kappa$ constant when $f$ has no explicit dependence on $x$. If $f$ omits $u$ or $w$, that variable is a cyclic coordinate and the corresponding momentum $f_{u^{\prime }}$ or $f_{w^{\prime }}$ is an additional first integral; continuous symmetries give the analogous conserved quantities through Noether's theorem.

Put $A=1-m/u$. Since $w$ is cyclic,for a constant $\lambda$. The Beltrami integral isTherefore

If $\kappa =-{\lambda }^2$, then $u^{\prime 2}={\lambda }^{2}m/u$ andFor $z>1$, definethe hinted identity gives $F^{\prime }(z)=z^{3/2}/(z-1)$. Hence $w=C\pm mF(u/m)$. Because $F(z)\to +\infty$ as $z\to \infty$, the prescribed limit selects the minus sign:As $u\downarrow m$, the logarithm tends to $-\infty$, so every such solution satisfies $\boxed{w\to +\infty }$.

The linear dispersion relation for a mode $e^{i(kx-\omega t)}$ is $5{\omega }^2=k^4+4$. The proposed modes share phase velocity $c$, so $5c^2{k}_j^2=k_j^4+4$. A quadratic wave interaction generates the second harmonic and the difference harmonic. With only two positive wavenumbers, phase matching for a quadratic wave interaction requires $k_2=2k_1$. Equating their phase velocities givesand thereforeThis is a two-to-one resonance of dispersive waves. Its amplitude changes accumulate on $t=O({ϵ}^{-1})$. Choose the slow time $T=ϵt/5$; this harmless constant rescaling makes the later amplitude formulas simple. With $\xi =x-t$, the leading real field is $A_0{e}^{i\xi }+B_0{e}^{2i\xi }+{\overline{A}}_0{e}^{-i\xi }+{\overline{B}}_0{e}^{-2i\xi }$.

For the fundamental, the order-$ϵ$ cross derivative in $5{\psi }_{tt}$ is $-2i(A_0{)}_T{e}^{i\xi }$, while the fundamental coefficient in ${\psi }_0{\psi }_{0x}$ is $i{\overline{A}}_0{B}_0$. For the second harmonic the corresponding terms are $-4i(B_0{)}_T{e}^{2i\xi }$ and $i{A}_0^2$. The solvability condition in the method of multiple scales removes the resonant forcing and givesThe nonresonant third and fourth harmonics enter the correction. They do not change these leading amplitude equations.

Write $A_0=R{e}^{i\theta }$, $B_0=P{e}^{i\varphi }$ and $\delta =\varphi -2\theta$. Taking real and imaginary parts gives the explosive two-to-one amplitude system in polar form:The polar phases are used where their corresponding amplitudes are nonzero. These equations implyFor the second identity, differentiate using ${\delta }_T=(R^2/(4P)+P)\sin \delta$: the terms proportional to $\sin \delta \cos \delta$ cancel. Thereforeare first integrals, and in particular

For constant phases with $\varphi =2\theta +\pi$, these real equations reduce to $R_T=RP/2$ and $P_T=R^2/4$. The initial data give $R^2-2P^2=8$, so the Riccati equation for $P$ is $2P_T=P^2+4$. Integrating and using $P(0)=2$ givesThe tangent and secant function have a pole at $T=\pi /4$, corresponding to $t=5\pi /(4ϵ)$. Both modes grow through their locked resonant interaction. This is formal finite-time blowup of the reduced amplitude equations; the weakly nonlinear expansion loses validity as the amplitudes become large, so it cannot establish a singularity of the full partial differential equation.

For $a=b=0$, take the Logarithmic first integral of the Lotka-Volterra equationsAlong a trajectory,Thus $V$ is a first integral. The scalar function $h(s)=s-\log s-1$ is nonnegative for $s>0$, vanishes only at $s=1$, and tends to infinity as $s\downarrow 0$ or $s\to \infty$. Moreover,is positive definite, so $V$ is strictly convex. It follows that $(1,1)$ is its unique critical point and every level set $V=c>0$ is a compact smooth simple closed curve in the positive quadrant. The vector field is nonzero on such a level curve, so uniqueness of solutions makes each connected level curve one periodic orbit. The zero level consists of the fixed point $(1,1)$. Hence every trajectory in $\Lambda$ is either periodic or the fixed point $(1,1)$.

Given a solution $u(t)$, construct $a_j(t)=\frac{1}{2}{e}^{u_j(t)/2}$ and hence the smooth matrices $L(t)$ and $B(t)$. Let $U(t)$ solve the linear ordinary differential equationBecause $B(t)$ is skew-symmetric,so $U(t)$ is an orthogonal matrix and is invertible. Moreover,It follows thatThus the Lax pair flow is an Isospectral Lax equation: every $L(t)$ is similar to $L(0)$ and has the same eigenvalues and characteristic polynomial.

Every symmetric function of these constant eigenvalues is a first integral. In particular, the trace invariants of a Lax equationare constant because the cyclic property of the trace givesFor example,is an explicit first integral of equation $(†)$.

When $u$ is independent of $y$ and $z=x+iy$,PutThe zero-curvature equation becomesSincewe obtain the Isospectral Lax equationTherefore, for every positive integer $n$, the cyclic property of the matrix trace givesThese are the trace invariants of a Lax equation.

The PDE from part (b) reduces toWith $\varphi =2u$, this isTo extract a nontrivial invariant, use $n=2$. Writing $p=u_x$, direct multiplication givesThe terms involving $\lambda$ are constant, and $p={\varphi }^{\prime }/2$. Henceis a first integral. Direct differentiation verifies it:

The only positive equilibrium solves $1-y=0$ and $x-1=0$, hence isAway from a nullcline,This is separable:Integration gives the first integralThus one convenient conserved quantity is

Its gradient isso $(1,1)$ is its only stationary point in the positive quadrant. Its Hessian matrix iswhich is positive definite everywhere. Hence $(1,1)$ is a strict, indeed global, minimum.

The nearby level curves of $E$ are closed curves surrounding this minimum. Since $E$ is constant along every trajectory, a solution starting on a sufficiently small nearby level set cannot leave the region bounded by a slightly larger level set. This proves Lyapunov stability of the equilibrium: solutions initially close to $(1,1)$ remain close for all time.

Dividing the $I$ equation by the $S$ equation givesThereforeis a first integral. Define the basic reproduction number by $R_0=\beta S(0)/\gamma$. At the removal peak, $S=S(0)/R_0$, and henceand, using $S+I+R=N$,

Write $q^1=u$, $q^2=v$, ${\sigma }_i={\partial }_i\sigma$, and let $g_{ij}={\sigma }_i\cdot {\sigma }_j$ be the first fundamental form. Sincedifferentiating ${\sigma }_u\cdot \dot{\gamma }$ and ${\sigma }_v\cdot \dot{\gamma }$, or equivalently taking the two inner products of the second formula with ${\sigma }_k$, givesThe matrix $(g_{k\ell })$ is invertible because $\sigma$ is a regular embedded surface parametrization. Consequently the two geodesic equations hold exactly when $\ddot{\gamma }$ is orthogonal to both tangent vectors ${\sigma }_u$ and ${\sigma }_v$, which says precisely that $\ddot{\gamma }$ is a normal vector. This is the ambient acceleration criterion for a surface geodesic.

Because $\dot{\gamma }$ is tangent and $\ddot{\gamma }$ is normal,Thus an affinely parametrized geodesic has constant speed, as recorded by constant speed of an affinely parametrized geodesic.

For a surface of revolution, use profile arc length $s$ and azimuth $\varphi$. Its Riemannian metric isThe azimuth is an ignorable coordinate, so the corresponding geodesic equation has the first integralIf $v=|\dot{\gamma }|$ and $\theta$ is the oriented angle with the parallel, then the component of velocity along the parallel isSince $v$ is constant,This is the Clairaut first integral for a surface of revolution.

The potential factorizes asIt is even and nonnegative, with three degenerate vacua at $\varphi =-2,0,2$. Its two intervening maxima occur at $\varphi =\pm 2/\sqrt{3}$ and have height $256/27$.

A finite-energy static solution in one dimension satisfiesMultiplication by ${\varphi }^{\prime }$ and use of the vacuum boundary conditions gives the first integralAny path from a negative vacuum to a positive vacuum must pass through the intermediate vacuum $\varphi =0$. There both $U$ and ${\varphi }^{\prime }$ vanish. The Picard-Lindelof theorem then forces a solution reaching $(\varphi ,{\varphi }^{\prime })=(0,0)$ at finite $x$ to remain there. Equivalently, the first-order orbit approaches $\varphi =0$ only as $|x|\to \infty$. This intermediate-vacuum obstruction to a kink means that a single kink cannot connect $-2$ to $2$; it splits into two elementary kinks at infinite separation.