Section1.3Geometric and Quantitative Analysis¶ permalink

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\text{,}\) then a solution curve will have a slope of \(2\) at the point \((1, 1)\text{.}\) We can use this information to obtain a geometric description of the solutions to the equation.

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.17). 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)\text{,}\) 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)\text{,}\) 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\text{,}\) and the voltage drop across the capacitor, \(E_C\text{.}\) According to Kirchhoff's Law, this is
\begin{equation*}
E_R + E_C = E.
\end{equation*}

Resistance, \(R\text{,}\) to the current is measured in ohms (\(\Omega\)). Ohm's Law tells us that the voltage drop across a resistor is given by
\begin{equation}
E_R = IR.\label{equation-firstlook03-ohms-law}\tag{1.5}
\end{equation}
Finally, capacitance, \(C\text{,}\) is measured in farads (F). Coulomb's Law tells us how current flows across a capacitor,
\begin{equation}
I = C \frac{d E_C}{dt}. \label{equation-firstlook03-coulombs-law}\tag{1.6}
\end{equation}
Thus, if we combine the equations (1.5) and (1.6), our equation \(E_R + E_C = E\) becomes
\begin{equation}
RC \frac{d E_C}{dt} + E_C = E(t).\label{equation-firstlook03-rc-circuit}\tag{1.7}
\end{equation}

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.7) for different voltage sources \(E(t)\text{.}\) If we view the differential equation (1.7) 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.

Any differential equation
\begin{equation*}
\frac{dy}{dt} = f(t, y)
\end{equation*}
can be viewed as a formula for the slope of a function \(y = y(t)\text{.}\) Geometrically, the equation tells us that, at any point \((t_0, y_0)\text{,}\) the slope of a solution curve is given by \(f(t_0, y_0)\text{.}\) Suppose that our differential equation is defined on the rectangle \(R = [a, b] \times [c, d]\text{,}\) and let \(y(t)\) be a solution curve for \(y' = f(t, y)\) that passes through the point \((t_0, y_0)\text{.}\) Then the differential equation tells us the slope of this solution curve at \((t_0, y_0)\text{.}\) 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)\text{.}\) 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\text{.}\) The direction field for this equation is given in Figure 1.18 along with several solution curves.

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.

Now let us return to our RC circuit and consider different functions \(E(t)\) for the differential equation
\begin{equation*}
\frac{d E_C}{dt} = \frac{E(t) - E_C}{RC},
\end{equation*}
with \(R = 1\) and \(C = 1\text{.}\) First suppose that there is no voltage source in the circuit. If we let \(E(t) = 0\) for all \(t \ge 0\text{,}\) we will get the direction field of given in Figure 1.19. The direction field agrees with our analytic solution
\begin{equation*}
E_C(t) = v_0 e^{-t},
\end{equation*}
where \(v_0 = E_C(0)\text{.}\)

If we assume that we have a nonzero constant source of voltage, \(E(t) = K\text{,}\) in our circuit such as a battery, then we obtain the separable differential equation
\begin{equation*}
\frac{d E_C}{dt} = K - E_C.
\end{equation*}
The direction field for this differential equation for \(K = 10\) is given in Figure 1.20.

If we attach a battery to our circuit at time \(t = 0\) and then disconnect the battery at \(t = 4\text{,}\) then we obtain a different solution. For example, if
\begin{equation*}
E(t) =
\begin{cases}
10 & 0 \leq t \leq 4, \\
0 & t \gt 4,
\end{cases}
\end{equation*}
we will obtain a different direction field (Figure 1.21).

If our voltage source emits a series of pulses, say
\begin{equation*}
E(t) =
\begin{cases}
10 & 0 \leq t \lt 4, \\
0 & 4 \leq t \lt 8, \\
10 & 8 \leq t \lt 12, \\
& \vdots
\end{cases}
\end{equation*}
then the direction field for our differential equation is given in Figure 1.22.

Finally, if we use a generator for a voltage source, the voltage source might be given by a function such as \(E(t) = \sin(\pi t /2)\text{.}\) The direction field for this circuit is given in Figure 1.23.

Subsection1.3.4Plotting Direction Fields with Sage¶ permalink

It is very tedious to plot direction fields by hand. Fortunately, we have technology to help. 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. Simply evaluate the Sage cell. You can then change each field in the cell output to solve a particular initial value problem.

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

An autonomous differential equation is one of the form
\begin{equation*}
\frac{dx}{dt} = f(x).
\end{equation*}
In other words, a differential equation is autonomous if it does not depend on $t$. Since an autonomous differential equation is time independent, its direction field is particularly easy to draw. The slope only depends on \(x\) and is the same for all values of \(t\text{.}\)

Example1.24A 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
\begin{equation*}
\frac{dP}{dt} = P\left(1 - \frac{P}{200} \right),
\end{equation*}
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
\begin{equation*}
\frac{dP}{dt} = P\left(1 - \frac{P}{200} \right) - 32.
\end{equation*}
The direction field for this equation is given in Figure 1.25.

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
\begin{equation*}
\frac{dP}{dt} = P\left(1 - \frac{P}{200} \right) - 32 = -\frac{1}{200}(P - 40)(P - 160)
\end{equation*}
is an autonomous differential equation, the direction field does not depend on \(t\text{.}\) 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.

Notice that \(dP/dt = 0\) when \(P= 40\) or \(P = 160\text{.}\) Thus, the two constant solutions \(P(t) = 40\) and \(P(t) = 160\) are the same for all values of the independent variable \(t\text{.}\) We say that such a solution is an equilibrium solution. Equilibrium solutions graph as horizontal lines. 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 (Figure 1.26). For values of \(P\) between 40 and 160, we know that \(dP/dt \gt 0\text{.}\) 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 \(dP/dt \lt 0\) when \(P \lt 40\) or \(P \gt 160\text{.}\) In this case any solution curve will be decreasing, and we will indicate this by a downward pointing arrow.

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.

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\text{,}\) we have
\begin{equation*}
\lim_{t \to \infty} y(t) = y_0.
\end{equation*}
We say that an equilibrium point is a source if all solutions that start sufficiently close to \(y_0\) tend toward \(y_0\) as \(t \to - \infty\text{.}\) An equilibrium solution that is neither a sink or a source is called a node (Figure 1.26).

