See attachement.

profilevx6
lab_2_cipro1.doc

AMTH 518 Lab 2: DDS (Cipro)

Scenario

Consider a drug dosage problem for combating "anthrax". The doctor prescribes 500 mg of Cipro as a daily dosage. We assume that the kidney has diluted the drug by 45% every day. Model as a DDS and determine the long-term behavior of the system.

Problem Statement: Determine the long-term amount of Cipro in the body.

Assumptions: Nothing changes about the way the body processes Cipro for the duration of the model.

Model:

D(n) = Amount of Cipro in the body after day n.

D(n+1) = D(n)-.45*D(n)+500, D(0)=0

or

D(n+1) = .55 D(n) + 500, D(0)=0

1. Open a new worksheet

2. Enter your name and a brief problem description.

3. Enter the model

4. Create column headings for n, D(n) (assume we put n in cell A10, and D(n) in cell B10.

5. In the next row, enter the initial condition D(0)=0: put 0 in cell A11 and 0 in Cell B11.

6. In cell A12 put =a11+1

7. In cell B12 put =.55*b11+500

8. Copy those formulas into the remaining cells below until you get to a value of n that you desire.

This looks like:

image4.wmf

Drug Problem

0

500

1000

1500

0

5

10

15

20

time

D(n)

Next, we might want to view a graph.

Highlight the columns and go to CHART WIZARD

Select SCATTERPLOT (x vs y)

image1We see a stable equilibrium at about 1111.

What is the true equilibrium value?

EXTENSION#1:

If you wanted to specify an equilibrium amount , E, in the blood, derive an equation that tells how many mg, K, would need to be in each daily dose?

HINT: Set D(n+1) = D(n) = E in the iteration equation, and replace the “500” with K. Solve for K in terms of E.

EXTENSION#2:

If you wanted to maintain an equilibrium of 1000 mg in the blood by giving two doses per day, how many mg would need to be in each dose?

HINT: Look at what happens if you start with 1000 mg in the blood with NO ADDITIONAL SHOTS. Does this look like exponential decay? What is the equation? If 0.55 of the original amount is left at the end of 24 hours, how much would be left at the end of just 12 hours? Now, set D(n+1/2) = D(n) = 1000 (so we get 1000 mg equilibrium with two doses per day), and solve for the dose (call it K).

� EMBED Excel.Sheet.8 ���

image2.wmf

Drug Problem

0

500

1000

1500

0

5

10

15

20

time

D(n)

image3.png

_1155371825.xls

Chart1

0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
time
D(n)
Drug Problem
0
500
775
926.25
1009.4375
1055.190625
1080.35484375
1094.1951640625
1101.8073402344
1105.9940371289
1108.2967204209
1109.5631962315
1110.2597579273
1110.64286686
1110.853576773
1110.9694672252
1111.0332069738

Sheet1

Name:
Drug Problem
D(n+1)=.55*D(n)+500, D(0)=0
n d(n)
0 0
1 500
2 775
3 926.25
4 1009.4375
5 1055.190625
6 1080.35484375
7 1094.1951640625
8 1101.8073402344
9 1105.9940371289
10 1108.2967204209
11 1109.5631962315
12 1110.2597579273
13 1110.64286686
14 1110.853576773
15 1110.9694672252
16 1111.0332069738

Sheet1

0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
time
D(n)
Drug Problem
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0

Sheet2

Sheet3