Past exam of the differential equations 1 course of the University of Oxford / 2021 / 3 / a

The equilibrium points satisfy $y=\beta$ and $x(x^2-1)=0$, so they areThe Jacobian matrix and its eigenvalue equation areAt $x=\pm 1$ the eigenvalues are $\pm \sqrt{2}$, so both points are saddle equilibria. Their stable and unstable tangent directions have slopes $-\sqrt{2}$ and $\sqrt{2}$, respectively.

At $x=0$ the eigenvalues are $\pm i$, giving a centre for the linearized equation. Purely imaginary eigenvalues alone do not classify the nonlinear system. The first integral established in part (ii) proves that $(0,\beta )$ is in fact a centre equilibrium. It is clockwise: immediately to its right, with $0<x<1$ and $y=\beta$, the vector field points downwards. Thus the final nonlinear classification is

Where $y\ne \beta$, taking the quotient of the two ordinary differential equations givesChoose the first integralThe chain rule checks conservation even at points excluded in taking that quotient:Thus each trajectory is contained in a connected component of a level set of $h$. At $(0,\beta )$ the gradient vanishes and the Hessian matrix is the identity, so $h$ has a strict local minimum there. To make the closed-trajectory argument explicit, put $z=y-\beta$ and $0<E<1/4$. The bounded component of $h=E$ isThe two arcs join at $(\pm r_E,0)$ to form a simple closed curve that is a regular level set. The gradient cannot vanish on it: the only critical point of a dynamical system in $|x|<1$ is the centre, whose level is $0$.

The vector field is tangent and nowhere zero on this compact set, so its speed has a positive minimum. It therefore traverses the curve in finite time and returns to its initial point, giving a periodic orbit. These ovals shrink to the centre as $E\downarrow 0$. This is the planar first integral criterion for a centre and proves the required nonlinear assertion; it also gives Lyapunov stability without asymptotic stability.

The nullclines areOn $y=\beta$, away from the equilibrium points, the vector field is vertical. Its crossing direction is downwards for $x<-1$, upwards for $-1<x<0$, downwards for $0<x<1$, and upwards for $x>1$, because $\dot{y}=x(x^2-1)$. On each of $x=-1,0,1$, the vector field is horizontal: it points right when $y>\beta$ and left when $y<\beta$.

Throughout each complementary region, $\dot{x}$ has the sign of $y-\beta$, and $\dot{y}$ has the four signs just listed. These signs orient the level sets from part (ii). In particular, the inner periodic orbits are clockwise. The saddle level is $h=1/4$, whose separatrices satisfyFor $-1<x<1$, the upper arc $y-\beta =(1-x^2)/\sqrt{2}$ goes from $(-1,\beta )$ to $(1,\beta )$; the lower arc goes back from $(1,\beta )$ to $(-1,\beta )$. These two heteroclinic orbits form a closed boundary around the family of periodic orbits. The remaining branches of the same saddle level extend outwards. The sketch shows the surrounding nonclosed trajectories as well as the nullclines and equilibrium points.

Set $z=y-\beta$. Then the autonomous planar system becomeswhich has no $\beta$ dependence. ThereforeMore precisely, replacing $\beta$ by $\beta +\delta$ sends every trajectory $(x(t),y(t))$ to $(x(t),y(t)+\delta )$, with the same time parametrization. The equilibrium points, nullclines, periodic orbits and separatrices all undergo this same translation. Their types and orientations are independent of $\beta$, so varying $\beta$ produces no bifurcation.