next up previous
Next: About this document ... Up: Labs and Projects for Previous: Labs and Projects for

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$.

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);

\begin{maplelatex}
\begin{displaymath}
{f} := {x} \rightarrow {\rm e}^{(\, - 2\,{x}\,)}
\end{displaymath}\end{maplelatex}
  > 
simplify(ln(3)+ln(9));

\begin{maplelatex}
\begin{displaymath}
3\,{\rm ln}(\,3\,)
\end{displaymath}\end{maplelatex}
  > 
ln(exp(x));

\begin{maplelatex}
\begin{displaymath}
{\rm ln}(\,{\rm e}^{{x}}\,)
\end{displaymath}\end{maplelatex}
  > 
simplify(ln(exp(x)),assume=real);

\begin{maplelatex}
\begin{displaymath}
{x}
\end{displaymath}\end{maplelatex}
The assume=real is needed in the command above, because Maple usually works with complex variables.
  > 
solve(exp(-3*x)=0.5,x);

\begin{maplelatex}
\begin{displaymath}
.2310490602
\end{displaymath}\end{maplelatex}
  > 
plot(log[10](x),x=0..100);

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);

\begin{maplelatex}
\begin{displaymath}
{\it k1} := .4479398673
\end{displaymath}\end{maplelatex}
  > 
C1 := t -> 1.5*exp(-k1*t);

\begin{maplelatex}
\begin{displaymath}
{\it C1} := {t} \rightarrow 1.5\,{\rm e}^{(\, - {\it k1}\,{t}\,)}
\end{displaymath}\end{maplelatex}
  > 
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 ( 1 micron is 0.001 mm) on a standard seismograph at the same distance from the epicenter.
    1. 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 great quake of 1906, which has been estimated to have had the magnitude 8.3. What is the ratio of the amplitude of the 1906 quake to the 1989 quake? The largest earthquake magnitude ever measured was $8.9$ for an earthquake in Japan in 1933. Determine the ratio of the amplitude of this earthquake to that of the 1906 San Francisco earthquake.

  2. Worldwide use of the Internet is increasing at an exponential rate. If traffic on the Internet is doubling every 150 days, find the exponential growth rate constant.

  3. Suppose that the last bit of ice in a picnic cooler has melted. How long will it take for the temperature inside to reach $50
\, ^{\circ} \mathrm{F}$? Use Newton's law of cooling to model this, using $T_{\mathrm{out}} = 90 \, ^{\circ} \mathrm{F}$, $T_0 = 32 \, ^{\circ}
\mathrm{F}$ and $k = 0.005 \, \mathrm{min}^-1$.

  4. The worst nuclear accident in history happened in 1986 at the Chernobyl nuclear plant near Kiev in the Ukraine. An explosion destroyed one of the plant's four reactors, releasing large amount of radioactive isotopes into the atmosphere.

    Consider 10 grams of the plutonium isotope Pu-239 released in the Chernobyl nuclear accident. This isotope has a half-life of $24,360$ years. How long will it take for the 10 grams to decay to 1 gram?


next up previous
Next: About this document ... Up: Labs and Projects for Previous: Labs and Projects for
Jane E Bouchard
2001-04-16