Lyapunov–Malkin theorem

In the system of differential equations,

${\displaystyle {\dot {x}}=Ax+X(x,y),\quad {\dot {y}}=Y(x,y)}$ 

where ${\displaystyle x\in \mathbb {R} ^{m}}$  and ${\displaystyle y\in \mathbb {R} ^{n}}$  are components of the system state, ${\displaystyle A\in \mathbb {R} ^{m\times m}}$  is a matrix that represents the linear dynamics of ${\displaystyle x}$ , and ${\displaystyle X:\mathbb {R} ^{m}\times \mathbb {R} ^{n}\to \mathbb {R} ^{m}}$  and ${\displaystyle Y:\mathbb {R} ^{m}\times \mathbb {R} ^{n}\to \mathbb {R} ^{n}}$  represent higher-order nonlinear terms. If all eigenvalues of the matrix ${\displaystyle A}$  have negative real parts, and X(x, y), Y(x, y) vanish when x = 0, then the solution x = 0, y = 0 of this system is stable with respect to (x, y) and asymptotically stable with respect to x. If a solution (x(t), y(t)) is close enough to the solution x = 0, y = 0, then

${\displaystyle \lim _{t\to \infty }x(t)=0,\quad \lim _{t\to \infty }y(t)=c.}$ 

Consider the vector field given by

${\displaystyle {\dot {x}}=-x+x^{2}y,\quad {\dot {y}}=xy^{2}}$ 

In this case, A = -1 and X(0, y) = Y(0, y) = 0 for all y, so this system satisfy the hypothesis of Lyapunov-Malkin theorem.

The figure below shows a plot of this vector field along with some trajectories that pass near (0,0). As expected by the theorem, it can be seen that trajectories in the neighborhood of (0,0) converges to a point in the form (0,c).