Unit 9, lesson 116

Direction Fields, Equilibria, and Uniqueness

How to read the behavior of solutions from a slope field or a phase line without solving the equation, how to classify equilibria as stable or unstable, and when an initial-value problem has exactly one solution.

  • Theory part 231
  • Practice part 232
  • Equations and Solutions

Most differential equations cannot be solved by a formula. The equation still tells you the slope of the solution at every point, and that is often enough to see where solutions go, which values they settle on and which ones they run away from. This lesson is about reading a differential equation as a picture.

Key terms

Slope fields

Definition Let ff be defined on a region RR of the plane. The slope field (or direction field) of y′=f(x,y)y' = f(x, y) assigns to each point (x,y)(x, y) of RR a short line segment through that point with slope f(x,y)f(x, y). In practice one draws the segments at the points of a grid.

In plain words At each point the equation tells you how steep a solution through that point must be. Draw that steepness as a tiny dash.

Example For y′=x−yy' = x - y, a few slopes are below.

Point (x,y)(x, y)Slope x−yx - y
(0,0)(0, 0)00
(1,0)(1, 0)11
(0,1)(0, 1)−1-1
(2,1)(2, 1)11
(1,2)(1, 2)−1-1

Along every line y=x−cy = x - c the slope is the same number cc, so the dashes on such a line are all parallel.

Solution curves

Definition A solution curve (integral curve) of y′=f(x,y)y' = f(x, y) is the graph of a solution y=φ(x)y = \varphi(x) on an interval. At each of its points (x,φ(x))(x, \varphi(x)) the curve has slope φ′(x)=f(x,φ(x))\varphi'(x) = f(x, \varphi(x)), so it is tangent to the segment of the slope field there.

In plain words A solution curve follows the dashes, like a leaf floating along the current.

Example For y′=x−yy' = x - y, the line y=x−1y = x - 1 has slope 11, and on it x−y=x−(x−1)=1x - y = x - (x - 1) = 1. So the line is a solution curve. The general solution is y=x−1+Ce−xy = x - 1 + Ce^{-x}, and every other solution curve approaches this line as x→∞x \to \infty.

Autonomous equations

Definition A first order equation is autonomous if the right side does not depend on the independent variable, that is, it has the form

y′=g(y)y' = g(y)

for a function gg of yy alone.

In plain words The rule for how yy changes depends only on where yy is now, not on the time. So the slope field looks the same along every horizontal line.

Example The equation y′=y2−4y' = y^2 - 4 is autonomous. The equation y′=x−yy' = x - y is not, because the slope at height y=0y = 0 is 00 at x=0x = 0 but 11 at x=1x = 1.

Equilibrium solutions

Definition A constant function y(x)=cy(x) = c is an equilibrium solution of y′=f(x,y)y' = f(x, y) if f(x,c)=0f(x, c) = 0 for every xx in the interval under study. For an autonomous equation y′=g(y)y' = g(y), the equilibrium solutions are exactly y=cy = c with g(c)=0g(c) = 0.

In plain words An equilibrium is a value where the rate of change is zero, so staying at that value forever is a solution.

Example For y′=y2−4y' = y^2 - 4 the equilibria are the roots of y2−4=0y^2 - 4 = 0, namely y=−2y = -2 and y=2y = 2.

Phase line

Definition The phase line of an autonomous equation y′=g(y)y' = g(y) with gg continuous is the yy-axis with the zeros of gg marked and an arrow on each open interval into which these zeros cut the axis (here gg is assumed to have finitely many zeros). The arrow points up where g(y)>0g(y) > 0 and down where g(y)<0g(y) < 0, because there every solution increases or decreases.

In plain words It is the slope field squeezed into a single vertical line. It shows which way yy moves from each starting value.

Example For y′=y2−4y' = y^2 - 4, test one value in each interval.

IntervalTest valueSign of y2−4y^2 - 4Arrow
y>2y > 233++up
−2<y<2-2 < y < 200−-down
y<−2y < -2−3-3++up

So solutions starting between −2-2 and 22 fall toward −2-2, those starting above 22 rise, and those below −2-2 rise toward −2-2.

Stability

Definition Let cc be an equilibrium of y′=g(y)y' = g(y) with gg continuous, and write the independent variable as tt. It is stable if for every ε>0\varepsilon > 0 there is a δ>0\delta > 0 such that every solution with ∣y(t0)−c∣<δ\lvert y(t_0) - c \rvert < \delta exists for all t≥t0t \ge t_0 and satisfies ∣y(t)−c∣<ε\lvert y(t) - c \rvert < \varepsilon for all t≥t0t \ge t_0. It is asymptotically stable if moreover y(t)→cy(t) \to c as t→∞t \to \infty for every solution starting close enough to cc. It is unstable if it is not stable. If gg is differentiable at cc, then

