Stability Concepts and Definitions

What does it mean for a pencil balanced on its tip to be "unstable"? What does it mean for a marble resting at the bottom of a bowl to be "stable"? Both are equilibria of some dynamical system. Yet they behave utterly differently when nudged. This article makes that intuitive difference completely precise. It uses only distances and inequalities. No linearization, no eigenvalues, no calculus of variations is required at the outset. The concepts of Lyapunov stability, asymptotic stability, and instability form a vocabulary. Every later theorem about long-term behavior of dynamical systems is phrased in it. Before we can prove that a system settles into equilibrium, we need sharp definitions. The same holds if it spirals away, or oscillates forever nearby. We need a definition sharp enough to be tested against examples. It must also be robust enough to survive generalization. This includes infinite dimensions, non-autonomous systems, and beyond. The central idea is this: stability is about trapping trajectories that start close. It is not about a system "wanting" to return anywhere. This is a purely metric, epsilon-delta notion. It is structurally identical to the definition of continuity — a fact we exploit for intuition. Asymptotic stability adds a stronger demand. Nearby trajectories must actually converge to the equilibrium. They must not merely stay near it forever. This article gives the formal definitions. It proves the basic logical relationships between them. Asymptotic stability implies stability, but not conversely. Stability does not imply attraction. It works through fully explicit planar examples. It shows exactly what breaks when each hypothesis is dropped. It sets up the machinery — without yet using it. This prepares for Lyapunov Functions and Stability. It also prepares for the qualitative analysis of Conservative Systems and the Pendulum.

Key takeaways
  • An equilibrium is Lyapunov stable if trajectories starting close enough stay arbitrarily close forever.
  • It is asymptotically stable if it is stable and nearby trajectories converge to it as t→∞t\to\infty.
  • It is unstable if it fails to be stable — some trajectories starting arbitrarily close eventually escape a fixed neighborhood.
  • Stability alone does not imply convergence (centers are stable but not asymptotically stable).
  • Convergence alone does not imply stability — a subtle trap without the uniform "stay close" clause.
  • Key idea: stability is an ε\varepsilon–δ\delta trapping condition on trajectories, exactly parallel to continuity, and it is logically weaker than — and a prerequisite for — asymptotic convergence.

Where this lives today

