Lyapunov function and an open disk inside the basin of $(0,0)$

Let's first start off with LaSalle's Invariance Principle. It says:

Let $X^*$ be an equilibrium point for $X'= F(X)$ and let $L: U \to \mathbb{R}$ be a Liapunov function for $X^*$, where $U$ is an open set containing $X^*$. Let $P \subseteq U$ be a neighborhood of $X^*$ that is closed and bounded. Suppose that $P$ is positively invariant, and that there is no entire solution in $P\setminus\{X^*\}$ on which $L$ is constant. Then $X^*$ is asymptotically stable and $P$ is contained in the basin of attraction of $X^*$.

So basically we have our system $X' = F(X)$ and we have an equilibrium $X^*$ with a closed neighborhood $P$ of that equilibrium. So we just have to check 3 things:

1) $P$ is closed and bounded

2) $P$ is positively invariant

3) $L$ is never constant along any non-equilibrium solution which resides in $P$

If these are satisfied, then we have 2 conclusions:

A) The equilibrium is asymptotically stable

B) The region $P$ is in the basin of attraction

So what they tried to show is that the circle of radius 2 centered at origin is closed and bounded, positively invariant, and $L$ isn't constant along any non-equilibrium solution in the circle.

I'm not entirely convinced of this because the Liapunov function may no longer be strict on the boundary, but I am more confident in saying that the circle of radius 1 centered at the origin is in the basin.

The circle of radius 1 centered at the origin is closed and bounded and since the Liapunov function is strict in this circle, $L$ will not be constant along any non-equilibrium solution curve. Furthermore, our Liapunov function corresponds to the norm of our solution squared (the square of the distance from the origin to a point on the solution). Thus if a solution curve starts inside this circle and later left, it would correspond to an increase in the the distance of the solution from the origin. However, our Liapunov function is strict, so it is strictly decreasing. Therefore, the circle is positively invariant. Now we can apply LaSalle's to conclude the circle of radius 1 centered at the origin is in the basin.


  1. I beleive the simplest way to show that $\frac{\cos\theta\sin\theta(\cos\theta+\sin\theta)}{1+\cos^2\theta}>-\frac12$ is to plot the left part of the inequality: $\frac{\cos\theta\sin\theta(\cos\theta+\sin\theta)}{1+\cos^2\theta}$ They didn't ask to prove the inequality, right?

  2. As one can see, the minimum value of $\frac{\cos\theta\sin\theta(\cos\theta+\sin\theta)}{1+\cos^2\theta}$ is not exactly $-1/2$ (it is approximately $-0.4941$), so it is not quite correct to say that the radius $\delta=2$ is as much as possible.

  3. In fact, the LaSalle's Invariance Principle is not needed here. It is intended to use for non strictly negative $\dot L$. Since $\dot L<0$ inside the disk, we can conclude that $L$ is decreasing along the trajectories and replicate the proof of the Lyapunov's theorem on asymptotic stability.

  4. The polar coordinates are not needed too. The general way to obtain the largest possible (for a fixed Lyapunov function) subset of the domain of attraction looks as follows.

Suppose we have the Lyapunov function $L(x)$, its derivative is $\dot L(x)$. Consider the set $$ S= \left\{ x: \dot L(x)=0\right\}\setminus \{0\}. $$ Suppose that we can solve the minimization problem $$\tag{1} C=\min_{x\in S} L(x). $$ (For this system, in particular, $C\approx 4.0955$). The level set $$ \Omega_C = \left\{ x: L(x)<C \right\} $$ is the largest possible positive invariant set contained in the domain of attraction of the origin that can be obtained for the fixed Lyapunov functon. This fact can be demonstrated by this picture: the invariant set The point of intersection of the green and red curves is the minimum point of (1). Indeed, since $L(x)<C$ for any point of the level set $\Omega_C$, this set does not contain points such that $\dot L\ge 0$; at the intersection point the value of $\dot L$ is zero, so we can not say that any level set that is greater than $\Omega_C$ is a positively invariant set (as it contains the points where $\dot L\ge0$). This metod also applies in higher dimensions, not only in 2D