Past exam of the further mathematical biology course of the University of Oxford 2023 2 b Solution 2026-10-25
There is a consistency issue in the printed model. If is imposed while retaining its evolution equation, then forces . Thus the full two-equation system has no nonconstant travelling wave with exactly constant . Its constant solution exists only when . The intended wave calculation must treat the recovery variable as externally frozen, leaving the scalar reaction-diffusion equationFor this frozen-recovery reduction, letThe cubic polynomial has a local minimum at and a local maximum at . Hence the conditiongives three distinct real roots of , withThe outer roots are stable for the scalar spatially uniform dynamics and the middle root is unstable. Seek a travelling wave , where . The travelling-wave reduction of a reaction-diffusion system givesLet and chooseIt has , andSubstitution leavesConsequently the frozen-recovery cubic front isReflection reverses the front and its speed. A nonzero speed requires ; equality gives a stationary front and is the Maxwell balance for a scalar reaction-diffusion front. Since the roots sum to , this stationary value is . For the simple admissible choice , the roots are andThis particular front moves when . At , other currents inside the displayed open interval give moving fronts.
Let denote the increasing-front speed from the preceding solution. Reversing that heteroclinic reverses its speed, so a decreasing back has speed . A rigidly travelling pulse has one speed throughout, requiringThis is the correct use of the integral formula: compare the oppositely oriented transitions, including their signs, rather than equating two increasing-front speeds.
For completeness, write , . The explicit increasing frozen-recovery cubic front isSubstitution, using , verifies the equation. Since , speed matching givesFor this selects the appropriate positive back level within the cubic's three-root range. The high branch connects to slowly, and the low branch supplies recovery after .
The stated positive-speed, rising-front orientation imposes at leading order. The printed range by itself does not guarantee that pulse construction: at the leading front is stationary, while for the calculated increasing-front speed is negative. Existence and stability of a full finite- pulse also require the slow connections; speed matching is its leading necessary construction, rather than a proof for every parameter in the printed range.