next up previous
Next: About this document ... Up: lab_template Previous: lab_template

Subsections


Applications of exponential and logarithmic functions

Purpose

The purpose of this lab is to use Maple to study applications of exponential and logarithmic functions. These are used to model many types of growth and decay, for example bacterial growth and radiaoctive decay. This lab also describes applications of exponential and logarithmic functions for heating and cooling and to medicine dosage

Background

Exponential growth

The simple model for growth is exponential growth, where it is assumed that $y'(t)$ is proportional to $y$. That is,

\begin{displaymath}\frac{dy}{dt} = ky\quad\hbox{where $k$\ is a positive constant.}\end{displaymath}

Separating the variables and integrating (see section 4.4 of the text), we have

\begin{displaymath}\int\frac{dy}{y} = \int k\,dt \end{displaymath}

so that

\begin{displaymath}
\ln \mid y \mid = kt + C
\end{displaymath}

In the case of exponential growth, we can drop the absolute value signs around $y$, because $y$ will always be a positive quantity. Solving for $y$, we obtain

\begin{displaymath}\mid y \mid = e^{kt + C} \end{displaymath}

which we may write in the form $y = Ae^{kt}$, where $A$ is an arbitrary positive constant.

Exponential Decay

In a sample of a radioactive material, the rate at which atoms decay is proportional to the amount of material present. That is,

\begin{displaymath}
\frac{dA}{dt} = -kA
\end{displaymath}

where $k>0$ is a constant. This is the same equation as in exponential growth, except that $-k$ replaces $k$. The solution is

\begin{displaymath}A(t) = A_0 e^{-kt} \end{displaymath}

where $A_0$ is a positive constant. Physically, $A_0$ is the amount of material present at $t=0$.

Radioactivity is often expressed in terms of an element's half-life. For example, the half-life of carbon-14 is 5730 years. This statement means that for any given sample of ${}^{14}\hbox{C}$, after 5730 years, half of it will have undergone decay. So, if the half-life is of an element Z is $c$ years, it must be that $e^{-kc}=\frac{1}{2}$, so that $kc=\ln 2$ and $k=\frac{\ln 2}{c}$.

Newton's law of cooling

What is usually called Newton's law of cooling is a simple model for the change in temperature of an object that is in contact with an environment at a different temperature. It says that the rate of change of the temperature of an object is proportional to the difference between the object's temperature and the temperature of the environment. Mathematically, this can be expressed as the differential equation

\begin{displaymath}\frac{dT}{dt} = -k \left( T-T_{\mathrm{out}} \right) \end{displaymath}

where $k$ is the constant of proportionality and $T_{\mathrm{out}}$ is the temperature of the environment. Using a technique called separation of variables it isn't hard to derive the solution

\begin{displaymath}T(t) = T_{\mathrm{out}} + \left( T_0 - T_{\mathrm{out}} \right)
e^{-kt} \end{displaymath}

where $T_0$ is the temperature of the object at $t=0$.

Effective medicine dosage

If a drug is administered to a patient intravenously, the concentration jumps to its highest level almost immediately. The concentration subsequently decays exponentially. If we use $C(t)$ to represent the concentration at time t, and $C_0$ to represent the concentration just after the dose is administered then our exponential decay model would be given by

\begin{displaymath}C(t) = C_0 e^{-kt} \end{displaymath}

A problem facing physicians is the fact that for most drugs, there is a concentration, $m$, below which the drug is ineffective and a concentration, $M$, above which the drug is dangerous. Thus the physician would like the have the concentration $C(t)$ satisfy

\begin{displaymath}m < C(t) < M \end{displaymath}

This means that the initial dose must not produce a concentration larger than $M$ and that another dose will have to be administered before the concentration reaches $m$.

Maple commands

The main functions you need are the natural exponential and natural logarithm. The Maple commands for these functions are exp and ln. Here are a few examples.
> f := x -> exp(-2*x);
> simplify(ln(3)+ln(9));
> ln(exp(x));
> simplify(ln(exp(x)),assume=real);
The assume=real is needed in the command above, because Maple usually works with complex variables. The command for base 10 logarithms is log10. Here are some examples. Note how Maple likes to convert base 10 logarithms to natural logarithms.
> log10(100);
> simplify(log10(100));
> simplify(log10(a)-log10(b),assume=real);

Sometimes you need to use experimental data to determine the value of constants in models. For example, suppose that for a particular drug, the following data were obtained. Just after the drug is injected, the concentration is 1.5 mg/ml (milligrams per milliliter). After four hours the concentration has dropped to 0.25 mg/ml. From this data we can determine values of $C_0$ and $k$ as follows. The value of $C_0$ is the initial concentration, so we have

\begin{displaymath}C_0 = 1.5\end{displaymath}

To find the value of $k$ we need to solve the equation

\begin{displaymath}0.25 = 1.5 e^{-4k} \end{displaymath}

which we get by plugging in $t=4$ and using the data $C(4)=0.25$. Maple commands for solving for $k$ and defining and plotting the function $C(t)$ are shown below.
> k1 := solve(0.25=1.5*exp(-4*k),k);
> C1 := t -> 1.5*exp(-k1*t);
> plot(C1(t),t=0..6);

Exercises

  1. Suppose that the population of a certain bacteria can be modeled by an exponential function. In a particular experiment, the number of bacteria was $10,000$ at $t=0$. Four hours later, the number of bacteria was $46,000$.
    1. Find the growth constant $k$.
    2. Predict the number of bacteria that would have been present after $6$ hours.

  2. A thermometer is taken from a room at $68 \, ^{\circ} \mathrm{F}$ to the outdoors where the temperature is $13 \, ^{\circ} \mathrm{F}$:
    1. Solve for k if the reading of the thermometer drops to $46 \,^{\circ} \mathrm{F}$ after one minute.
    2. Find the reading of the thermometer after five minutes.

  3. Suppose that for a certain drug, the following results were obtained. Immediately after the drug was administered, the concentration was 13 mg/ml. Five hours later, the concentration had dropped to 3.7 mg/ml.
    1. Solve for the value of $k$ for this drug.
    2. Find the concentration of the drug after eight hours.
    3. Find the number of hours needed to drop the concentration of the drug to effectively $zero$ (a concentration $\leq 0.001$).


next up previous
Next: About this document ... Up: lab_template Previous: lab_template
Jane E Bouchard
2009-02-18