next up previous
Next: Exercises Up: Background Previous: Maximum time between doses

Equal, regularly-spaced doses

We next consider what happens if equal doses of the drug are given at regular time intervals. Recall that a drug has a maximum safe concentration, M, and a minimum effective concentration, m. We say that a treatment program of equal, regularly-spaced doses is safe and effective if the concentration C(t) of the drug satisfies


during the treatment.

In the first part of this lab, we presented the expression


for the concentration of the drug after the first dose. This expression is valid as long as only a single dose is given. However, suppose that at t=L a second dose is given and that the amount of the drug administered is the same as the first dose. According to our model, the concentration will jump immediately by an amount equal to tex2html_wrap_inline110 when the second dose is given. However, when the second dose is given, there is still some of the drug in the bloodstream remaining from the first dose. This means that to compute the concentration just after the second dose, we have to add the value tex2html_wrap_inline110 to the concentration remaining from the first dose. During the time between the second and third doses, the concentration decays exponentially from this value. To find the concentration after the third dose, we would have to repeat this process, but now we have contributions from the first and second doses to include.

Below, we describe how one can find the concentration at any later time, assuming that the doses are equal and that the time between doses is fixed, but it turns out to be rather complicated. To help you visualize what is going on, a Maple procedure called drug has been written. If you give this procedure values for tex2html_wrap_inline110 , the concentration jump caused by each dose, L, the time between doses, and k, the decay rate for the drug, it returns the concentration as a function of t. Before you can use this command, you must first load it with the first command shown below. If you get an error message, get help from your IA. Once you've successfully loaded the drug command, the second command below will bring up the help page.

  > readlib(drug):

  > ?drug

To see how to use the drug command, consider the drug from the exercises in part one of this lab. Use the values tex2html_wrap_inline170 and L=6. That is, the initial dose produces a concentration of 7 mg/ml and the time between doses is 6 hours. Recall that the maximum safe dose is M=8 and the minimum effective dose is m=1.8 Then the following commands will define the concentration C(t), evaluate it just after the first dose and just before and just after the second dose, and, finally, plot the concentration over a period of 24 hours.

  > C := drug({C0=7,L=6,k=k2});


  > C(0);


  > C(5.999);


  > C(6);


  > plot(C(t),t=0..24);

Notice that the concentration just after a dose increases with each successive dose. Notice also that the concentration soon exceeds the maximum safe concentration M=8. This means that the dose is too high, and this is not a safe treatment program.

You can go back and plot the concentration for a longer time interval. You should notice that the values of the peak concentrations just after a dose rise rapidly for the first few doses and then appear to level off. It can be shown that this always happens, so the maximum concentration eventually levels off as more doses are given. The limiting value for an infinite number of doses can be calculated explicitly, using information on geometric series, but that is a topic for Calc 3. To help you understand what is involved in the drug program, we describe how to obtain formulas for the concentration after the first few doses below.

We can calculate the concentration just before the second dose is administered by setting t=L in our equation


to get


Now, when the second dose is administered the concentration jumps by an increment tex2html_wrap_inline110 so that the concentration just after the second dose is given is


The concentration then decays from this value according to our exponential decay rule, but with a slight twist. The twist is that the ``inital'' concentration is at t=L, instead of the more familiar case of t=0. One way to handle this is to write the exponential term as


so that at t=L, the exponent is 0. If we do this, then we can write the concentration as a function of time as


This function is only valid after the second dose is administered and before the third dose is given.

Now, suppose that a third dose of the drug is given at t=2L. The concentration just before the third dose is given is C2(2L), which is


which we can also write as


When the third dose is given, the concentration would jump again by tex2html_wrap_inline110 and the concentration just after the third dose would be


Finally, an expression for the concentration that would be valid after the third dose is administered would be


This process can be continued to get expressions for the concentration for any number of equal, regularly-spaced doses.

next up previous
Next: Exercises Up: Background Previous: Maximum time between doses

William W. Farr
Fri Feb 7 13:14:44 EST 1997