g′(c)<0  ⟹  c is asymptotically stable,g′(c)>0  ⟹  c is unstable,g'(c) < 0 \implies c \text{ is asymptotically stable}, \qquad g'(c) > 0 \implies c \text{ is unstable},

and g′(c)=0g'(c) = 0 decides nothing.

In plain words A stable equilibrium is a valley. Nudge the solution a little and it rolls back. An unstable one is a hilltop, where the smallest nudge sends it away.

Example For g(y)=y2−4g(y) = y^2 - 4 we have g′(y)=2yg'(y) = 2y. Since g′(−2)=−4<0g'(-2) = -4 < 0, the equilibrium −2-2 is asymptotically stable, and since g′(2)=4>0g'(2) = 4 > 0, the equilibrium 22 is unstable. This matches the phase line. For y′=y2y' = y^2 the test fails at 00 since g′(0)=0g'(0) = 0. The phase line has arrows up on both sides, so solutions below 00 approach it and solutions above leave it. This equilibrium is unstable and is often called semistable.

Existence and uniqueness

Definition Let RR be an open set in the plane and ff continuous on RR. Suppose ff is locally Lipschitz in yy, uniformly in xx, meaning every point (a,b)(a, b) of RR has a rectangle ∣x−a∣≤α\lvert x - a \rvert \le \alpha, ∣y−b∣≤β\lvert y - b \rvert \le \beta inside RR and a constant KK with

∣f(x,y1)−f(x,y2)∣≤K∣y1−y2∣\lvert f(x, y_1) - f(x, y_2) \rvert \le K \lvert y_1 - y_2 \rvert

for all (x,y1)(x, y_1) and (x,y2)(x, y_2) in that rectangle. Then for every (x0,y0)(x_0, y_0) in RR there is an h>0h > 0 such that the problem y′=f(x,y)y' = f(x, y), y(x0)=y0y(x_0) = y_0 has exactly one solution on (x0−h,x0+h)(x_0 - h, x_0 + h). Moreover two solutions on any common interval containing x0x_0 agree there. The Lipschitz condition holds automatically when ∂f∂y\frac{\partial f}{\partial y} exists and is continuous on RR.

In plain words If the equation is continuous and does not change too violently in the yy direction, then through each starting point there is one solution curve and no other. In particular two solution curves cannot cross or touch.

Example Uniqueness can fail without the Lipschitz condition. For y′=3y2/3y' = 3y^{2/3}, y(0)=0y(0) = 0, both y=0y = 0 and y=x3y = x^3 are solutions on R\mathbb{R}, since

(x3)′=3x2=3(x3)2/3.\left(x^3\right)' = 3x^2 = 3\left(x^3\right)^{2/3}.

Here the difference quotient

∣f(x,y)−f(x,0)∣∣y−0∣=3∣y∣−1/3,y≠0,\frac{\lvert f(x, y) - f(x, 0) \rvert}{\lvert y - 0 \rvert} = 3\lvert y \rvert^{-1/3}, \qquad y \neq 0,

is unbounded near y=0y = 0, so no Lipschitz constant works there. By contrast, for y′=y2−4y' = y^2 - 4 the partial derivative 2y2y is continuous everywhere, so solutions never cross the equilibria. A solution starting at y(0)=0y(0) = 0 therefore stays strictly between −2-2 and 22 for all time.

Check yourself

  1. Find the equilibria of y′=y(3−y)y' = y(3 - y) and classify each one.
Show the answer

The function g(y)=3y−y2g(y) = 3y - y^2 has zeros 00 and 33, and g′(y)=3−2yg'(y) = 3 - 2y. Since g′(0)=3>0g'(0) = 3 > 0, the equilibrium 00 is unstable. Since g′(3)=−3<0g'(3) = -3 < 0, the equilibrium 33 is asymptotically stable.

  1. For y′=y2−4y' = y^2 - 4 with y(0)=1y(0) = 1, what does y(t)y(t) approach as t→∞t \to \infty?
Show the answer

It approaches −2-2. The starting value lies between −2-2 and 22, where y2−4<0y^2 - 4 < 0, so the solution decreases, cannot cross −2-2 by uniqueness, and approaches the stable equilibrium −2-2.

  1. Does the existence and uniqueness theorem apply to y′=yy' = \sqrt{y} at the starting point (0,0)(0, 0)? At (0,4)(0, 4)?
Show the answer

Not at (0,0)(0, 0), because ff is not even defined for y<0y < 0 and ∂f∂y=12y\frac{\partial f}{\partial y} = \frac{1}{2\sqrt{y}} is unbounded near y=0y = 0. In fact y=0y = 0 and y=x24y = \frac{x^2}{4} for x≥0x \ge 0 (with y=0y = 0 for x<0x < 0) both solve it. At (0,4)(0, 4) it applies, using the open set y>0y > 0 where ∂f∂y\frac{\partial f}{\partial y} is continuous.

How this lesson connects

It also uses slope fields.