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. An example of common logarithms is the decibel scale, particularly used for measuring loudness. (The decibel unit is named in honor of Alexander G. Bell (1847-1922), inventor of the telephone.) If $I$ is the intensity of sound in watts per square meter, the decibel level of the sound is

    \begin{displaymath}d = 10 \log_{10} (\frac{I}{I_0}) \end{displaymath}

    where $I_0$ is an intensity corresponding roughly to the faintest sound that can be heard.
    When tuning the rock band's equipment before the concert in a big concert hall, an audio engineer finds that in order to maintain appropriate loudness in this hall, he needs to increase the power of the amplifiers in comparison with the level used for the previous concert in a hall of smaller size.
    1. What does the doubling the intensity add to the level of loudness in decibels?

    2. By what factor $k$ will the engineer have to multiply the intensity of the sound to add 20 to the sound level for the next concert of the band on the stadium?

  2. The cry of a blue whale can be heard nearly 500 miles away, reaching an intensity of (6.3 x $\displaystyle 10^{18}$)$I_0$. Using the same formula for loudness given in exercise 1, determine the decibel level of this sound.

  3. A thermometer registered $-20 \, ^{\circ} \mathrm{C}$ outside and then was brought into the house where the temperature was $24 \, ^{\circ} \mathrm{C}$. After 5 minutes, it registered $0 \, ^{\circ} \mathrm{C}$. When will it register $20 \, ^{\circ} \mathrm{C}$?

  4. Suppose that for a certain drug, the following results were obtained. Immediately after the drug was administered, the concentration was 6.9 mg/ml. Four hours later, the concentration had dropped to 2.4 mg/ml. Determine the value of $k$ for this drug.


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