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. The loudness of sound is measured in decibels in honor of Alexander Graham Bell (1847-1922), inventor of the telephone. If the variation in pressure is $P$ pounds per square inch, then the loudness in decibels is

    \begin{displaymath}L = 20 \log_{10} (121.3P) \end{displaymath}

    Find how many times greater is the variation in pressure caused by a rock band at 115 decibels than a soft rock concert at 85 decibels.

  2. All living things contain carbon-12, which is stable, and carbon-14, which is radioactive. While a plant or animal is alive, the ratio of these two isotopes of carbon remains unchanged since the carbon-14 is constantly renewed; after death, no more carbon-14 is absorbed. The half life of carbon-14 is 5730 years. Human hair from a grave in Africa proved to have 57% of the carbon-14 of living tissue. When was the body buried?

  3. A 5-lb roast initially at a temperature of $50 \, ^{\circ} \mathrm{F}$ is put into an oven at $375 \, ^{\circ} \mathrm{F}$ when $t=0$. It is found that the temperature of the roast is $125 \, ^{\circ} \mathrm{F}$ 75 minutes later. How many more minutes will it take for the roast to be medium rare, a temperature of $150 \, ^{\circ} \mathrm{F}$ ?

  4. 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. Determine the value of $k$ for this drug. What is the concentration of the drug 8 hours later?


next up previous
Next: About this document ... Up: lab_template Previous: lab_template
William W. Farr
2003-02-09