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. In 1935 Charles F. Richter of Cal Tech developed a scale for measuring the magnitude of earthquakes. The Richter Scale formula is given by

    \begin{displaymath}M = \log_{10} \left( \frac{x}{c} \right) \end{displaymath}

    where $M$ is the magnitude of the earthquake, $x$ is the amplitude of the largest seismic wave as measured on a standard seismograph 100 kilometers from the epicenter and $c$ is the amplitude of a reference earthquake of amplitude 1 micron on a standard seismograph at the same distance from the epicenter.
    1. When the amplitude of an earthquake is tripled, by how much does the magnitude increase?
    2. In 1989, the San Francisco Bay area suffered severe damage from an earthquake of magnitude 7.1. However, the damage was not nearly as extensive as that caused by the Alaska earthquake of 1964, which had magnitude 8.6. What is the ratio of the amplitude of the 1964 quake to the 1989 quake?

    3. The 1906 San Francisco earthquake had an intensity 16 times that of the 1989 quake. What was the Richter scale magnitude of the 1906 quake?

  2. A metal plate is heated to a temperature of $425 \, ^{\circ} \mathrm{C}$ then allowed to cool in a room where the temperature is kept at $25 \, ^{\circ}
\mathrm{C}$. The plate cools so that the difference between the plate temperature and the room temperature decreases by half of the difference every 1/4 hour. What is the temperature of the plate after 1/4 hour? After 1/2 hour? After 1 hour?

  3. Suppose that for a certain drug, the following results were obtained. Immediately after the drug was administered, the concentration was 8 mg/ml. Four hours later, the concentration had dropped to 2.3 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
2002-12-03