Parameter Values:
nr = 30.000000
kao_f = 1.000000
ka_f = 0.003300
nc = 2.000000
kr_f = 0.500000
kmdiff_r = 0.010000
kc_f = 0.050000
t = 0.000000
ng = 2.000000
kmdiff_f = 1.000000
kecd = 0.005000
ko = 0.050000
ko_r = 1.000000
kecdiff = 1.000000
kao_r = 1.000000
ka = 0.250000
ka_r = 1.000000
kb = 0.000100
np = 10.000000
kr_r = 1.000000
time = 0.000000
kc_r = 1.000000
kd = 0.007500
ko_f = 0.033000

Initial State Vector:
CI = 0.000000
LacI = 0.000000
TetR = 0.000000
P0 = 2.000000
P1 = 2.000000
P2 = 2.000000

Initial Reaction Rate Array:
0.000000
0.000000
0.000000
0.049749
0.049749
0.049749

Reaction Rate Equation Array:
kd * CI
kd * LacI
kd * TetR
P0 * ko * ko_f / ko_r * nr / ( 1 + ko_f / ko_r * nr + ( kr_f / kr_r * CI ) ^ nc )
P1 * ko * ko_f / ko_r * nr / ( 1 + ko_f / ko_r * nr + ( kr_f / kr_r * LacI ) ^ nc )
P2 * ko * ko_f / ko_r * nr / ( 1 + ko_f / ko_r * nr + ( kr_f / kr_r * TetR ) ^ nc )

Reactant Matrix:
1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 

Delta Matrix:
-1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 -1.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 -1.000000 0.000000 0.000000 0.000000 
0.000000 10.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 0.000000 10.000000 0.000000 0.000000 0.000000 
10.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
