Past exam of the topics in applied mathematics course of the University of Oxford 2008 B 7 Solution 2026-10-31
For , differentiate and use the two differential equations:Substituting , gives the required equations in polar coordinates. The origin is handled in Cartesian coordinates, where it is an equilibrium point; an angular coordinate is not defined there.
A Lyapunov function near an equilibrium point is a continuously differentiable positive-definite scalar function whose derivative along solutions is nonpositive. This implies local Lyapunov stability under the usual comparison bounds; if its derivative is negative definite, the Second Lyapunov theorem gives local asymptotic stability. Here takeIf , continuity gives a small ball and such that inside it. Consequently . A sufficiently small initial sublevel set stays inside the ball, and . The origin is asymptotically stable, in fact locally exponentially stable. If , continuity instead gives in a small ball. Every nonzero initial point there exits that ball in forward time, proving instability of the origin. These inequalities hold along all local solutions; continuity of alone need not supply uniqueness for a general vector field.
If , a nonzero constant radius satisfying gives a periodic orbit of period , since . Linearizing the scalar radial equation there givesThus the periodic orbit is orbitally asymptotically stable for and unstable for . The direction along the periodic orbit is neutral, as expected for an autonomous periodic orbit; orbital asymptotic stability means attraction in a dynamical system to the orbit with a possible phase shift.
For the specified radial response,The origin is therefore stable for and unstable for . The nonzero-radius equation has roots of an equationretaining only positive radii. For , the square root exceeds one, so only exists and : the periodic orbit is unstable. For , both radii are positive; gives and a stable periodic orbit, while remains unstable.
As increases through , the stable small-radius periodic orbit emerges from the origin while the origin loses Lyapunov stability. This is the supercritical Hopf-type transition with a nonsmooth radial response requested. It follows directly from the explicit radial dynamics. There is a regularity qualification: near zero, so the vector field in Cartesian coordinates contains and terms. It is continuously differentiable but not generally twice differentiable at the origin, and the periodic orbit radius grows as rather than the square-root law of a generic smooth Hopf normal form. One should therefore use the explicit computation here rather than invoke a smooth Hopf bifurcation theorem with unmet hypotheses. At itself the origin is still weakly attracting, because for small positive .