ODE-Project Geometric and Quantitative Analysis

Section 1.3 Geometric and Quantitative Analysis

Objectives

To understand that direction fields are a useful way of analyzing a differential equation from a geometric point of view, especially since not all differential equations can be solved analytically.
To understand that an autonomous equation is a differential equation of the form $y^{'} = f (y)$ and that phase lines can be used to analyze autonomous differential equations.
To understand equilibrium solutions to a differential equation $y^{'} = f (y) .$ are those solutions given by $f (y) = 0$ for all $y .$ In particular, an equilibrium solution is either a sink, source, or node.

If we view the differential equation

y^{'} = f (t, y)

as a formula for the slope of a tangent line to a solution curve, we can approximate the graph of a solution curve. For example, if we consider the equation

y^{'} = t + y,

then a solution curve will have a slope of

2

at the point

(1, 1) .

We can use this information to obtain a geometric description of the solutions to the equation.

🔗

Subsection 1.3.1 RC Circuits

Suppose that we wish to analyze how an electric current flows through a circuit. An RC circuit is a very simple circuit might contain a voltage source, a capacitor, and a resistor (Figure 1.3.1). A battery or generator is an example of a voltage source. The glowing red heating element in a toaster or an electric stove is an example of something that provides resistance in a circuit. A capacitor stores an electrical charge and can be made by separating two metal plates with an insulating material. Capacitors are used to power the electronic flashes for cameras. Current,

I (t),

is the rate at which a charge flows through this circuit and is measured in amperes or amps (A). We assign a direction to the current. A current flowing in the opposite direction will be given negative values.

🔗

The source voltage,

E (t),

is measured in volts (V). Kirchhoff’s Second Law tells us that the impressed voltage in a closed circuit is equal to the sum of the voltage drops in the rest of the circuit. Thus, we need only compute the voltage drop across the resistor,

E_{R},

and the voltage drop across the capacitor,

E_{C} .

According to Kirchhoff’s Law, this is

🔗

E_{R} + E_{C} = E .

🔗

Resistance,

R,

to the current is measured in ohms (

Ω

). Ohm’s Law tells us that the voltage drop across a resistor is given by

🔗

\begin{matrix} (1.3.1) & E_{R} = I R . \end{matrix}

Finally, capacitance,

C,

is measured in farads (F). Coulomb’s Law tells us how current flows across a capacitor,

🔗

\begin{matrix} (1.3.2) & I = C \frac{d E_{C}}{d t} . \end{matrix}

Thus, if we combine the equations (1.3.1) and (1.3.2), our equation

E_{R} + E_{C} = E

becomes

🔗

\begin{matrix} (1.3.3) & R C \frac{d E_{C}}{d t} + E_{C} = E (t) . \end{matrix}

🔗

We will now investigate how our circuit reacts under different voltage sources. For example, we might have a zero voltage source (the capacitor could still hold a charge). We could also have a constant nonzero source of voltage such as a battery or a fluctuating source of voltage such as a generator. We might even have a series of pulses of voltage where the current is periodically turned on and off. We would like to be able to understand the solutions to the equation (1.3.3) for different voltage sources

E (t) .

If we view the differential equation (1.3.3) as an expression for computing how fast current is flowing across the capacitor, we can analyze our circuit from a geometric point of view and can actually say a great deal about circuits without solving a differential equation.

🔗

Subsection 1.3.2 Direction Fields

Any differential equation

🔗

\frac{d y}{d t} = f (t, y)

can be viewed as a formula for the slope of a function

y = y (t) .

Geometrically, the equation tells us that, at any point

(t_{0}, y_{0}),

the slope of a solution curve is given by

f (t_{0}, y_{0}) .

Suppose that our differential equation is defined on the rectangle

R = [a, b] \times [c, d],

and let

y (t)

be a solution curve for

y^{'} = f (t, y)

that passes through the point

(t_{0}, y_{0}) .

Then the differential equation tells us the slope of this solution curve at

(t_{0}, y_{0}) .

We can indicate this on the

(t, y)

-plane by drawing a short line segment at the point

(t, y)

with slope

f (t_{0}, y_{0}) .

Thus, we can obtain a direction field or slope field for the differential equation. A solution curve must be tangent to its direction field at every point.

🔗

For example, consider the differential equation

y^{'} = y^{2} / 2 - t .

The direction field for this equation is given in Figure 1.3.2 along with several solution curves.

🔗

Figure 1.3.2. The direction field for $y^{'} = y^{2} / 2 - t$
🔗

Although direction fields can be tedious to compute using pencil and paper, we can easily generate direction fields for any differential equation with the use of computer software. Most computer algebra systems, including Sage, have facilities for generating and graphing direction fields. For example, the following Sage cell will generate the direction field for

y^{'} = y^{2} / 2 - t

using the command plot_slope_field. We will learn how to add solution curves at the end of this section.

🔗

xxxxxxxxxx
 
t, y = var('t, y')
f(t, y) = y^2/2 - t
plot_slope_field(f, (t, -1, 5), (y, -5, 10), headaxislength=3, headlength=3, axes_labels=['$t$','$y(t)$'])

Activity 1.3.1. Plotting Direction Fields.

Plot the direction field and the solution curves for each of the following initial conditions.

🔗

(a)

y^{'} = y - 2 x,

y (0) = 4,

y (0) = 1 .

🔗

(b)

y^{'} = y (1 + x),

y (0) = 1,

y (0) = 0,

y (0) = - 1 .

🔗

(c)

\frac{d x}{d t} + 2 t x = x,

x (0) = 5 .

🔗

(d)

\frac{d y}{d t} = 2 y (1 - y / 4),

y (0) = 1,

y (0) = 4,

y (0) = 5 .

🔗

(e)

y^{'} = - 1 - y^{4},

y (0) = 4,

y (0) = 0,

y (0) = - 4 .

🔗

In each case comment on anything that you notice about the direction field and the solutions.

🔗

Subsection 1.3.3 RC Circuits Revisited

Now let us return to our RC circuit and consider different functions

E (t)

for the differential equation

🔗

\frac{d E_{C}}{d t} = \frac{E (t) - E_{C}}{R C},

with

R = 1

and

C = 1 .

First suppose that there is no voltage source in the circuit. If we let

E (t) = 0

for all

t \geq 0,

we will get the direction field of given in Figure 1.3.3. The direction field agrees with our analytic solution

🔗

E_{C} (t) = v_{0} e^{- t},

where

v_{0} = E_{C} (0) .

🔗

Figure 1.3.3. The direction field for no current
🔗

If we assume that we have a nonzero constant source of voltage,

E (t) = K,

in our circuit such as a battery, then we obtain the separable differential equation

🔗

\frac{d E_{C}}{d t} = K - E_{C} .

The direction field for this differential equation for

K = 10

is given in Figure 1.3.4.

🔗

Figure 1.3.4. The direction field for a constant current
🔗

If we attach a battery to our circuit at time

t = 0

and then disconnect the battery at

t = 4,

then we obtain a different solution. For example, if

🔗

E (t) = {\begin{cases} 10 & 0 \leq t \leq 4, \\ 0 & t > 4, \end{cases}

we will obtain a different direction field (Figure 1.3.5).

🔗

Figure 1.3.5. The direction field for a single pulse
🔗

If our voltage source emits a series of pulses, say

🔗

E (t) = {\begin{cases} 10 & 0 \leq t < 4, \\ 0 & 4 \leq t < 8, \\ 10 & 8 \leq t < 12, \\ ⋮ \end{cases}

then the direction field for our differential equation is given in Figure 1.3.6.

🔗

Figure 1.3.6. The direction field for a series of pulses
🔗

Finally, if we use a generator for a voltage source, the voltage source might be given by a function such as

E (t) = \sin (π t / 2) .

The direction field for this circuit is given in Figure 1.3.7.

🔗

Figure 1.3.7. The direction field for an oscillating voltage
🔗

🔗

Subsection 1.3.4 Plotting Direction Fields with Sage

The Sage interact below will help you plot direction fields for a given initial value problem. A Sage interact is a menu driven Sage applet. All of the Sage code is “under the hood.” You can then change each field in the cell output to solve a particular initial value problem.

🔗

Figure 1.3.8. A Sage applet for plotting direction fields

🔗

Of course, it is quite easy to give the interact an initial value problem that yield unitelligible results.

🔗

Subsection 1.3.5 Autonomous Differential Equations

An autonomous differential equation is one of the form

🔗

\frac{d x}{d t} = f (x) .

In other words, a differential equation is autonomous if the variable

t

does not appear on the righthand side of the equation. Since an autonomous differential equation

d x / d t = f (x)

only depends on the variable

x,

its direction field is particularly easy to graph. The slope only depends on

x

and is the same for all values of

t .

🔗

Example 1.3.9. A Logistic Model with Harvesting.

Let us consider a trout pond that has a carrying capacity of 200 fish. Suppose that the trout population can be modeled according to the logistic equation

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}),

where

t

is the time in years. If we allow the fish to be harvested at a constant rate of 32 per year, our equation becomes

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) - 32.

The direction field for this equation is given in Figure 1.3.10.

🔗

Figure 1.3.10. Logistic growth with harvesting

One of the basic questions that we can ask of our model is whether or not we have a sustainable population in our trout pond given this harvest rate. If so, under what conditions for sustainablility?

🔗

Since

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) - 32 = - \frac{1}{200} (P - 40) (P - 160)

is an autonomous differential equation, the direction field does not depend on

t .

Consequently, we need only keep track of what happens on the vertical axis. We can do this with a phase line. Instead of drawing the entire direction field, we can draw a single line containing the same information (Figure 1.3.11).

🔗

Notice that

d P / d t = 0

when

P = 40

P = 160 .

Thus, the two constant solutions

P (t) = 40

and

P (t) = 160

are the same for all values of the independent variable

t .

We say that such a solution is an equilibrium solution. Equilibrium solutions graph as horizontal lines on the direction field. We can identify equilibrium solutions by setting the derivative of the function equal to zero. On our phase line we will represent these solutions as equilibrium points. For values of

P

between 40 and 160, we know that

d P / d t > 0 .

Thus, any solution curve must be increasing. We denote this property on the phase line by drawing an upward pointing arrow. On the other hand, we know that

d P / d t < 0

when

P < 40

P > 160 .

In this case any solution curve will be decreasing, and we will indicate this by a downward pointing arrow.

🔗

Let

y^{'} = f (y)

and suppose that

y = y_{0}

is an equilibrium solution. We say this solution is a sink if for any solution

y (t)

with initial condition sufficiently close to

y_{0},

we have

🔗

lim_{t \to \infty} y (t) = y_{0} .

We say that an equilibrium point is a source if all solutions that start sufficiently close to

y_{0}

tend toward

y_{0}

t \to - \infty .

An equilibrium solution that is neither a sink or a source is called a node (Figure 1.3.11). When

P = 40,

we have a source, and when

P = 160,

we have a sink.

🔗

An equilibrium solution is stable if a small change in the initial conditions gives a solution which tends toward the equilibrium as the independent variable tends towards positive infinity. An equilibrium solution is unstable if a small change in the initial conditions gives a solution which veers away from the equilibrium as the independent variable tends towards positive infinity.

🔗

Consider the differential equation

🔗

\begin{matrix} (1.3.4) & \frac{d y}{d t} = y^{4} - 4 y^{2} = y^{2} (y + 2) (y - 2) . \end{matrix}

The graph of

f (y) = y^{4} - 4 y^{2}

is given in Figure 1.3.12. If

y = - 2,

we have a sink. If

y = 2,

we have a source. Finally, if

y = 0,

we have a node.

🔗

Figure 1.3.12. Sinks, sources, and nodes
🔗

It is easy to generate a phase line diagram for equation (1.3.4) from the graph of

f (y) = y^{2} (y + 2) (y - 2)

(Figure 1.3.12). If the graph is above the

y

-axis, then

y

is increasing. If the graph is below the

y

-axis, then

y

is decreasing. Therefore, the phase line is easy to sketch (Figure 1.3.13).

🔗

Figure 1.3.13. Phase line diagram for $y^{'} = y^{2} (y + 2) (y - 2)$
🔗

Activity 1.3.2. Autonomous Equations and Phase Lines.

For each of the differential equations below, draw the phase line and classify each equilibrium solution as a sink, a source, or a node.

🔗

(a)

y^{'} = y (y - 2) (y + 3) .

🔗

(b)

y^{'} = y^{2} (y - 2) (y + 3) .

🔗

(c)

y^{'} = \cos y .

🔗

(d)

y^{'} = \cos^{2} y .

🔗

In each case comment on anything that you notice about the phase line and the equilibrium solutions.

🔗

One of the reasons why autonomous equations are so important is Taylor’s theorem, which tells us that any function

f (x)

can be approximated near a point

x_{0}

by an

n

th degree polynomial,

🔗

f (x) \approx f (x_{0}) + f^{'} (x_{0}) (x - x_{0}) + \frac{f^{″} (x_{0})}{2!} (x - x_{0})^{2} + \dots + \frac{f^{(n)} (x_{0})}{n!} (x - x_{0})^{n}

near

x_{0} .

For example, if

🔗

\frac{d x}{d t} = f (x) = \cos (x^{2} + π)

with

x (0) = x_{0},

then we may approximate this initial value problem near

x_{0}

with

🔗

\begin{aligned} \frac{d x}{d t} & = f (x_{0}) + f^{'} (x_{0}) (x - x_{0}) = \cos (x_{0}^{2} + π) + 2 x_{0} \sin (x_{0}^{2} + π) (x - x_{0}) \\ x (0) & = x_{0} . \end{aligned}

Of course, this strategy might not work very well if

f (x_{0}) = \cos (x_{0}^{2} + π) = 0

f^{'} (x_{0}) = 2 x_{0} \sin (x_{0}^{2} + π) = 0 .

🔗

Subsection 1.3.6 Important Lessons

Direction fields and phase lines are a useful way of analyzing a differential equation from a geometric point of view, especially since not all differential equations can be solved analytically.
An autonomous equation is a differential equation of the form $y^{'} = f (y) .$ We can use a phase line to analyze autonomous differential equations.
Equilibrium solutions to a differential equation $y^{'} = f (y)$ are those solutions given by $f (y) = 0$ for all $y .$ In this case, any solution must be constant. We can classify equilibrium solutions according to whether they are stable or unstable. In particular, an equilibrium solution is either a sink, source, or node.

🔗

Reading Questions 1.3.7 Reading Questions

1.

Explain why we do not expect solution curves to intersect for a differential equation.

🔗

2.

Explain in your own words what an autonomous differential equation is.

🔗

Exercises 1.3.8 Exercises

Plotting direction fields by hand.

For each of the differential equations in Exercise Group 1.3.8.1–6, plot the direction field on the integer coordinates

(t, x)

of the rectangle

- 2 < t < 2

and

- 2 < x < 2

by drawing a short line of the appropriate slope.

🔗

1.

x^{'} = x + t

🔗

2.

x^{'} = x t

🔗

3.

x^{'} = x^{2} + t^{2}

🔗

4.

x^{'} = t + \tan (x)

🔗

5.

x^{'} = (x + t) / (x^{2} + t^{2})

🔗

6.

x^{'} = x - t + 1

🔗

7.

Use Sage to plot the direction fields of each differential equation in Exercise Group 1.3.8.1–6.

🔗

xxxxxxxxxx
 
t, x = var('t, x')
f(t,x) = x + t
plot_slope_field(f, (t, -2, 2), (x, -2, 2), headaxislength=3, headlength=3, axes_labels=['$t$','$x(t)$'])

🔗

Equilibrium solutions and phase lines.

Find the equilibrium solutions for each of the differential equations in Exercise Group 1.3.8.8–13. Draw the phase line for each equation and classify each equilibrium solution as a sink, a source, or a node.

🔗

8.

y^{'} = 2 y - 5

🔗

9.

\frac{d x}{d t} = (x - 1) (x + 2)

🔗

10.

\frac{d x}{d t} = (x^{2} - 1) (x - 2)

🔗

11.

\frac{d y}{d x} = \sin 2 y

🔗

12.

x^{'} = (x^{2} + 1) (x - 1)

🔗

13.

x^{'} = x^{2} + x + 1

🔗

Sketching solutions.

Each of the differential equations in Exercise Group 1.3.8.14–19 has several initial conditions specified. Sketch the solution curves that satisfy the initial conditions. Sketch your solutions for each equation on the same pair of axes.

🔗

14.

y^{'} = 2 y - 5,

y (0) = 4,

y (0) = 2,

y (0) = 0 .

🔗

15.

\frac{d x}{d t} = (x - 1) (x + 2),

x (0) = 2,

x (0) = 1,

x (0) = 0 .

🔗

16.

\frac{d x}{d t} = (x^{2} - 1) (x - 2),

x (0) = 1,

x (0) = 0,

x (0) = - 1 .

🔗

17.

\frac{d y}{d x} = \sin 2 y,

y (0) = 2,

y (0) = 1,

y (0) = 0 .

🔗

18.

x^{'} = (x^{2} + 1) (x - 1),

x (0) = 2,

x (0) = 1,

x (0) = 0 .

🔗

19.

x^{'} = x^{2} + x + 1,

x (0) = 1,

x (0) = 0,

x (0) = - 1 .

🔗

Phase lines from graphs of the derivative.

Consider the differential equation

y^{'} = f (y),

where the graph of

f (y)

is given in Exercise Group 1.3.8.20–23. Draw the phase line for each equation and classify each equilibrium solution as a sink, a source, or a node.

🔗

24.

What happens if we increase the harvest rate to 100 in Example 1.3.9? What should be our strategy to maintain a viable population in the trout pond and still permit fishing?

🔗

Subsection 1.3.9 Sage—Plotting Direction Fields and Solutions

Subsubsection 1.3.9.1 Plotting direction fields

If we wish to plot a direction field for a differential equation, we can use the command plot_slope_field. Let us plot the direction field for the equation

y^{'} = y^{2} / 2 - t .

🔗

xxxxxxxxxx
 
t, y = var('t, y')
f(t, y) = y^2/2 - t
v = plot_slope_field(f, (t,-1,5), (y,-5,10), headaxislength=3, headlength=3, axes_labels=['$t$','$y(t)$'])
v

There are a few extra commands to specify the size of the arrows in the plot and to label the axes. Try changing or omitting these commands and see what happens.

🔗

xxxxxxxxxx
 
t, y = var('t, y')
f(t, y) = y^2/2 - t
v = plot_slope_field(f, (t,-1,5), (y,-5,10), headaxislength=3, headlength=3, axes_labels=['$t$','$y(t)$'])
v

For more examples and options, see doc.sagemath.org/html/en/reference/plotting/sage/plot/plot_field.html

🔗

Subsubsection 1.3.9.2 Plotting solutions

Now let us find a numerical solution to the equation using the command desolve_rk4. This is a fourth-order Runge-Kutta method, and returns a numerical solution (a table of values). Here, we must supply the dependent variable and initial conditions.

🔗

xxxxxxxxxx
 
t, y = var('t, y')
f(t, y) = y^2/2 - t
p = desolve_rk4(f, y, ics=[-1,0], ivar=t, output='plot', end_points=[-1,5], thickness=2)
p

Of course, we can combine the two plots.

🔗

xxxxxxxxxx
 
t, y = var('t, y')
f(t, y) = y^2/2 - t
p = plot_slope_field(f, (t,-1,5), (y,-5,10), headaxislength=3, headlength=3, axes_labels=['$t$','$y(t)$'], fontsize=12)
p += desolve_rk4(f, y, ics=[-1,0], ivar=t, output='plot', end_points=[-1,5], thickness=2)
p.show(xmin = -1, xmax = 5, ymin = -5, ymax = 10)  #set the size of the plot window

There are many other commands and packages to solve ordinary differential equations using Sage. For more information, see www.sagemath.org/doc/reference/calculus/sage/calculus/desolvers.html. Below is an empty Sage cell, where you can practice.

🔗

xxxxxxxxxx
 
​

🔗

Exercises 1.3.9.3 Sage Exercises

1.

Suppose that the population of a trout pond can be accurately modeled by the logistic equation

🔗

\frac{d p}{d t} = 0.4 p (1 - \frac{p}{500}) .

At time

t = 30,

a disease is introduced into the population that kills 10% of the population per year. To see how the disease affects the fish population, we will change our original model to the following:

🔗

\frac{d p}{d t} = {\begin{cases} 0.4 p (1 - \frac{p}{500}) & for 0 \leq t < 30; \\ 0.4 p (1 - \frac{p}{500}) - 0.1 p & for t > 30. \end{cases}

Plot the direction field for this equation using Sage.
xxxxxxxxxx
Language:
Messages
Plot the graphs of two or three representative solutions to this equation on the direction field.
Find formulas for the solutions of this equation for initial conditions $p (0) = 30 .$
Give a qualitative description of how the disease affects the population.

🔗

Prev Top Next

Section 1.3 Geometric and Quantitative Analysis

Objectives

Subsection 1.3.1 RC Circuits

Subsection 1.3.2 Direction Fields

Activity 1.3.1. Plotting Direction Fields.

(a)

(b)

(c)

(d)

(e)

Subsection 1.3.3 RC Circuits Revisited

Subsection 1.3.4 Plotting Direction Fields with Sage

Subsection 1.3.5 Autonomous Differential Equations

Example 1.3.9. A Logistic Model with Harvesting.

Activity 1.3.2. Autonomous Equations and Phase Lines.

(a)

(b)

(c)

(d)

Subsection 1.3.6 Important Lessons

Reading Questions 1.3.7 Reading Questions

1.

2.

Exercises 1.3.8 Exercises

Plotting direction fields by hand.

1.

2.

3.

4.

5.

6.

7.

Equilibrium solutions and phase lines.

8.

9.

10.

11.

12.

13.

Sketching solutions.

14.

15.

16.

17.

18.

19.

Phase lines from graphs of the derivative.

20.

21.

22.

23.

24.

Subsection 1.3.9 Sage—Plotting Direction Fields and Solutions

Subsubsection 1.3.9.1 Plotting direction fields

Subsubsection 1.3.9.2 Plotting solutions

Exercises 1.3.9.3 Sage Exercises

1.