{UTF-8}
input=
	rate constant in*PULSE( 0 , duration )
	~	
	~		|

duration=
	20
	~	 [0,75,1]
	~		|

absorption= INTEG (
	ka1-ka2,
		0)
	~	
	~		|

drug level= INTEG (
	input-output-ka1+ka2,
		0)
	~	
	~		|

ka1=
	k1a*drug level
	~	
	~		|

k1a=
	0.1
	~	 [0,1,0.1]
	~		|

k2a=
	0.1
	~	 [0,1,0.1]
	~		|

ka2=
	k2a*absorption
	~	
	~		|

rate constant in=
	0.2
	~	 [0,1,0.1]
	~		|

output=
	rate constant*drug level
	~	
	~		|

rate constant=
	0.2
	~	 [0,1]
	~		|

therapeutic level=
	2
	~	 [0,10,1]
	~		|

toxic level=
	4
	~	 [0,10,1]
	~		|

********************************************************
	.Control
********************************************************~
		Simulation Control Parameters
	|

FINAL TIME  = 100
	~	Hour
	~	The final time for the simulation.
	|

INITIAL TIME  = 0
	~	Hour
	~	The initial time for the simulation.
	|

SAVEPER  = 
        TIME STEP 
	~	Hour [0,?]
	~	The frequency with which output is stored.
	|

TIME STEP  = 1
	~	Hour [0,?]
	~	The time step for the simulation.
	|

\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*View 1
$192-192-192,0,Comic Sans MS|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|96,96,100,0
10,1,drug level,416,136,74,24,3,131,0,1,0,3,0,0,0-0-255,0-0-0,|12||0-0-0
12,2,48,197,136,10,8,0,3,0,0,-1,0,0,0
1,3,5,1,4,0,0,22,0,0,0,-1--1--1,,1|(311,136)|
1,4,5,2,100,0,0,22,0,0,0,-1--1--1,,1|(237,136)|
11,5,48,274,136,6,8,34,3,0,0,1,0,0,0
10,6,input,274,157,21,13,40,3,0,0,-1,0,0,0
12,7,48,661,136,10,8,0,3,0,0,-1,0,0,0
1,8,10,7,4,0,0,22,0,0,0,-1--1--1,,1|(613,136)|
1,9,10,1,100,0,0,22,0,0,0,-1--1--1,,1|(527,136)|
11,10,48,570,136,6,8,34,3,0,0,1,0,0,0
10,11,output,570,157,26,13,40,3,0,0,-1,0,0,0
1,12,1,10,1,0,0,0,0,64,0,-1--1--1,,1|(562,117)|
12,13,0,419,349,304,139,3,188,0,0,1,0,0,0
graph
10,14,therapeutic level,87,36,20,20,2,3,0,1,1,3,0,0,0-128-0,0-0-0,|12||0-0-0
10,15,toxic level,746,31,20,20,2,3,0,1,1,3,0,0,255-0-0,0-0-0,|12||0-0-0
10,16,rate constant in,128,159,14,14,2,131,0,0,1,0,0,0
1,17,16,5,1,0,0,0,0,64,0,-1--1--1,,1|(153,111)|
10,18,rate constant,671,77,16,16,2,131,0,0,1,0,0,0
1,19,18,10,0,0,0,0,0,64,0,-1--1--1,,1|(622,105)|
10,20,absorption,414,40,62,19,3,131,0,0,0,0,0,0
1,21,23,20,4,0,0,22,0,0,0,-1--1--1,,1|(379,68)|
1,22,23,1,100,0,0,22,0,0,0,-1--1--1,,1|(379,101)|
11,23,268,379,84,8,6,33,3,0,0,4,0,0,0
10,24,ka1,401,84,14,13,40,3,0,0,-1,0,0,0
1,25,27,1,4,0,0,22,0,0,0,-1--1--1,,1|(443,100)|
1,26,27,20,100,0,0,22,0,0,0,-1--1--1,,1|(443,68)|
11,27,284,443,83,8,6,33,3,0,0,4,0,0,0
10,28,ka2,467,83,16,13,40,3,0,0,-1,0,0,0
1,29,1,23,1,0,0,0,0,64,0,-1--1--1,,1|(369,85)|
1,30,20,27,1,0,0,0,0,64,0,-1--1--1,,1|(489,61)|
10,31,k1a,267,41,16,16,2,131,0,0,1,0,0,0
10,32,k2a,563,35,15,15,2,131,0,0,1,0,0,0
1,33,31,23,0,0,0,0,0,64,0,-1--1--1,,1|(319,60)|
1,34,32,27,0,0,0,0,0,64,0,-1--1--1,,1|(507,56)|
10,35,duration,66,136,15,15,2,131,0,0,1,0,0,0
1,36,35,5,1,0,0,0,0,64,0,-1--1--1,,1|(167,106)|
///---\\\
:GRAPH GRAPH
:TITLE GRAPH
:SCALE
:VAR drug level
:Y-MIN 0
:LINE-WIDTH 3
:VAR toxic level
:LINE-WIDTH 3
:VAR therapeutic level
:LINE-WIDTH 3
:L<%^E!@
1:Current.vdf
9:Current
23:0
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
4:Time
5:duration
24:0
25:100
26:100