L20-exponential.mws
Calculus I
Lesson 20: Exponential Growth and Decay
Suppose we model the growth or decline of a population with the following differential equation. That is, the rate of growth is proportional to the amount present. Let's solve this equation for y.
.
Then,
=
=> ln(y) =
.
Hence,
=
and setting
we have
.
Notice that
.
Conclusion: The solution to the differential equation is
Example 1
A bacteria culture starts with 500 bacteria and after 3 hours there are 8000
bacteria (assume exponential grwoth model).
a) Find an expression for the number of bacteria after t hours.
b) Find the number of bacteria after 4 hours.
c) When will the population reach 30,000?
d) Plot the expression from part (a)
e) from the Plot in (d) what happens to the population as times increases?
We have that
and need to find k. We have:
and hence
. Thus,
.
>
k:= ln(80/5)/3;
>
evalf(%);
>
f1:= t -> 500 * exp(.9241962406 *t);
>
plot(f1(t), t = 0..5, color = red);
>
f1(4);
(b) Population after 4 hours is:
>
fsolve(500 * exp(.9241962406 *t)= 30000,t);
(c) From Maple we have a numerical estimate that at t about
,
the population is 30,000. More precisely:
30,000 =
=>
=>
.
Thus at time
we have the population being 30,000.
Remember that k is about
. Using Maple we have:
>
evalf(ln(60)/ .9241962406);
(e) As times increases the population size goes to infinity.
Example 2
Polonium-214 has a half life of 1.4 x
seconds.
a) If a sample has a mass of 50 mg, find a formula for the mass that remains
after t seconds.
b) How long would it take for the mass to decay to 40 mg?
c) Will Maple plot the formula from (a); if not, why?
We have:
25 =
=>
=> ln(.5) = .00014k.
Thus,
.
>
k:= ln(.5)/(.00014);
a) Formula is: y = 50*e^(
*t)
b) We want: 40 = 50*e^(
*t) Thus,
t = ln(4/5) / (
).
>
evalf(ln(4/5) / (-4951.051290));
Hence it takes about
seconds for 40 mg to remain.
>
f2:= t -> 50 * e^(-4951.051290 * t);
>
plot(f2(t), t = 0..5);
>
evalf(e^(-4951.051290));
The graph of f2 very quickly approaches the x-axis. Hence we do not
get a good plot from Maple (for the x values chosen).
Example 3
Scientists can determine the age of ancient objects by a method called
radiocarbon dating. The bombardment of the upper atmosphere by cosmic
rays converts bitrogen to a radioactive isotope of carbon, C-14, with a half
life of about 5730 years. Vegetation absorbs carbon dioxide through the
atmosphere and animal life assimilates C-14 through food chains.
When a plant or animal dies it stops replacing its carbon and the amount
of C-14 begins to decrease through radioactive decay. Therefore, the level
of radioactivity must also exponentially decay. A parchment fragment was
discovered that had about 74% as much C-14 radioactivity as does plant material
on the earth today. Estimate the age of the parchment.
Let y(t) denote the level of radioactivity. Then
.
Using the half-life information we have:
=>
.
Thus we have:
=> t =
.
Thus the age of the parchment is:
.
Using Maple for a numerical estiamte of the age we obtain:
>
>
evalf( ln(.74) / ( ln(.5)/5730 ) );
Hence the parchment is about 2489 years old.
Example 4
Newtons Law of Cooling states that the rate of cooling of an object is
proportional to the temperature difference between the object and its surroundings.
A roast turkey is taken from the oven when tis temperature has reached 185 F
and is placed on a table in a room where the temperature is 75 F.
a) If the temperature of the turkey is 150 F after half an hour,
what is the temperature after 45 minutes?
b) When will the turkey have cooled to 100F?
c) Find a formula for the temperature of the turkey at time t.
d) Plot the formula from (c).
Let y(t) be the temperature of the turkey after t minutes. Then
Newtons Law of Cooling states that
.
Then
and hence
.
Thus,
We have that
y(30) = 150 =
=> l
.
Thus, the temperature of the turkey after t minutes is:
, where
.
>
evalf(ln(75/110)/30);
(a) We want y(45).
>
y:= t -> 75 + 110 * exp(-.01276640841 * t);
>
y(45);
Thus the turkey is about
F after 45 minutes.
(b) We want to know when y(t) = 100.
>
fsolve( 75 + 110 * exp(-.01276640841 * t) = 100, t);
Numerically we have that it takes 1 hour and 56 minutes for the turkey to cool to 100F.
Without Maple we can solve by:
=>
=>
.
Thus it takes
minutes for the turkey to cool to 100 F, recall
.
>
evalf( ln(25/110) / ( ln(75/110) / 30));
Note that our answers agree.
>
plot(y(t),t = 0..120, color = red);