No Title                     |POLVER05_3
d(Cvm)/d(t) =   (vm/Vm)*(Cvc-Cvm)-kv*Cvm*Cs-kas*Cvm*Cas   #Venom concentration in muscle compartment
d(Cam)/d(t) = (vm/Vm)*(Cac-Cam)-ka*(Cam*Cs-Cas/Ka)-kas*Cam*Cs #Antivenin concentration in muscle compartment
d(Cavm)/d(t) =    (vm/Vm)*(Cavc-Cavm)+(Ct)*(kas*fas*Cvm+kvs*fvs*Cam)    #Venom-antivenin product concentration in muscle compartment
d(fvs)/d(t) =     kv*fs*Cvm-kvs*fvs*Cam     #Fraction of sites occupied by venom in muscle compartment
d(fas)/d(t) =     ka*(Cam*fs-fas/Ka)-kas*fas*Cvm     #Fraction of sites occupied by antivenin in muscle compartment
d(Cvc)/d(t) = (vm/Vc)*(Cvm-Cvc)-(ve/Vc)*Cvc-kpc*Cac*Cvc #Venom concentration in central compartment
d(Cac)/d(t) = (vm/Vc)*(Cam-Cac)-(ve/Vc)*Cac-kpc*Cac*Cvc #Antivenin concentration in central compartment
d(Cavc)/d(t) =     (vm/Vc)*(Cavm-Cavc)-(ve/Vc)*Cavc+kpc*Cac*Cvc     #Venom-antivenin product concentration in central compartment
fs =     1-fvs-fas     #Fraction of free sites in muscle compartment
Ka =     ka/k_a     #Equilibrium constant for antivenin adsorption reaction in muscle compartment
ka = 3e6 #Forward rate constant for Antivenin Adsorption in muscle compartment (L/mol hr)
k_a = 4.5 #Backward rate constant for antivenin adsorption in muscle compartment (1/hr)
kv = 5e8 #Rate constant for adsorption of venom in muscle compartment (L/mol hr)
kvs = 1.8e7 #Rate constant for reaction of antivenom with venom on site in muscle compartment (L/mol hr)
kas = 14 #Rate constant for reaction of venom with antivenin on site in muscle compartment (L/mol hr)
Cas =     fas*Ct     #Concentration of sites occupied by antivenin 
Cs =     fs*Ct     #Concentration of free sites 
Ct = 2.2e-9 #Total concentration of sites
vm = 11.16 #Blood flow from muscle compartment to central compartment (L/hr)
Vm = 1.62 #Tissue water volume of muscle compartment (L)
Vc = .798 #Tissue water volume of central compartment (L)
ve = 7.5e-3 #Elimination flow rate (L/hr)
kpc = 8.6e4 #Rate constant for formation of AV product in central compartment (L/mol hr)
t(0)=0
Cvm(0)=0
Cam(0)=0
Cavm(0)=0
fvs(0)=0
fas(0)=0
Cvc(0)=1.35e-8
Cac(0)=2.5e-8
Cavc(0)=0
t(f)=8