Consider the differential equation
\begin{equation}
\frac{dy}{dt} = y^4 - 4y^2 = y^2(y + 2)(y - 2).\label{equation-firstlook03-sinks-sources-nodes}\tag{1.8}
\end{equation}
The graph of \(f(y) = y^4 - 4y^2\) is given in Figure 1.27. If \(y = -2\text{,}\) we have a sink. If \(y = 2\text{,}\) we have a source. Finally, if \(y = 0\text{,}\) we have a node.

It is easy to generate a phase line diagram for equation (1.8) from the graph of \(f(y) = y^2(y + 2)(y - 2)\) (Figure 1.27). 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.28).

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,
\begin{equation*}
f(x) \approx f(x_0) + f'(x_0)(x - x_0) + \frac{ f''(x_0)}{2!} (x - x_0)^2 + \cdots + \frac{ f^{(n)}(x_0)}{n!} (x - x_0)^n
\end{equation*}
near \(x_0\text{.}\) For example, if
\begin{equation*}
\frac{dx}{dt} = f(x) = \cos(x^2 + \pi).
\end{equation*}
with \(x(0) = x_0\) then we may approximate this initial value problem near \(x_0\) with
\begin{align*}
\frac{dx}{dt} & = f(x_0) + f'(x_0)(x - x_0) = \cos(x_0^2 + \pi) + 2 x _0 \sin(x_0^2 + \pi)(x - x_0)\\
x(0) & = x_0.
\end{align*}
Of course, this strategy might not work very well if \(f(x_0) = \cos(x_0^2 + \pi) = 0\) or \(f'(x_0) = 2 x _0 \sin(x_0^2 + \pi) = 0\text{.}\)

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)\text{.}\) 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\text{.}\) 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.

Find all of the equilibrium solutions for each of the following differential equations. Use Sage to graph the direction field of each equation and then superimpose a plot of the equilibrium solution(s) on the direction field. Classify each equilibrium solution as a sink, a source, or a node.

\(y' = 2y(1 - y)\)

\(\dfrac{dP}{dt} = 0.1P(P - 100)\)

\(\dfrac{dx}{dt} = (2 - x)\cos x\)

\(\dfrac{dy}{dx} = (2 - y)\cos x\)

2

Consider the differential equation \(y' = f(y)\text{,}\) where the graph of \(f(y)\) is given below. Identify any equilibrium solutions and classify each equilibrium solution as a sink, source, or node.

3

What happens if we increase the harvest rate to 100 in Example 1.24.? What should be our strategy to maintain a viable population in the trout pond?

Subsection1.3.8Sage—Plotting Direction Fields and Solutions¶ permalink

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\text{.}\)

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.

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.

Suppose that the population of a trout pond can be accurately modeled by the logistic equation
\begin{equation*}
\frac{dp}{dt} = 0.4p \left( 1 - \frac{p}{500}\right).
\end{equation*}
At time \(t = 30\text{,}\) 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 folllowing:
\begin{equation*}
\frac{dp}{dt} =
\begin{cases}
0.4p \left( 1 - \dfrac{p}{500}\right) & \text{for } 0 \leq t \lt 30; \\
0.4p \left( 1 - \dfrac{p}{500}\right) - 0.1p & \text{for } t \gt 30.
\end{cases}
\end{equation*}

Plot the direction field for this equation using Sage.

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\text{.}\)

Give a qualitative description of how the disease affects the population.