where \(F\) is a function of \(n + 2\) variables. A solution to this equation on an interval \(I = (a, b)\) is a function \(u = u(t)\) such that the first \(n\) derivatives of \(u\) are defined on \(I\text{,}\) and

\begin{equation*}
F(t, u, u', u'', \ldots, u^{(n)}) = 0.
\end{equation*}

We will concentrate on first-order differential equations in this chapter. That is, we will consider equations of the form

For some first-order differential equations such as

\begin{equation*}
\frac{dx}{dt} = k x
\end{equation*}

the solution is given by an explicit formula. In this particular case, the solution is \(x(t) = Ce^{kt}\text{.}\) In general, we cannot generally find such a formula for an arbitrary first-order differential equation. We can, however, solve a differential equation \(y' = f(x, y)\) if we can write the equation in the form

Such equations are called separable. We can solve separable equations by integrating the first term with respect to \(x\) and the second term with respect to \(y\text{.}\)

Example1.2.1

Suppose that we wish to solve the initial value problem

where \(C\) is an arbitrary constant. Using the initial condition, \(y(0) = 1\) to find \(C\text{,}\) we see that

\begin{equation*}
0 = \ln|1| = 0 + C.
\end{equation*}

Thus, the solution to our initial value problem can be given implicitly by \(\ln y = x^2/2\text{.}\) In this example, we can actually write down an explicit solution that is defined everywhere,

\begin{equation*}
y = e^{x^2/2}.
\end{equation*}

The Sage commands for solving our initial value problem are below.

If we ask Sage to solve a differential equation that is impossible to solve analytically, the computation will return an error. Sage will not solve the initial value problem

Notice that the solution does not make sense for all values of \(t\text{.}\) In fact, the solution is only defined on the interval \(-1 \lt t \lt 1\text{,}\) if we require that our solution be continuous. Let us see what Sage has to say.

Sage does return a solution even if it looks a bit different than the one that we arrived at above. Notice that we have an imaginary term in our solution, where \(i^2 = -1\text{.}\) We will examine the role of complex numbers and how useful they are in the study of ordinary differential equations in a later chapter, but for the moment complex numbers will just muddy the situation.

Example1.2.3

The initial value problem in Example 1.1.2 is a good example of a separable differential equation,

Since \(C\) is an arbitrary constant, we know that \(e^C\) is an arbitrary postive constant, which we will also call \(C\text{.}\) So we can rewrite this last equation as

Separable equations arise in a wide range of application problems. One does not have to watch too many crime dramas to realize that the time of death of a murder victim is an important question in many criminal investigations. How does a forensic scientist or a medical examiner determine the time of death? Human beings have a temperature of \(98.6^\circ\)F. If the surrounding temperature is cooler, then the body will cool down after death. Eventually, the temperature of the body will match the temperature of the environment. We should not expect the body to cool at a constant rate either. Think of how a hot cup of coffee or tea cools. The liquid will cool quite quickly during the first few minutes but will remain relatively warm for quite a long period.

The answer to our forensic question can be found by using Newton's law of cooling, which tells us that the rate of change of the temperature of an object is proportional to the difference between the temperature of the object and the temperature of the surrounding medium. Newton's law of cooling can be easily stated as a differential equation,

where \(T\) is the temperature of the object, \(T_m\) is the temperature of the surrounding medium, and \(k\) is the proportionality constant.

Suppose that the temperature of the surrounding environment is \(70^\circ\)F, and we know from experience that a body under these conditions cools off approximately \(2^\circ\)F during the first hour after death. In order to determine a formula for the time of death, we must solve the initial value problem

we see that this equation is separable. Integrating both sides of the last equation, we obtain \(\ln | T - 70| = kt + C\text{.}\) Since we are assuming that \(T \gt 70\text{,}\) we can write \(T - 70\) instead of \(|T - 70|\text{.}\) Thus, we have

There is a large class of problems in modeling known as mixing problems. These problems refer to situations where two or more substances are mixed together in a container or containers. For example, we might wish to model how chemicals are mixed together in a refinery, how pollutants are mixed together in a pond or a lake, how ingredients are mixed together when brewing beer, or even how various greenhouse gases mix together across different layers of the atmosphere.

Suppose that we have a large tank containing 1000 gallons of pure water and that water containing 0.5 pounds of salt per gallon flows into the tank at a rate of 10 gallons per minute. If the tank is also draining at a rate of 10 gallons per minute, the water level in the tank will remain constant. We will assume that the water in the tank is constantly stirred so that the mixture of salt and water is uniform in the tank.

We can model the amount of salt in the tank using differential equations. If \(x(t)\) is the amount of salt in the tank at time \(t\text{,}\) then the rate at which the salt is changing in the tank is the difference between the rate at which salt is flowing into the tank and the rate at which it is leaving the tank, or

Of course, the salt flows into the tank at the rate of \(10 \cdot 0.5 = 5\) pounds of salt per minute. However, the rate at which the salt leaves the tank depends on \(x(t)\text{,}\) the amount salt in the tank at time \(t\text{.}\) At time \(t\text{,}\) there is \(x(t)/1000\) pounds of salt in one gallon. Therefore, salt flows out of the tank at a rate of \(10x(t)/1000 = x(t)/100\) pounds per minute. Equation (1.2.2) now becomes

Differential equations have many applications in economics and finance. For example, Dr. J., a college professor, wisely started saving for his retirement as soon as he entered the workforce, and he now has $500,000 in a retirement account earning an interest of 5% compounded continuously. The initial value problem,

provides a nice model of Dr. J.'s investment, where \(P(t)\) is the amount in thousands of dollars in the fund at time \(t\text{.}\) The solution to our initial value problem is

If Dr. J. plans to retire in 10 years, he can expect a nest egg of \(P(10) \approx 824.360635350064\) or about $824,360.

Of course, Dr. J. still plans to make contributions to his retirement fund during his next ten years of employment. His annual contribution will be $5,000, which his employer will generously match. If we assume that these contributions will spread out evenly over the course of the year, we can incorporate this information into our original initial value problem,

Dr. J.'s nest egg is now \(P(10) \approx 954.104889490090\) or about $954,105.

Once Dr. J. retires, he will need to begin withdrawing money from his account. He estimates that he will need to withdraw $60,000 a year for living expenses if he wishes to travel and enjoy his golden years. Of course, whatever remains in his account at any given time will still collect interest. We describe J.'s retirement situation with the initial value problem,

\begin{align*}
\frac{dP}{dt} \amp = \begin{cases} 0.05 P + 10, \amp t \leq 10 \\ 0.05 P - 60, \amp t \gt 10 \end{cases}\\
P(0) \amp = 500.
\end{align*}

Hence, the rate of withdrawal exceed the rate at which Dr. J.'s account is earning interest. Eventually, Dr. J.'s retirement fund will be disappear. This may pose a problem, if Dr. J. plans to retire early and live a long life.

Again, the differential equation \(dP/dt = 0.05 P - 60\) is separable, and we have

We now give a theoretical basis for solving first-order separable differential equations. A differential equation \(y' = F(x, y)\) is called separable if it can be written in the form

We now will prove that such an equation can solved by integrating the first term with respect to \(x\) and the second term with respect to \(y\text{.}\) If

\begin{equation*}
\arctan y + \ln|y| = x + C.
\end{equation*}

However, we have no method of solving this last equation explicitly for \(y\text{.}\)

Example1.2.6

Another difficulty arises if we consider the equation

\begin{equation*}
y' = t e^{-y^2}.
\end{equation*}

This equation is separable since we can rewrite it in the form

\begin{equation*}
e^{y^2} \, dy = t \, dt.
\end{equation*}

Although the Fundamental Theorem of Calculus guarantees that every continuous function has an antiderivative, we cannot find an antiderivative for the function \(e^{y^2}\) in terms of elementary functions. Thus, we are forced to write our solution as

Even if we have a separable differential equation, we are not guaranteed a unique solution. Consider the initial value problem \(y' = y^{1/3}\) with \(y(0) = 0\) and \(t \geq 0\text{.}\) Separating the variables,

\begin{equation*}
y^{-1/3} \, dy = dt.
\end{equation*}

Thus,

\begin{equation*}
\frac{3}{2} y^{2/3} = t + C
\end{equation*}

or

\begin{equation*}
y = \left( \frac{2}{3} ( t + C)\right)^{3/2}.
\end{equation*}

If \(C = 0\text{,}\) the initial condition is satisfied and

\begin{equation*}
y = \left( \frac{2}{3} t \right)^{3/2}
\end{equation*}

is a solution for \(t \geq 0\text{.}\) However, we can find at least two additional solutions for \(t \geq 0\text{:}\)

\begin{align*}
y & = - \left( \frac{2}{3} t \right)^{3/2},\\
y & \equiv 0.
\end{align*}

In Section 1.5 we will learn sufficient conditions for a first-order initial value problem to have a unique solution.

Example1.2.8

Suppose that \(y' = y^2\) with \(y(0) = 1\text{.}\) Separating the variables,

\begin{equation*}
\frac{1}{y^2} \, dy = dt,
\end{equation*}

we see that

\begin{equation*}
y = - \frac{1}{t + C}
\end{equation*}

or

\begin{equation*}
y = \frac{1}{1-t}.
\end{equation*}

Therefore, a continuous solution also exists on \((-\infty, 1)\) if \(y(0) = 1\text{.}\) In the case that \(y(0) = -1\text{,}\) the solution is

\begin{equation*}
y = - \frac{1}{t + 1},
\end{equation*}

and a continuous solution exists on \((-1, \infty)\text{.}\)

where \(F\) is a function of \(n + 2\) variables. A solution to the equation on an interval \(I = (a, b)\) is a function \(u = u(t)\) such that the first \(n\) derivatives of \(u\) are defined on \(I\text{,}\) and

\begin{equation*}
F(t, u, u', u'', \ldots, u^{(n)}) = 0.
\end{equation*}

A first-order differential equation is an equation that can be written in the form

Use the fact that this new equation is separable to solve for \(y\text{.}\)

Show that any homogeneous equation \(y' = f(x, y)\) can be transformed into a separable differential equation by making the substitution \(y = vx\text{.}\)

A function \(f\) is said to homogeneous of degree \(n\) if \(f(t x, t y)= t^n f(x, y)\) for \(n = 1, 2, \ldots\text{.}\) Show that differential equation

where \(P\) and \(Q\) are both homogeneous of degree \(n\text{,}\) can be transformed into a separable differential equation using the substitution \(y = vx\text{.}\)

Mr. Ratchett, an elderly American, was found murdered in his train compartment on the Orient Express at 7 A.M. When his body was discovered, the famous detective Hercule Poirot noted that Ratchett had a body temperature of 28 degrees. The body had cooled to a temperture of 27 degrees one hour later. If the normal temperature of a human being is 37 degrees and the air temperature in the train is 22 degrees, estimate the time of Ratchett's death using Newton's Law of Cooling.

Subsection1.2.8Sage—Quick Start Guide to Solving Ordinary Differential Equations

Sage has powerful algorithms for finding exact and numerical solutions of differential equations. In addition, we can plot solutions and direction fields. Although some differential equations have an exact solution and can be solved using analytic techniques with calculus, many differential equations can only be solved using numerical techniques. This should not be too surprising if we consider how we solve polynomials. It is quite easy to find the roots of any equation of the form \(ax^2 + bx + c = 0\) by either factoring or using the quadratic equation, but solving an equation such as

is a much more difficult problem. Unlike the situation for quadratic equations, there does not exist a general formula for solving seventh degree equations. We can even encounter difficulties when using a numerical method such as the Newton-Raphson algorithm.

In general, Sage needs three things to solve a differential equation:

An abstract function

A differential equation

A Sage command to solve the equation.

Suppose we wish to solve the equation

\begin{equation*}
\frac{dy}{dx} = x + y.
\end{equation*}

We can use the following sequence of Sage commands.

The first command defines the abstract function. The second describes the actual differential equation. Finally, we use the Sage command desolve to find the actual solution. Try replacing the h command with show(h) or show(expand(h)).

We can also specify an initial condition for our differential equation, say \(y(0) = 1\text{.}\)