Video Transcript
In many natural phenomena,
quantities grow or decay at a rate proportional to their size. For example, investing money in a
bank account which pays compound interest is a good example of exponential
growth. Though perhaps not quite at the
rate most savers would like. At the other end, atmospheric
pressure decreases exponentially with increasing height above sea level, an example
of exponential decay. In this video, weβll learn how to
model exponential growth and decay arising from differential equations and use these
models to make predictions about future behavior.
In general, we say that if π¦ of π‘
is the value of a quantity π¦ at time π‘, and if the rate of change of π¦ with
respect to π‘ is proportional to its size, π¦ of π‘, then dπ¦ by dπ‘ is equal to π
times π¦ for real constants π. This is sometimes called the law of
natural growth, for values of π greater than zero, and the law of natural decay,
for values of π less than zero. And we can solve this differential
equation generally. Remember, whilst dπ¦ by dπ‘ isnβt
really a fraction, under certain circumstances we do treat it a little like one, and
we rewrite our differential equation as one over π¦ dπ¦ equals π dπ‘.
We now integrate both sides of this
equation. Using the standard result that the
integral of one over π¦ is the natural log of the absolute value of π¦, we see that
we have the natural log of the absolute value of π¦ plus π΄ equals π times π‘ plus
π΅, where π΄ and π΅ are constants of integration. We combine these constants by
subtracting π΄ from both sides. And letβs call this new constant
πΆ. So, the natural logarithm of the
absolute value of π¦ is equal to π times π‘ plus πΆ. We now raise both sides of this
equation as a power of π, such that π to the power of the natural log of the
absolute value of π¦ equals π to the power of ππ‘ plus πΆ.
Well, π to the natural log of π¦
is π¦. And we donβt need the absolute
values here because weβre talking about a population. It must be greater than zero by
definition. Then, using the laws of exponents,
we see that we can rewrite π to the power of ππ‘ plus πΆ as π to the power of
ππ‘ times π to the power of πΆ. But of course, π to the power of
πΆ, given that πΆ is already a constant, is itself a constant. So, letβs call that π. And we find that the general
solution to the differential equation dπ¦ by dπ‘ equals π times π¦ is π¦ equals π
times π to the power of ππ‘.
Now, in fact, this π has a special
significance. We can identify what this
significance is by imagining the scenario when π‘ is equal to zero. Clearly, since π‘ equals zero is
the initial time, we know that π¦ equals π times π to the power of zero to be the
initial value of π¦. However, π to the power of zero is
one. So, we find that when π‘ is equal
to zero, π¦ is equal to π. This means that π is the initial
value of the function. Now, we redefine a little and we
say that, in general, the solution to the initial value problem dπ by dπ‘ equals π
times π for π of zero equals π nought is π of π‘ equals π nought times π to
the ππ‘th power. In order to see how to work with
these equations, letβs have a look at a specific example.
The population of rabbits on a farm
grows exponentially. If there are currently 245 rabbits
and the relative growth rate is 23 percent, find a function π of π‘ to describe the
number of rabbits after π‘ years.
In this question, we have our
variables defined. Weβre told that π‘ is the number of
years, and π of π‘ is the function that describes the number of rabbits. We know that since the population
grows exponentially, we can use the equation dπ by dπ‘ equals π times π to
describe it. In fact though, weβll rewrite this
using our variables such that dπ by dπ‘ equals π times π. We could have written this using
function notation such that π prime of π‘ equals π times π. But personally, I prefer to work
with Leibniz notation.
Now, in fact, we have a value for
π. We know that the relative growth
rate is 23 percent. Thatβs 0.23 as a decimal. And we can therefore say that dπ
by dπ‘ equals 0.23 times π. Now remember, once dπ by dπ‘ is
absolutely not a fraction, there are certain circumstances where we treat it a
little like one. Here, we can rewrite our
differential equation as one over π dπ equals 0.23 dπ‘. Weβre then going to integrate both
sides of this equation. The integral of the left-hand side
is the natural log of π plus some constant of integration π. And the integral of the right-hand
side is 0.23π‘ plus a different constant of integration π.
Now, note that usually the integral
of one over π would be the natural log of the absolute value of π. But we know the number of rabbits
canβt be negative, so we donβt actually need to include that here. Combining our constants of
integration by subtracting π from both sides, and we see the natural log of π is
equal to 0.23π‘ plus π. We then raise both sides to the
power of π to find that π to the power of the natural log of π is equal to π to
the power of 0.23π‘ plus π. Well, π to the natural log of π
is simply π.
And then, we use the laws of
exponents to rewrite the right-hand side as π to the power of 0.23π‘ times π to
the power of π. But of course, π to the power of
π is itself a constant. Letβs call it π. So, π is equal to π times π to
the power of 0.23π‘. Now, weβre almost done. There were initially 245 rabbits on
the farm. In other words, when π‘ is equal to
zero, π is equal to 245. So, 245 equals π times π to the
power of zero because π to the power of zero is one. So, we find π to be equal to
245. And we found our function for
π. We can say that π of π‘ is equal
to 245π the power of 0.23π‘.
Now, itβs probably quite clear that
we didnβt need to set up and solve a differential equation here. We could have simply stated the
general form for an equation that describes exponential growth and decay. Itβs the general solution to dπ by
dπ‘ equals π times π for the initial value of π is equal to π nought. Itβs π of π‘ equals π nought
times π to the power of ππ‘. Essentially, π nought is the
initial value of the function, and π is the growth rate. Of course, itβs always useful to
know how to derive this formula, when necessary.
Weβll now look at an example of a
question where we can use the general solution to the differential equation without
fully deriving it.
A doctor injected a patient with 13
milligrams of radioactive dye that decays exponentially. A doctor injected a patient with 13
milligrams of radioactive dye that decays exponentially. After 12 minutes, there were 4.75
milligrams of dye remaining in the patientβs system. Which of the following is an
appropriate model for the situation? Is it A) π of π‘ equals 13 times
0.0805 to the power of π‘? B) π of π‘ equals 13 times π to
the power of 0.9195π‘. Is it C) π of π‘ equals 4.75 over
one plus 13 times π to the power of negative 0.83925π‘? Or D) π of π‘ equals 13 times π
to the power of negative 0.0839π‘.
Remember, we can model exponential
growth and decay by using the formula π of π‘ equals π of nought times π to the
power of ππ‘, where π nought is the initial value of π, and π is the rate of
growth or decay. So, letβs list what we actually
have. Since the answers in this question
are listed in terms of π, weβre going to change this equation to π of π‘ equals π
nought times π to the power of ππ‘, where π is the amount of dye thatβs actually
in the patientβs system after π‘ minutes. Now, weβre told that the doctor
injected the patient with 13 milligrams of radioactive dye. In other words, thatβs the initial
amount of dye in the patientβs system. So, π nought is 13. And we can therefore say that π of
π‘ is equal to 13 times π to the power of ππ‘.
But we also know that after 12
minutes, there were 4.75 milligrams of dye in the patientβs system. In other words, when π‘ is equal to
12, π is equal to 4.75. So, letβs substitute these values
into our function and solve for π. Thatβs 4.75 equals 13π to the
power of 12π. Weβll solve for π by first
dividing both sides by 13. And we see that 4.75 over 13 must
be equal to π to the power of 12π. Next, we find the natural logarithm
of both sides of this equation. This is useful as itβs the inverse
operation of raising something to the power of π. And so, these two operations will
essentially cancel one another out. In other words, the natural log of
π to the power of 12π is just 12π.
So, we find that the natural log of
4.75 over 13 equals 12π. Weβll solve now by dividing both
sides by 12. And so, π is equal to a 12 of the
natural log of 4.75 over 13. Typing this into our calculator and
we obtain π to be equal to negative 0.083900, and so on. Correct to three significant
figures, thatβs π is approximately equal to negative 0.0839. And substituting this back into our
equation gives us that π of π‘ equals 13 times π to the power of negative
0.0839π‘. And we see that that corresponds to
the answer D. π of π‘ is 13 times π to the
power of negative 0.0839π‘.
Once we have these models, we can
use them to make predictions about the future behavior. Letβs have a look at an example of
this form.
At the start of an experiment, a
scientist has a sample which contains 250 milligrams of a radioactive isotope. The radioactive isotope decays
exponentially at a rate of 1.3 percent per minute. a) Write the mass of the isotope
in milligrams, π, as a function of the time in minutes, π‘, since the start of the
experiment. And b) Find the half-life of the
isotope, giving your answer to the nearest minute.
Remember, we can model exponential
growth and decay using the formula π of π‘ equals π nought times π to the power
of ππ‘, where π nought is the initial value of π, and π is the rate of growth or
decay. In our case, there are 250
milligrams of the sample at the start of the experiment. So, π nought must be equal to
250. We said that π is the rate of
growth or decay. Since the isotope is decaying, this
value is going to be negative. And 1.3 percent as a decimal is
0.013. So, π is negative 0.013. We can therefore say that the
function that describes the mass in terms of minutes is π equals 250 times π the
power of negative 0.013π‘.
Weβre now going to look at part
b). Here, it helps us to understand the
definition of the phrase half-life. Itβs the time taken for the
radioactivity of the isotope to fall to half of its original value. In our case, thatβs half of
250. Thatβs 125 milligrams. We need to work out for what value
of π‘ π is equal to 125. We therefore set π equal to 125
and solve for π‘.
We see that 125 equals 250 times π
to the power of negative 0.013π‘. We then divide both sides of this
equation by 250. 125 divided by 250 is 0.5. So, we see that 0.5 equals π to
the power of negative 0.013π‘. We then take the natural logarithm
of both sides of this equation to get the natural log of 0.5 equals the natural log
of π to the power of negative 0.013π‘. But of course, the natural log of
π to the power of negative 0.013π‘ is just negative 0.013π‘. And so, our final step is to divide
both sides of this equation by negative 0.013. That gives us π‘ equals 53.319
minutes which, correct to the nearest minute, is 53 minutes.
In this video, weβve learned that
we can solve differential equations of the form dπ by dπ‘ equals ππ by writing
them in the form one over π dπ equals π dπ‘ and integrating both sides. When we do, we see that the
solution is of the form π equals π nought π to the ππ‘, where π of zero equals
π nought. We saw that whilst itβs useful to
be able to derive the formula, itβs absolutely fine to quote it and apply it to
problems involving exponential growth and decay.