VII.H Solving Differential Equations by Separation of Variables
© 2014/2023 M. Flashman
Motivation: Our model for the growth of a population $P$ in Chapter VI.A
led us to the differential equation $ \frac {dP}{dt} = P$ with the initial condition
that $P(0) = 1$. Although this model might seem reasonable over a
short period of time, with many situations the relation of the
population size to some aspects of the environment that cannot change
such as physical space and food resources can restrain the growth
substantially. Our initial model may not be very accurate for any long
period of time. In fact, if the population is greater than the capacity
of the environment, one might expect to see the population size diminish
rather than grow.
What kind of growth rate would be appropriate for a situation with a
population growing in a way that depends both on the size of the
population and a population capacity fixed by factors other than the
population size.
One differential equation that might reflect this situation
is $\frac {dP}{dt} = P ( A - P )$ with the initial condition $P(0)
= 1$. [This
equation is called the logistic equation.] Here the constant $A$
represents the population capacity, so that $A-P$ is the excess (or
deficit ) of capacity over the current population at time $t$. This
factor
is a control on the rate of population growth which otherwise would
depend solely on $P$, the population size at time $t$.
Figure VII.H.1.S Evaluate Sage to see the tangent field
for $\frac {dP}{dt} = P ( 2 - P )$ with the integral curve through
$(0,1)$.
|
Figure VII.H.1 The tangent field
for $\frac {dy}{dx} = y ( 2 - y )$ with the integral curve through
$(0,1)$.
|
By routine analysis of $\frac {dP}{dt}$, the differential equation
seems to
make sense. When $A - P > 0$ the population will increase, while when
$A
- P < 0$ the population will decrease. We notice also that when $P =
0$ or $P = A$ the population will not grow without some external
disturbance. We would like to answer the usual questions for this model
as well, namely what do we expect for $P(1)$ and what if any is
the
long run population limit?
For a visual approach to this
differential equation we can consider the tangent field in Figures
VII.H.1S created with SAGE and VII.H.1 created with GeoGebra where $A =
2$ and the integral curve drawn through the initial
point $(0,1)$.
Now what we would like best is to have a formula for the
population as an explicit function of time.
In this section we
will develop a simple technique for solving this and many other
differential equations called separation of variables.
We'll begin with
some simpler examples to understand this technique before we return to
solve our motivating example.
Example VII.H.1. Suppose $P'(t) = 3 t^2$. There is no difficulty in
finding the general solution to this differential equation.
$P(t) = \int P'(t) dt = \int 3 t^ 2 dt = t^ 3
+ C$.
This seems simple enough. To find $P$ we merely integrate $P'(t)$ with
respect to $t$.
Here's another way to look at this: $\frac {dP}{dt} = 3 t^
2$, so
$P = \int \frac {dP}{dt}dt = \int 3 t ^2
dt = t^ 3 + C$.
Recalling the differential notation here we can replace $\frac {dP}{dt}
dt$ with the differential $dP$ and the last equation becomes
$P
= \int dP = \int 3 t^2 dt = t^ 3 + C$.
Now at least formally in terms of manipulating the symbols one might follow what we've done as follows:
● we have taken the original equation $\frac {dP}{dt} = 3 t^ 2$
● multiplied both sides by dt to obtain $\frac {dP}{dt}dt = 3 t^2 dt$ ,
● cancelled the factor $dt$ (even though it really isn't
a factor) resulting in the differential equation $dP = 3 t^ 2 dt$.
● Now integrating the differentials gives the final equation:
$P = \int 1 dP = \int 3 t ^2 dt = t^ 3 + C.$
Because in the differential equation the variable $P$ appears on one side
while the variable $t$ appears on the other, we say that in this equation
we have separated the variables.
Comment: Don't be fooled by the notation here. Though the differential
$dt$ seems to be cancelled what has happened is that the differential $dP$
is hiding a rather trivial substitution for the sake of simplifying the
integral of the expression $P'(t)$.
Here's an example where we truly need to separate the variables.
Figure VII.H.2.S Evaluate Sage to see the tangent field
for $\frac {dP}{dt} = \frac{3 t^
2} {2 P}$ with the integral curve through
$(0,1)$.
|
Example VII.H.2. Suppose $P'(t) = \frac {3 t ^2 } {2 P(t)}$ where
$P(t) \ne 0$ for any $t$ and $P(0) = 1$. To solve this differnetial equation, we begin by expressing the
equation in the Liebniz differential notation as $\frac {dP}{dt} = \frac{3 t^
2} {2 P}$ . Now multiply both sides of the equation by $P
dt$ to obtain the differential equation $2 P dP = 3 t^2 dt$. See!
We've separated the variables. Now we integrate both sides of the
equation to find
$ P^ 2 = \int 2 P dP = \int 3 t ^2 dt = t ^3 + C$.
This gives us an implicit equation for $P$ which we can solve for two
possible continuous solutions, $P = \sqrt{t^ 3 + C}$ and $P = - \sqrt{t^ 3 +
C}$ . Finally the initial condition that $P(0) = 1$ means that $C =1$
and therefore $P = \sqrt{t^ 3 + 1}$.
For a visual approach to this
differential equation we can consider the tangent field in Figures
VII.H.2S created with SAGE and VII.H.2 created with GeoGebra with the integral curve drawn through the initial
point $(0,1)$.
Figure VII.H.2 The tangent field
for $\frac {dy}{dx} = \frac{3x^2}{2y}$ with the integral curve through
$(0,1)$.
|
Let's summarize the process used in the last example.
Step 0. To use the technique called separation of variables the
differential equation must have the form $P'(t) = \frac{g(t)}{r(P(t))}$. This is sometimes not so easy to recognize.
Step 1. Write the equation in Leibniz form as $ \frac {dP}{dt} = \frac{g(t)}{r(P)}$.
Step 2. Multiply both sides $r(P) dt$ to obtain the differential
equation $r(P) dP = g(t) dt$. See! You've separated the
variables.
Step 3. Integrate both sides to obtain an equation giving an implicit
relation between $P$ and $t$. [This may not be so easy. Remember that
integration is sometimes just not possible in elementary terms.]
Step 4. If possible solve the equation of step 4 explicitly for $P$ as a function of $t$.
Example VII.H.3. Let's apply this technique to solving the
differential equation for a population model where the rate of growth is
proportional to the population $P$ at any time $t$ and where the population
initially is $1000$ while after one hour the population is $1200$.
Solution: First we express the rate statement as a differential
equation, namely that $ \frac {dP}{dt} = k P$ where $k$ is a constant
of proportionality. Now we separate the variables to obtain the
differential equation $\frac1P dP = k dt$. Now integrate:
$\ln│P│= \int \frac1P dP = \int k dt = k t +
C$.
This gives us an implicit relation between $P$ and $t$. We'll use this to
find the constants before going any further to solve explicitly for $P$.
When $t = 0$ we were told $P = 1000$. Thus $\ln(1000) = k 0 + C$.
When $t = 1$ we have $P = 1200$ , so $\ln(1200) = k⋅1 + ln(1000)$.
Thus $k = \ln(1200) - \ln(1000) = \ln(\frac {1200}{1000} = \ln (\frac 65 )$.
So we have $\ln │P│ = \ln(\frac 65 t + \ln(1000)$.
Writing this relation
using the natural exponential function we have $P = e^{ \ln(\frac 65)t +
ln(1000)} = e^{ \ln(\frac 65)t} ⋅ e^{ln(1000)}= 1000 (\frac65)^ t$.
So we have found that then $P = 1000 (\frac65)^ t$.
It is not hard to see that this
solution matches the initial conditions and the differential equation , $P
' = 1000 (\frac65)^t ⋅ ln(\frac65) = ln(\frac65) P$.
Comment: The solution to the differential equation $\frac {dP}{dt} =
kP$ with $P(0) = P_0$ may be found similarly to be $P = P_0e^{kt}$ . This
has application to numerous other contexts besides population which
are explored further in the
exercises.
Now we return to our initial motivation for this section.
Example VII.H.4. Suppose $\frac {dP}{dt} = (2 - P)P = 2P - P ^2$ with $P(0) = 1$.
Find $P(1)$ and discuss the value of $P$ when $t > > 0$.
Solution: We begin by separating the variable to obtain the differential equation
$\frac1{ (2 - P) P } dP = 1 dt$ .
Integrating the left side of this equation uses the techniques for
integrating rational functions.
It is not difficult to check
that $\int \frac 1{ (2 - P)
P }dP = \frac12 ln│P│ - \frac12 ln│2-P│+ C = \frac12 ln│\frac P{2-P}│+C$.
Since the integral of $1 dt$ is $t + C$, we have the implicit relation $\frac12 ln│\frac P{2-P}│ = t + C$.
Using the initial condition that when $t =
0, P = 1$, we find that $\frac12 ln(1) = 0 + C$ so $C=0$ and $ln│\frac P{2-P}│ = 2t$.
Now we can transform the equation using the natural
exponential function to find that $\frac P{2-P} = e^{2t}$.
Now a little algebra will allow us to solve for $P$.
Multiply both sides
by $2-P$ so $P = (2-P) e^{2t}$ and thus $2 e ^{2t} = P + P e^{2t} = P
(1 + e ^{2t})$.
Therefore $P = \frac{2e^{2t}} {1+e^{2t}}= \frac2 {1 + e^{ -2t}}$ .(Did you
follow the last step in algebra?)
Now we can answer our questions.
$P(1) = \frac 2{1 + e ^{-2}}$ and when $t > > 0 , P(t) ≈ 2$, i.e. $\lim_{ t→∞} P(t) = 2 $
.
Exercises VII.H:
In exercises 1 - 5, give (at least) an equation that implicitly solves
the differential equation. Assume $P(0) = 3$ and either find
exactly or estimate $P(1)$. Draw a tangent field for the differential
equation and illustrate on that figure the graph of the particular
solution where $P(0) = 3$.
- $P'(t) = \frac1{3 P(t)^2} P(t) \ne 0$.
- $P'(t) = \frac t{3 P(t)^2} P(t) \ne 0$.
- $P'(t) = P(t)^2$ .
- $P'(t) = t P(t)^2$ .
- $P'(t) = 2(P(t) - 1)$
- Using the technique of separation of variables solve the logistic
equation : $\frac {dP}{dt} = 3 P ( 10 - P )$ with a)$P(0) = 5$
and b) $P(0) = 15$.
In each case find exactly and estimate P(1) and P(100).
- Using the technique of separation of variables solve the general
logistic equation : $\frac {dP}{dt} = k P ( A - P )$ with $P(0) = P_0$.
- Suppose a population satisfies the logistic equation and $P(0) = 100$
and $P(1) = 200$. Discuss briefly the solution to the logistic equation
with relation to this information and what, if any, additional
information you would need to estimate $P(100)$.
- Suppose a population satisfies the differential equation
$\frac {dP}{dt} = k P$ with $P(0)=100 and P(1)=200$. Solve this
differential
equation and estimate $P(100)$.
- (Project) Go to the library and look up the populations of the
United States, China , India, and Egypt for the years 1920, 1940,
and 1980. Based on this data for each country make estimates for the
population in 1990 based on an exponential growth model and the logistic
model. Which estimate was closer to the correct population
according to your source. Draw a graph illustrating the information you
have assembled and the models you are using for your estimates. Based on
this information give an estimate for the populations of these
countries in the year 2000.
- Use separation of variables to solve the differential equation $\frac {dA}{dt}=kA$ assuming that $A(0)=A_0$.
Radioactive Decay: Radioactive materials decay with time in the sense
that in a given amount of a certain matter the mass of radioactive
matter it contains will decrease. One can presume that this decay
is due in large part to the interaction of the radioactive matter. The
radioactive atoms emit subatomic particles of matter as they decay
triggering further radioactive reactions in other radioactive atoms.
This suggests that the rate of decay for a radioactive mass will be
proportional to the amount of radioactive material present. This can be
verified experimentally for most radioactive matter. If we let $A(t)$
denote the amount of radioactive matter present at time $t$, then the
rate of decay would be $\frac{dA}{dt} = A'(t)$ and the equation that expresses the
relation of the rate of decay to the material present is $\frac{dA}{dt} = k A$
where $k$ is a constant of proportionality. Notice that for this to make
sense in terms of the value of $A$ decreasing, the constant $k$ should be a
negative number.
Half Life: The half life of a radioactive material is the time it takes
for one half of the material to decay. Thus if $A(t)$ is the amount of
material present at time $t$, and $t_h$ is the half life,
then
$A(t_h) = \frac12 A(0)$.
- Suppose $A(t)$ is the amount of a radioactive isotope measured in
grams present at time t hours and initially there are 100 grams
present but that after one hour only 75 grams remain.
a) Estimate the amount of material that will remain after 50, 100 and 500 years.
b) Estimate the number of years before there will be 100, 10, and 1 kilogram of the material remaining.
- Suppose the half life of a radioactive material is 20 years and we have 1000 kilograms of the material presently.
a) Find a formula expressing $A(t)$ explicitly as a function of time.
b) Estimate how many grams of the isotope will remain after 2, 5, and 10 hours.
c) Estimate when there will be exactly 50 , 10 and 1 gram of the isotope remaining.
- Suppose the half life of a radioactive isotope is $t_h$. Give a
simple exponential formula (without the number $e$ ) for $A(t)$ in terms of
$t, t_h$ and $A_0 = A(0)$.
Radioactive Isotope Dating: Many organic and inorganic materials contain
small but measurable amounts of radioactive substances which at some
point in time will not increase. [For example, traces of the radioactive
isotope carbon-14 are found in human bones at a fairly constant
proportion of all the carbon in the body. These amount of these traces
in the body change only by ingestion of carbon-14 present in food.
After death the only changes in the carbon-14 present in a body result
from radioactive decay. ] Measuring the presence of these isotopes in an
object where the decay is presumed and knowing the levels of these
isotopes when the decaying started allows a scientist to estimate
how long it has been since the moment when the decay began. The next
problems though fictitious in facts illustrate how this dating is done.
- Suppose that the amount of the radioactive isotope USH-92 present in
a skeleton was found to be 20 milligrams whereas at death a normal
skeleton would have 85 milligrams . The half-life of USH-92 has been
found to be 5000 years.
a) Give an estimate for the age of the skeleton, i.e. if $A(t)$ is the
amount of USH-92 present in a skeleton $t$ years after death, find a value
$t*$ so that using $A(0) = 85$ as an initial condition, $A(t*) = 20$.
b) Our devices for measuring USH-92 are not very accurate. It has been
found that any measurement can be as much as 2 milligrams away from the
correct measurement. What is the possible error in the dating of
the skeleton's age resulting from an inaccurate measurement of the
USH-90 ?
c) If we purchased a more accurate USH measuring device, we could
improve the accuracy so the error would be no more than $.5$ milligrams.
How would this improve the accuracy of our age estimates?
d) Current technology for detecting the presence of USH-92 in bones does
not allow us to measure accurately less than 5 milligrams of
USH-92. What is the oldest possible bone that can be detected
accurately using USH-92 dating.
e) The half-life of USH-90 has recently been reviewed and one group of
scientists now believe its half life is 5500 years. Redo the
previous parts of this problem using this new half life for USH-92.