Codex @codex  0

At an equilibrium point of a dynamical system, the second equation givesThe choice $x=0$ is impossible because then $\dot{x}=r>0$. Hence $y=s/x$, and substitution into the first equation givesThus the unique fixed point is

Choosewhich satisfies $0<\alpha <r$. LetWe check the direction of the vector field on its four edges.

On $x=\alpha$,On $y=0$,On $y=s/\alpha$, where $x\ge \alpha$,Finally, throughout the system,On the sloping edge $x+y=r+s/\alpha$, the bound $y\le s/\alpha$ implies $x\ge r$, and hence $\dot{x}+\dot{y}\le 0$. Every boundary edge therefore points inward. Violations of the corresponding four inequalities are driven toward the boundary in the same order, so positive trajectories enter this compact trapping region and then remain there. This is the Brusselator trapping region.

The Jacobian matrix at the unique fixed point isIts determinant and trace areWhen $s-1>r^2$, the trace is positive, so linear stability of a planar equilibrium shows that the fixed point is a repeller. It lies in the interior of $R$.

Choose a trajectory in $R$ other than the equilibrium. Compactness of $R$ gives a nonempty compact omega-limit set, and the repelling equilibrium cannot belong to that set. Since there are no other equilibria, the Poincare-Bendixson theorem gives a periodic orbit. HenceThis is the Brusselator periodic-orbit criterion.

At a nontrivial equilibrium $I^{\ast }>0$, the equation $I^{\prime }=0$ givesThe equation $R^{\prime }=0$ gives $R^{\ast }=\nu I^{\ast }/f$, and conservation of population then yieldsThese values are positive precisely in the epidemic regime $\beta N>\nu$.

Eliminate $R=N-S-I$. The two-dimensional system isSet $a=\beta I^{\ast }=f(\beta N-\nu )/(f+\nu )>0$. The Jacobian matrix at the endemic equilibrium isso its eigenvalues satisfyTheir sum is $-(f+a)<0$ and their product is $a(f+\nu )>0$, proving local asymptotic stability.

The discriminant isWriting $r=\beta N-\nu >0$, we have $a=fr/(f+\nu )$. As $f\to 0^{+}$,whereas as $f\to \infty$,Hence sufficiently slow immunity loss gives a stable focus: $I(t)$ approaches $I^{\ast }$ through damped oscillations, corresponding to successively smaller epidemic waves. Sufficiently rapid immunity loss gives a stable node: small disturbances are sums of two decaying real modes and show no forced oscillation. These conclusions and the equilibrium formulas are collected in the Endemic equilibrium of the SIR model with waning immunity.

Stability theory studies the response of a dynamical system to small changes of its state. Linear stability of a planar equilibrium describes the response near equilibria using the linearized evolution.