#ratushny2



#reactions
___r1:
 $pool = {1.0}P
  ks * (k0 + pow(dspspkd / ka, h)) / (1 + pow(dspspkd / ka, h))
___r2:
 {1.0}P = $pool
  __RATE__ * P


#parameters
Kd = 1e-05
Ksp = 0.001
h = 2.0
k0 = 0.1
ka = 40.0
ks = 10.0
ku = 0.1
s = 1000.0
___r2___RATE__ = 0.1
univ = 1.0


#initial values
P = 0.0


#assignment rules
!F dspspkd = Kd / 2.0 * (1.0 + 2.0 * dspksp / Kd - pow(pow(1.0 + 2.0 * dspksp / Kd, 2.0) - 4.0 * pow(dspksp, 2.0) / pow(Kd, 2.0), 0.5))
!F dspksp = Ksp / 2.0 * (1.0 + (s + P * univ) / Ksp - pow(pow(1.0 + (s + P * univ) / Ksp, 2.0) - 4.0 * s * P * univ / pow(Ksp, 2.0), 0.5))