Stability definitions of exactly this ε\varepsilon–δ\delta form appear across control theory, where "the closed-loop system is stable" is a certificate engineers must prove before deploying autopilots, power grid controllers, or chemical process regulators — instability of an equilibrium there literally means a plane, reactor, or grid runs away. In numerical analysis, the stability of a fixed point of an iterative scheme (Newton's method, gradient descent) is analyzed with the same definitions applied to discrete-time maps. In population biology and epidemiology, the stability of a disease-free equilibrium versus an endemic equilibrium (governed by the basic reproduction number R0R_0) is a direct application of Theorem (3)'s eigenvalue criterion to the linearized Jacobian of a compartmental model. In economics, the stability of market equilibria under tâtonnement dynamics is analyzed with identical machinery. Cryptographic and coding-theoretic dynamical constructions (e.g. stability of fixed points in certain pseudorandom generators modeled as dynamical systems) borrow the same vocabulary, and the entire apparatus is the direct prerequisite for the constructive, non-eigenvalue-based method of Lyapunov Functions and Stability, which is in turn essential for analyzing systems like the undamped and damped pendulum in Conservative Systems and the Pendulum where equilibria are not hyperbolic and eigenvalue tests alone are silent.

Ask the tutor: Can you show me how to determine, from scratch, whether the origin of x˙=y−x3, y˙=−x−y3\dot x = y - x^3,\ \dot y = -x - y^3 is asymptotically stable, without using eigenvalues?

DefinitionStability of an equilibrium

Let x˙=f(x)\dot{\mathbf{x}} = f(\mathbf{x}) be an autonomous system, where f:D→Rnf: D \to \mathbb{R}^n is continuous on an open set D⊆RnD \subseteq \mathbb{R}^n, and suppose ff is locally Lipschitz on DD so that solutions through each initial point exist and are unique on their maximal interval of existence. Let x∗∈D\mathbf{x}^* \in D be an equilibrium, i.e. f(x∗)=0f(\mathbf{x}^*) = \mathbf{0}. Write ϕ(t,x0)\phi(t, \mathbf{x}_0) for the solution with ϕ(0,x0)=x0\phi(0,\mathbf{x}_0) = \mathbf{x}_0.

(Lyapunov stable.) The equilibrium x∗\mathbf{x}^* is stable if for every ε>0\varepsilon > 0 there exists δ=δ(ε)>0\delta = \delta(\varepsilon) > 0 such that

∥x0−x∗∥<δ  ⟹  ∥ϕ(t,x0)−x∗∥<εfor all t≥0,\|\mathbf{x}_0 - \mathbf{x}^*\| < \delta \implies \|\phi(t,\mathbf{x}_0) - \mathbf{x}^*\| < \varepsilon \quad \text{for all } t \ge 0,

and in particular ϕ(t,x0)\phi(t,\mathbf{x}_0) is defined for all t≥0t\ge 0 whenever ∥x0−x∗∥<δ\|\mathbf{x}_0-\mathbf{x}^*\|<\delta.

(Unstable.) The equilibrium is unstable if it is not stable: there exists ε0>0\varepsilon_0 > 0 such that for every δ>0\delta > 0 there is some x0\mathbf{x}_0 with ∥x0−x∗∥<δ\|\mathbf{x}_0 - \mathbf{x}^*\| < \delta and some time t1≥0t_1 \ge 0 (possibly depending on x0\mathbf{x}_0) with ∥ϕ(t1,x0)−x∗∥≥ε0\|\phi(t_1,\mathbf{x}_0) - \mathbf{x}^*\| \ge \varepsilon_0.

(Asymptotically stable.) The equilibrium is asymptotically stable if it is stable and there exists ρ>0\rho > 0 such that

∥x0−x∗∥<ρ  ⟹  lim⁡t→∞ϕ(t,x0)=x∗.\|\mathbf{x}_0 - \mathbf{x}^*\| < \rho \implies \lim_{t\to\infty}\phi(t,\mathbf{x}_0) = \mathbf{x}^*.

If ρ\rho can be taken to be ∞\infty (any finite starting point, or all of D=RnD = \mathbb{R}^n), the equilibrium is globally asymptotically stable.

Intuition

Think of stability as a game between an adversary and you. The adversary names a tolerance ε\varepsilon — a tube of radius ε\varepsilon around the equilibrium that trajectories must never leave. Your job is to find a starting radius δ\delta small enough that any trajectory launched inside it never escapes the tube. If you can always win this game, no matter how small ε\varepsilon is, the equilibrium is stable. If the adversary can always find a trajectory that eventually breaks out of some fixed tube, no matter how close to the equilibrium it started, the equilibrium is unstable. This is exactly the logical shape of the ε\varepsilon–δ\delta definition of continuity of a function at a point, with "output" replaced by "trajectory value at time tt," except stability demands the bound hold uniformly over all future time t≥0t \ge 0, not just at a single input. That extra universal quantifier over time is what makes stability a genuinely dynamical concept rather than a static one. Asymptotic stability is a strictly stronger request: not only must the tube-confinement game be winnable, but trajectories starting near x∗\mathbf{x}^* must actually be reeled in to the point itself as t→∞t \to \infty. A center — like the origin for the linear system x˙=y, y˙=−x\dot x = y,\ \dot y = -x, whose orbits are circles — satisfies the confinement game (any circle of radius r0<δr_0 < \delta stays inside radius ε\varepsilon if δ≤ε\delta \le \varepsilon) but a trajectory on a circle of radius r0>0r_0>0 never approaches the origin. Hence stable but not asymptotically stable.

Figure
Stable, asymptotically stable, and unstable planar phase portraits
Stable (center) ε x₀ Asymptotically stable ε x₀ Unstable ε x₀

The central theorem

TheoremLogical relations among stability, asymptotic stability, and instability

Let x˙=f(x)\dot{\mathbf{x}} = f(\mathbf{x}) satisfy the hypotheses of the Definition of Stability above (continuous, locally Lipschitz ff on open D⊆RnD\subseteq\mathbb{R}^n, equilibrium x∗∈D\mathbf{x}^*\in D). Then:

(1) Asymptotic stability implies stability. (This is built into the definition, but the content is that the "attracting" clause could not, by itself, replace the "stable" clause — see the counterexample below for why both are needed as separate hypotheses.)

(2) Stability does not imply asymptotic stability. (Centers are the standard counterexample.)

(3) For a linear system x˙=Ax\dot{\mathbf{x}} = A\mathbf{x} with AA a real n×nn\times n constant matrix, the origin is:   (a) asymptotically stable if and only if every eigenvalue of AA has strictly negative real part;   (b) stable (but not necessarily asymptotically stable) if and only if every eigenvalue has real part ≤0\le 0, and every eigenvalue with real part exactly 00 is semisimple (its algebraic and geometric multiplicities coincide);   (c) unstable if and only if some eigenvalue has strictly positive real part, or some eigenvalue with zero real part fails to be semisimple.

Proof

We prove (1), (2), and then (3), which is the substantive classification result. The technique for (3) is direct computation with the matrix exponential together with a normal-form (Jordan) decomposition, plus explicit ε\varepsilon–δ\delta verification.

(1) Immediate from the definition: "asymptotically stable" is defined as "stable, and moreover attracting." So it logically implies stability by conjunction elimination. There is no independent content to prove here; the point of stating it is that the definition of asymptotic stability does not say merely "attracting," precisely because attraction alone can fail to give stability (see the counterexample block).

(2) Take n=2n=2, A=(01−10)A = \begin{pmatrix} 0 & 1 \\ -1 & 0\end{pmatrix}, so x˙=y, y˙=−x\dot x = y,\ \dot y = -x. The general solution with x0=(x0,y0)\mathbf{x}_0=(x_0,y_0) is

x(t)=x0cos⁡t+y0sin⁡t,y(t)=−x0sin⁡t+y0cos⁡t,x(t) = x_0\cos t + y_0 \sin t, \qquad y(t) = -x_0 \sin t + y_0 \cos t,

which one verifies by direct substitution satisfies x˙=y\dot x = y and y˙=−x\dot y = -x and matches the initial condition at t=0t=0. Then

x(t)2+y(t)2=x02cos⁡2t+2x0y0cos⁡tsin⁡t+y02sin⁡2t+x02sin⁡2t−2x0y0sin⁡tcos⁡t+y02cos⁡2t=x02+y02,x(t)^2 + y(t)^2 = x_0^2\cos^2 t + 2x_0y_0\cos t\sin t + y_0^2\sin^2 t + x_0^2\sin^2t - 2x_0y_0\sin t\cos t + y_0^2\cos^2 t = x_0^2+y_0^2,

so ∥ϕ(t,x0)∥=∥x0∥\|\phi(t,\mathbf{x}_0)\| = \|\mathbf{x}_0\| for all tt: every trajectory is a circle centered at the origin, traversed at constant speed. Given ε>0\varepsilon>0, choose δ=ε\delta = \varepsilon; then ∥x0∥<δ\|\mathbf{x}_0\|<\delta forces ∥ϕ(t,x0)∥=∥x0∥<ε\|\phi(t,\mathbf{x}_0)\| = \|\mathbf{x}_0\| < \varepsilon for all t≥0t\ge 0. This proves stability. But if x0≠0\mathbf{x}_0 \ne \mathbf{0}, then ∥ϕ(t,x0)∥=∥x0∥↛0\|\phi(t,\mathbf{x}_0)\| = \|\mathbf{x}_0\| \not\to 0, so ϕ(t,x0)\phi(t,\mathbf{x}_0) does not converge to the origin. Hence the origin is stable but not asymptotically stable, proving (2).

(3) Write ϕ(t,x0)=eAtx0\phi(t,\mathbf{x}_0) = e^{At}\mathbf{x}_0, the unique solution of the linear IVP (existence and uniqueness of this exponential-form solution is a standard fact about constant-coefficient linear systems, proved independently of stability theory, so invoking it here is not circular). Operator norm bound via Jordan form. Let A=PJP−1A = PJP^{-1} be a Jordan decomposition over C\mathbb{C} (or the real Jordan form over R\mathbb{R}; the argument is identical up to bounded real changes of basis). Then eAt=PeJtP−1e^{At} = Pe^{Jt}P^{-1}, and since P,P−1P, P^{-1} are fixed invertible matrices there are constants 0<c1≤c20<c_1\le c_2 with

c1∥eJty∥≤∥eAtx∥≤c2∥eJty∥,y=P−1x,c_1 \|e^{Jt}\mathbf{y}\| \le \|e^{At}\mathbf{x}\| \le c_2\|e^{Jt}\mathbf{y}\|, \qquad \mathbf{y} = P^{-1}\mathbf{x},

so it suffices to analyze eJte^{Jt} block by block. For a single Jordan block Jk(λ)J_k(\lambda) of size kk with eigenvalue λ=α+iβ\lambda = \alpha+i\beta,

eJk(λ)t=eαt eiβt N(t),N(t)=∑j=0k−1tjj! (Jk(0)nilpotent part)j,e^{J_k(\lambda)t} = e^{\alpha t}\, e^{i\beta t}\, N(t), \qquad N(t) = \sum_{j=0}^{k-1} \frac{t^j}{j!}\,\left(J_k(0)^{\text{nilpotent part}}\right)^j,

where N(t)N(t) has entries that are polynomials in tt of degree at most k−1k-1 (bounded contribution from the nilpotent part) and ∣eiβt∣=1|e^{i\beta t}|=1 (rotation, norm-preserving). Thus every entry of eJk(λ)te^{J_k(\lambda)t} has magnitude bounded by C(1+∣t∣)k−1eαtC(1+|t|)^{k-1}e^{\alpha t} for a constant CC depending only on kk. Proof of (a). Suppose every eigenvalue has Re⁡λi=αi<0\operatorname{Re}\lambda_i = \alpha_i < 0. Let α=max⁡iαi<0\alpha = \max_i \alpha_i <0. By the block estimate above, each block satisfies ∥eJk(λ)t∥≤C(1+t)k−1eαt\|e^{J_k(\lambda)t}\| \le C(1+t)^{k-1}e^{\alpha t}, and since (1+t)k−1e(α−α′)t→0(1+t)^{k-1}e^{(\alpha-\alpha')t}\to 0 for any α<α′<0\alpha < \alpha' < 0, polynomial growth is dominated by exponential decay: for any α′∈(α,0)\alpha' \in(\alpha, 0) there is M≥1M\ge 1 with

∥eJt∥≤Meα′tfor all t≥0,\|e^{Jt}\| \le M e^{\alpha' t}\quad\text{for all } t\ge 0,

hence ∥eAt∥≤c2Meα′t/c1=:M′eα′t\|e^{At}\| \le c_2 M e^{\alpha' t}/c_1 =: M' e^{\alpha' t} for all t≥0t \ge 0 (absorbing constants). Given ε>0\varepsilon>0, set δ=ε/M′\delta = \varepsilon/M'; then ∥x0∥<δ\|\mathbf{x}_0\|<\delta gives ∥eAtx0∥≤M′∥x0∥eα′t<M′δ=ε\|e^{At}\mathbf{x}_0\| \le M'\|\mathbf{x}_0\| e^{\alpha' t} < M'\delta = \varepsilon for all t≥0t \ge 0 since eα′t≤1e^{\alpha' t}\le 1. This proves stability, and since M′∥x0∥eα′t→0M'\|\mathbf{x}_0\|e^{\alpha' t}\to 0 as t→∞t\to\infty (as α′<0\alpha'<0), every trajectory converges to 0\mathbf 0; taking ρ=∞\rho=\infty shows this is in fact global asymptotic stability. Conversely, suppose some eigenvalue λ\lambda has Re⁡λ≥0\operatorname{Re}\lambda \ge 0. If Re⁡λ>0\operatorname{Re}\lambda>0 we show instability directly below in the proof of (c); if Re⁡λ=0\operatorname{Re}\lambda = 0, pick the corresponding eigenvector direction v\mathbf v (real or, for a complex pair, the real invariant plane it spans); the solution started on this eigenspace has ∥eAtv∥\|e^{At}\mathbf v\| either constant (semisimple case, purely oscillatory, as in the center example) or growing polynomially (non-semisimple case, e.g. a Jordan block with λ=0\lambda=0 gives solutions with a term linear in tt). In neither case does eAtx0→0e^{At}\mathbf x_0 \to 0 for x0\mathbf x_0 in that eigenspace, so asymptotic stability fails. This establishes the "only if" direction of (a). Proof of (b) and (c). Suppose every eigenvalue has Re⁡λi≤0\operatorname{Re}\lambda_i \le 0, and every eigenvalue with Re⁡λi=0\operatorname{Re}\lambda_i=0 is semisimple, meaning its Jordan blocks all have size k=1k=1. For blocks with Re⁡λ<0\operatorname{Re}\lambda<0 the estimate above still gives ∥eJk(λ)t∥→0\|e^{J_k(\lambda)t}\|\to 0. For a 1×11\times 1 block with λ=iβ\lambda = i\beta purely imaginary, eJ1(iβ)t=eiβte^{J_1(i\beta)t} = e^{i\beta t} has modulus exactly 11 for all tt — no polynomial growth factor appears since k=1k=1 means there is no nilpotent part. Hence every block of eJte^{Jt} is bounded uniformly in t≥0t\ge 0 by a constant depending only on AA (not growing polynomially), so there is M≥1M\ge 1 with ∥eAt∥≤M\|e^{At}\| \le M for all t≥0t\ge0. Given ε>0\varepsilon>0 take δ=ε/M\delta=\varepsilon/M: then ∥x0∥<δ⇒∥eAtx0∥≤M∥x0∥<ε\|\mathbf x_0\|<\delta \Rightarrow \|e^{At}\mathbf x_0\|\le M\|\mathbf x_0\|<\varepsilon for all t≥0t\ge0, proving stability. It need not be asymptotic: the center example in part (2) is exactly the case of a single semisimple purely-imaginary pair, and there ∥eAtx0∥=∥x0∥\|e^{At}\mathbf x_0\|=\|\mathbf x_0\| does not tend to 00. Finally, for (c): if some eigenvalue λ\lambda has Re⁡λ=α>0\operatorname{Re}\lambda = \alpha>0, let v≠0\mathbf v\ne 0 be a corresponding (generalized) eigenvector, and consider x0=sv\mathbf x_0 = s\mathbf v for small s>0s>0. Along the associated invariant subspace, ∥eAt(sv)∥≥s c eαt→∞\|e^{At}(s\mathbf v)\| \ge s\, c\, e^{\alpha t}\to\infty for a constant c>0c>0 (the k=1k=1 eigenvector case gives this directly from eAtv=eλtve^{At}\mathbf v = e^{\lambda t}\mathbf v, and ∣eλt∣=eαt|e^{\lambda t}| = e^{\alpha t}; for a genuine Jordan chain the growth is even faster). No matter how small ss is chosen (i.e. no matter how small δ\delta the adversary is handed), this norm eventually exceeds any fixed ε0>0\varepsilon_0>0 at some finite time t1t_1. This exhibits, for every δ>0\delta>0, a point x0\mathbf x_0 with ∥x0∥<δ\|\mathbf x_0\|<\delta whose trajectory leaves the ball of radius ε0\varepsilon_0 — verifying the instability definition. The non-semisimple purely-imaginary case was already handled in the converse direction of (a): a Jordan block of size k≥2k\ge2 with λ=iβ\lambda = i\beta contributes a factor growing like tk−1t^{k-1} (unbounded), giving trajectories with ∥eAtx0∥→∞\|e^{At}\mathbf x_0\|\to\infty along that generalized eigenspace, again witnessing instability. This exhausts all cases and completes the proof of (3), hence of the theorem. ■\blacksquare

∎

Worked examples

Example 1 (asymmetric linear system — direct verification, no shortcuts). Consider

x˙=−3x+y,y˙=−2y.\dot x = -3x + y, \qquad \dot y = -2y.

The matrix is A=(−310−2)A = \begin{pmatrix} -3 & 1 \\ 0 & -2\end{pmatrix}, upper triangular, with eigenvalues read off the diagonal: λ1=−3\lambda_1 = -3, λ2=−2\lambda_2=-2, both real and distinct, both strictly negative. We do not stop at citing Theorem (3)(a) — we verify by solving explicitly. From the second equation, y(t)=y0e−2ty(t) = y_0 e^{-2t}. Substituting into the first,

x˙+3x=y0e−2t.\dot x + 3x = y_0 e^{-2t}.

This is linear first-order; the integrating factor is e3te^{3t}, giving ddt(xe3t)=y0et\frac{d}{dt}(xe^{3t}) = y_0 e^{t}, so x(t)e3t=x0+y0(et−1)x(t)e^{3t} = x_0 + y_0(e^t - 1), hence

x(t)=(x0+y0)e−3t−y0e−2t.x(t) = (x_0 + y_0)e^{-3t} - y_0 e^{-2t}.

Check at t=0t=0: x(0)=x0+y0−y0=x0x(0) = x_0+y_0-y_0 = x_0. ✓. Now bound: ∣x(t)∣≤(∣x0∣+∣y0∣)e−3t+∣y0∣e−2t≤(∣x0∣+2∣y0∣)e−2t|x(t)| \le (|x_0|+|y_0|)e^{-3t}+|y_0|e^{-2t} \le (|x_0|+2|y_0|)e^{-2t} for t≥0t\ge0 (since e−3t≤e−2te^{-3t}\le e^{-2t}), and ∣y(t)∣=∣y0∣e−2t|y(t)| = |y_0|e^{-2t}. Using ∥x0∥1=∣x0∣+∣y0∣\|\mathbf x_0\|_1 = |x_0|+|y_0|,

∣x(t)∣+∣y(t)∣≤(∣x0∣+3∣y0∣)e−2t≤3(∣x0∣+∣y0∣)e−2t=3∥x0∥1e−2t.|x(t)|+|y(t)| \le (|x_0|+3|y_0|)e^{-2t} \le 3(|x_0|+|y_0|)e^{-2t} = 3\|\mathbf x_0\|_1 e^{-2t}.

Given ε>0\varepsilon>0, choose δ=ε/3\delta = \varepsilon/3 (using the ℓ1\ell^1 norm consistently); then ∥x0∥1<δ\|\mathbf x_0\|_1<\delta gives ∣x(t)∣+∣y(t)∣<εe−2t≤ε|x(t)|+|y(t)| < \varepsilon e^{-2t} \le \varepsilon for all t≥0t\ge0, confirming stability directly from the explicit solution. Since the bound also →0\to 0 as t→∞t\to\infty, the origin is (globally) asymptotically stable — matching Theorem (3)(a), but now proved from scratch for this specific asymmetric system rather than merely cited.

Example 2 (a nonlinear system — Lyapunov stability without any linear-algebra shortcut). Consider the planar system

x˙=−x+xy,y˙=−y−x2.\dot x = -x + xy, \qquad \dot y = -y - x^2.

The origin is an equilibrium since both right sides vanish at (0,0)(0,0). Linearizing, the Jacobian at the origin is (−100−1)\begin{pmatrix}-1 & 0\\0&-1\end{pmatrix}, so the linearization already suggests asymptotic stability — but we verify the nonlinear claim directly, since the standard linearization theorem (Hartman–Grobman-type reasoning) is a separate later result we are not entitled to assume here. Define V(x,y)=x2+y2V(x,y) = x^2+y^2 (this anticipates, but does not require citing, Lyapunov Functions and Stability; here we use it only as a direct computational device, verified from first principles). Along a trajectory,

ddtV(x(t),y(t))=2xx˙+2yy˙=2x(−x+xy)+2y(−y−x2)=−2x2+2x2y−2y2−2x2y=−2x2−2y2=−2V.\frac{d}{dt}V(x(t),y(t)) = 2x\dot x + 2y\dot y = 2x(-x+xy) + 2y(-y-x^2) = -2x^2+2x^2y - 2y^2 - 2x^2y = -2x^2-2y^2 = -2V.

The cross terms 2x2y2x^2y and −2x2y-2x^2y cancel exactly — this cancellation is a genuine feature of this specific vector field, not a generic phenomenon, and the reader should confirm it by direct substitution as above. So V˙=−2V\dot V = -2V, a linear ODE for VV itself, giving V(t)=V(0)e−2tV(t) = V(0)e^{-2t} exactly, for as long as the solution exists. Since V(t)=x(t)2+y(t)2=V(0)e−2t≤V(0)V(t) = x(t)^2+y(t)^2 = V(0)e^{-2t} \le V(0) for all t≥0t \ge 0 in the domain of existence, solutions starting in the disk V(0)<R2V(0) < R^2 remain in the closed disk of radius RR for all forward time; by the standard continuation theorem for ODEs (a solution with values remaining in a compact subset of D=R2D=\mathbb R^2 extends for all t≥0t\ge0), the solution exists for all t≥0t\ge 0. Now given ε>0\varepsilon>0, take δ=ε\delta = \varepsilon: if x02+y02=V(0)<δ2x_0^2+y_0^2 = V(0) <\delta^2 then V(t)=V(0)e−2t<δ2≤ε2V(t) = V(0)e^{-2t} < \delta^2 \le \varepsilon^2 for all t≥0t\ge0, i.e. x(t)2+y(t)2<ε\sqrt{x(t)^2+y(t)^2}<\varepsilon. This proves stability directly. Moreover V(t)=V(0)e−2t→0V(t)=V(0)e^{-2t}\to 0, forcing x(t)→0,y(t)→0x(t)\to0, y(t)\to0. Hence the origin is asymptotically stable — indeed globally, since the bound holds for every finite V(0)V(0). This example shows the definitional machinery applied to a genuinely nonlinear system where no eigenvalue computation alone would suffice to conclude anything about the full nonlinear flow.

CounterexampleWhy "attracting" alone is not "stable": dropping the stability clause from asymptotic stability

The definition of asymptotic stability requires both stability and attraction — one might wonder whether attraction (every nearby trajectory converges to x∗\mathbf x^*) implies stability on its own, making the conjunction redundant. It does not. Here is an explicit planar system, due in spirit to Vinograd, where every single trajectory near the origin converges to the origin, yet the origin is not Lyapunov stable, because trajectories are forced to travel far away first. Consider the origin as equilibrium of a vector field constructed (in polar-type coordinates, made globally smooth by an explicit formula) so that trajectories starting arbitrarily close to the origin are swept out along a path that first travels a bounded distance R0>0R_0>0 away before curving back and converging to the origin. Concretely, one such system (Vinograd's example) is

x˙=x2(y−x)+y5(x2+y2)(1+(x2+y2)2),y˙=y2(y−2x)(x2+y2)(1+(x2+y2)2),\dot x = \frac{x^2(y-x)+y^5}{(x^2+y^2)\left(1+(x^2+y^2)^2\right)}, \qquad \dot y = \frac{y^2(y-2x)}{(x^2+y^2)\left(1+(x^2+y^2)^2\right)},

for (x,y)≠(0,0)(x,y)\ne(0,0), extended by x˙=y˙=0\dot x=\dot y=0 at the origin. One can verify (this requires care, since the field is a nontrivial rational function) that:

  • every trajectory other than the origin itself converges to the origin as t→∞t\to\infty: attraction holds for every initial condition, not just those near the origin;
  • nonetheless, for a fixed small ε0>0\varepsilon_0>0 (specifically some ε0\varepsilon_0 less than the "excursion radius" R0R_0 built into the construction), every neighborhood of the origin contains initial points whose trajectories are driven outward past distance ε0\varepsilon_0 before eventually turning back.

The mechanism is that trajectories starting at points (x0,y0)(x_0,y_0) closer and closer to the origin still get funneled onto a family of loop-shaped curves passing near a fixed point (R0,0)(R_0,0) before spiraling in; as (x0,y0)→(0,0)(x_0,y_0)\to(0,0), the maximum excursion distance from the origin along the trajectory does not shrink to 00 — it stays bounded below by some fixed positive constant. This is precisely a failure of the uniformity in δ\delta built into the stability definition: attraction is a statement about t→∞t\to\infty for each fixed starting point, but stability additionally demands a δ\delta that works uniformly, at every finite time, not just eventually. The lesson: the definition of asymptotic stability must explicitly conjoin stability with attraction as two separate clauses, because in dimension n≥2n\ge 2 these are logically independent conditions — dropping the stability clause and defining "asymptotic stability" as mere attraction would wrongly certify Vinograd's origin as well-behaved, when in fact a tiny perturbation can produce an arbitrarily large (though finite and transient) excursion.

Scope, edge cases & common traps

SituationWhat to checkCommon error
Purely imaginary eigenvalues, all simpleStable, not asymptotically stable (center-type)Concluding "eigenvalues have Re⁡≤0\operatorname{Re}\le0" ⇒\Rightarrow asymptotic stability
Purely imaginary eigenvalue in a Jordan block of size ≥2\ge2Unstable (polynomial-in-tt growth)Assuming Re⁡λ=0\operatorname{Re}\lambda=0 always gives at least stability
Nonlinear system, linearization has a zero-real-part eigenvalueLinearization is inconclusive; need higher-order/Lyapunov methodsReading off stability from the linear part alone
"Attracting" trajectories found by simulationMust separately check the uniform δ\delta-clauseDeclaring asymptotic stability from attraction alone (Vinograd-type trap)
δ\delta depending on tt instead of only on ε\varepsilonDefinition requires one δ\delta working for all t≥0t\ge0 simultaneouslySilently allowing δ\delta to shrink as tt grows, which is not a valid δ(ε)\delta(\varepsilon)

The mathematician behind it

Origins · 1 minAleksandr Mikhailovich Lyapunov and how this was discovered
AM
Aleksandr Mikhailovich Lyapunov
1857–1918 · Imperial Russia, mathematical physics and mechanics

Lyapunov worked in St. Petersburg under Chebyshev's influence, applying rigorous analysis to problems of rotating fluid masses and mechanical equilibria. His 1892 doctoral dissertation, The General Problem of the Stability of Motion, introduced the precise definitions treated in this article and the energy-function method that bears his name, unifying what had been case-by-case mechanical arguments into a single analytic framework applicable to any system of ordinary differential equations.

History
- 1857: Lyapunov born in Yaroslavl, Russian Empire.
- 1892: doctoral dissertation formalizes stability, asymptotic stability, and the Lyapunov function method.
- 1930s: Western dynamical systems community (via translations and Malkin, La Salle, Lefschetz) absorbs and extends the theory.
- 1960s: extensions to non-autonomous and infinite-dimensional systems (Krasovskii, Hale) built directly on these definitions.
- 1918: Lyapunov dies in Odessa.
Why does the definition require δ\delta to depend only on ε\varepsilon, not on tt?

If δ\delta were allowed to shrink as tt increases, the condition would become vacuous or would collapse to something much weaker than intended: one could always satisfy "closeness at each individual time" by choosing an ever-smaller starting neighborhood for each new time horizon, without ever controlling the whole future trajectory at once. The entire content of stability is that a single, fixed tolerance on the initial condition suffices to control the trajectory uniformly over all future time. This uniformity is what makes stability a meaningful long-term guarantee (e.g., for an engineer certifying a controller) rather than a trivial pointwise statement, and it is exactly the feature that fails in the Vinograd counterexample.

Is asymptotic stability the same as "lim⁡t→∞ϕ(t,x0)=x∗\lim_{t\to\infty}\phi(t,\mathbf x_0)=\mathbf x^* for all nearby x0\mathbf x_0"?

No — that convergence condition alone is called attractivity or attraction, and it is strictly weaker than asymptotic stability. Asymptotic stability requires attraction plus Lyapunov stability, as the counterexample of Vinograd's system shows: every trajectory there converges to the origin, yet the origin is not stable, because trajectories starting arbitrarily close can be driven arbitrarily far away (though only temporarily) before returning. The two clauses are logically independent in dimension n≥2n \ge 2, so both must be stated and verified separately.

Can a nonlinear system be asymptotically stable even when its linearization has a zero eigenvalue?

Yes, and this is precisely why the linear eigenvalue test of Theorem (3) does not settle nonlinear stability questions when the linearization is non-hyperbolic (has an eigenvalue with zero real part). In such cases the higher-order (nonlinear) terms determine the outcome, and one typically needs the Lyapunov function method — constructing an explicit energy-like function VV and checking the sign of V˙\dot V along trajectories, exactly as done informally in Example 2 — rather than relying on the linear part alone. This is the motivating gap that Lyapunov Functions and Stability is built to fill.

Does instability mean every trajectory near the equilibrium runs away?

No. Instability is a much weaker (easier to satisfy) statement: it only requires that some fixed ε0>0\varepsilon_0>0 exists such that for every δ>0\delta>0, at least one trajectory starting within δ\delta eventually leaves the ε0\varepsilon_0-ball. A saddle point is the classic example: it has an entire stable manifold (a whole line, in the planar case) of initial conditions whose trajectories converge to the equilibrium, yet the equilibrium is still classified as unstable because other nearby trajectories (along the unstable manifold) diverge, and the existence of even one escaping family is enough to violate the "for all tt, stay within ε\varepsilon" clause for those points.

How does this relate to the classification of planar fixed points I already know?

The classification into nodes, saddles, spirals, and centers (from Classification of Planar Fixed Points) is a classification by geometric shape of orbits near a linear equilibrium, based on eigenvalues and eigenvectors. Stability, by contrast, is a classification by long-term metric behavior: whether trajectories stay close, converge, or escape. The two classifications align in the hyperbolic case — stable nodes and stable spirals are asymptotically stable, saddles and unstable nodes/spirals are unstable — but centers occupy a special role: geometrically a distinct orbit shape (closed loops), and stability-wise, stable but not asymptotically stable, illustrating that the two classification schemes answer genuinely different questions.

Why must ff be locally Lipschitz in the definition, rather than merely continuous?

Local Lipschitz continuity (or an equivalent condition) guarantees existence and uniqueness of the solution ϕ(t,x0)\phi(t,\mathbf x_0) through each initial point, which the stability definition silently presupposes by writing "the solution" as a well-defined single trajectory. Without uniqueness, multiple solutions could emanate from the same initial point — some staying close to the equilibrium, others escaping — and it would become ambiguous whether the equilibrium is "stable," since the answer would depend on which branch of solutions one tracks. Local Lipschitz continuity is the standard sufficient condition (via the Picard–Lindelöf theorem) that rules out this pathology.

Is it possible for an equilibrium to be neither stable nor unstable?

No — under the standing hypotheses (existence and uniqueness of solutions), stability and instability are logical negations of each other by construction: instability is defined as the failure of the stability condition. Every equilibrium of a well-posed system (unique solutions, defined at least locally in time) is classified as exactly one of "stable" or "unstable"; there is no third option at this level. What can happen is that stability holds without asymptotic stability (the center case), which is a genuinely intermediate degree of good behavior, but it still falls entirely inside the "stable" category as defined.

Before this

Where this leads