Section7 4notes
Section7 4notes
In the next two sections, we examine how population growth can be modeled using
differential equations. We start with the basic exponential growth and decay models.
Before showing how these models are set up, it is good to recall some basic background
ideas from algebra and calculus.
dP
= kP .
dt
1 dP
k=
P dt
dP
= kP .
dt
dP
= kdt (Separate the variables)
P
1
∫ P dP = ∫ kdt (Integrate both sides)
dP
= kP, P (0) = P0
dt
P(t ) = P0 e kt
Notes[
1. When modeling a population with an exponential growth model, if the relative
growth rate k is unknown, it should be determined. This is usually done using the
known population at two particular times.
2. Exponential growth models are good predictors for small populations in large
populations with abundant resources, usually for relatively short time periods.
3. The graph of the exponential equation P(t ) = P0 e kt has the general form
P (t )
P0
Example 1: Solve a certain organism develops with a constant relative growth of 0.2554
per member per day. Suppose the organism starts on day zero with 10 members. Find the
population size after 7 days.
Solution:
█
4
Example 2: A population of a small city had 3000 people in the year 2000 and has grown
at a rate proportional to its size. In the year 2005 the population was 3700.
a. Find an expression for the number of people in the city t years after the year 2000.
b. Estimate the population of the city in 2006. In 2010.
c. Find the rate of growth of the population in 2006.
d. Assuming the growth continues at the same rate, when will the town have 25000
people?
Solution:
5
█
6
Exponential Decay
When a population decays exponentially, it decreases at a rate that is proportional to its
size at any time t. The model for exponential decay is
dP
= −kP, P(t ) = P0
dt
1 dP
Here, k = − is called the relative decay constant. Note that k > 0 since, because
P dt
dP 1 dP
the population is decreasing, < 0 and k = − ⋅ > 0 . Using separation of
dt N P dt
N
negative negative
variables in a process similar to exponential growth, it can be shown that the solution to
the initial value problem is P(t ) = P0 e − kt . Summarizing, we have the following:
dP
= −kP, P(0) = P0
dt
P(t ) = P0 e − kt
Notes
1. Many times the rate of decay is expressed in terms of half-life, the time it takes for
half of any given quantity to decay so that only half of its original amount remains.
2. Radioactive elements typically decay exponentially.
7
Solution: (Part a) Since this is an exponential decay problem, we will use the formula
P(t ) = P0 e − kt . Since we start with 800 mg, then we know that P0 = 800 . Thus the
formula becomes
P(t ) = 800e − kt
To complete the equation that models this population, we need to find the relative decay rate k.
We can use the half life of the substance to do this. The half life of Bismuth-210 is 5 days. This
says that after t = 5, the original population of 800 mg has decay to half of its original amount,
1
or (800) = 400 mg. Mathematically, since P (t ) represents that amount of population of the
2
substance left after time t, this says that P(5) = 400 . Using the decay equation, we have
or rearranging, we have
We must solve this equation for k. We proceed with the following steps.
400
e −5 k = (Divide by sides by 800)
800
e −5k = 0.5 (Simplify)
ln e −5k = ln(0.5) (Take ln of both sides)
− 5k ln e = ln(0.5) (Use property ln b u = u ln b)
− 5k (1) = ln(0.5) (Recall ln e = 1)
ln(0.5)
k= (Divide both sides of - 5)
−5
k ≈ 0.1386 (Use calculator and round to 4 decimail places)
Substituting k = 0.1386 and P0 = 800 gives a formula for finding the remaining mass.
P (t ) = 800e −0.1386t
(continued on next page)
8
Part b.) Using the formula P (t ) = 800e −0.1386t found in part a, we see that
Mass remaining
= P(30) = 800e −0.1386(30) = 800e −4.158 ≈ 12.5 grams
after t = 30 days
Part c.) In this problem, we want the time t it takes for the mass to have reduced down to
1 mg. That is, we want t when P(t ) = 1 . We perform the following steps using
P (t ) = 800e −0.1386t to solve for t.
Thus, it takes approximately t = 48.2 days for the substance to decay to 1 mg.
> P := 800*exp(-0.1386*t);
( −0.1386 t )
P := 800 e
> plot(P, t = 0..60, color = blue, thickness = 2, view = [-
1..60, -10..1000], title = "Graph of 800e^(-0.1386t) for t
= 0..60 for Bismuth-210");
█
10
Example 4: Radiocarbon Dating. Scientists can determine the age of ancient objects (fossils,
for example) using radiocarbon dating. The bombardment of the upper atmosphere by cosmic
rays converts nitrogen to a radioactive isotope of carbon, 14C , with a half life of about 5730
years. Vegetation absorbs carbon dioxide through the atmosphere and animal life assimilates
14
C through food chains. When a plant or animal dies, it stops replacing its carbon and the
amount of 14C begins to decrease through radioactive decay. Therefore, the level of
radioactivity must also decay exponentially. Suppose a fossil found has about 35 % as much
14
C radioactivity as normal animals do on Earth today. Estimate the age of the fossil.
Solution:
11
█
12
dT
= k (T − Ts )
dt
Suppose we let y = T − Ts . Then, taking the derivative of both sides with respect to the
dy dT dT dT dy
time t gives = −0 = (remember, Ts is constant). Substituting = and
dt dt dt dt dt
y = T − Ts into the Newton Law of cooling model gives the equation
dy
=k y.
dt
This is just the basic exponential growth model. The solution of this differential equation is
y (t ) = y 0 e kt ,
T (t ) − Ts = (T0 − Ts )e kt
Solving for T (t ) gives the solution to the Newton Law of Cooling differential equation:
T (t ) = Ts + (T0 − Ts )e kt
The rate that the temperature T of an object that is cooling is given by the initial value
problem
dT
= k (T − Ts ), T (0) = T0 .
dt
The particular solution of this initial value problem describing the objects temperature
after a particular time t is given by
T (t ) = Ts + (T0 − Ts )e kt
where
Example 5: Suppose you cool a pot of soup in a 75 0 F room. Right when you take the
soup off the stove, you measure its temperature to be 220 0 F. Suppose after 20 minutes,
the soup has cooled to 170 0 F.
a. What will be the temperature of the soup in 30 minutes.
b. Suppose you can eat the soup when it is 130 0 F. How long will it take to cool to this
temperature?
dT
Solution: (Part a) The solution to initial value problem = k (T − Ts ), T (0) = T0
dt
describing Newton’s Law of Cooling is
T (t ) = Ts + (T0 − Ts )e kt
Since the temperature of the room is 75 0 , Ts = 75 . Since the soup when taken off the stove is
220 0 , T0 = 220. Thus, the formula for describing the temperature after time t is given by
T (t ) = 75 + (220 − 75)e kt
or, when simplified,
T (t ) = 75 + 145e kt
As with any exponential model, if the proportionality constant k, is not given, we must
find it. To do this, we use the fact that after t = 20 minutes, the temperature is 170 0 .
Mathematically, this says that T (20) = 170 . We substitute this fact into the above
equation and solve for k using the following steps.
T (t ) = 75 + 145e −0.0211t
Thus, we find the temperature after t = 30 substituting into this equation. Hence,
Temperature of Soup
= T (30) = 75 + 145e −0.0211(30)
after t = 30 minutes
= 75 + 145e −0.633
= 152 0 F
Part b) We want the time t it takes for the soup to cool to 130 0 , that is, the time t when
T (t ) = 130 . Taking the equation we found in part a T (t ) = 75 + 145e −0.0211t , we solve this
equation as follows: