#ratushny3



#reactions
___r1:
 $pool = {1.0}P1
  ks * (k0 + pow(dsp1p2kd / ka, h)) / (1 + pow(dsp1p2kd / ka, h))
___r2:
 {1.0}P1 = $pool
  __RATE__ * P1
___r3:
 $pool = {1.0}P2
  ks * (k0 + pow(dsp1p2kd / ka, h)) / (1 + pow(dsp1p2kd / ka, h))
___r4:
 {1.0}P2 = $pool
  __RATE__ * P2
___r5:
 $pool = {1.0}Target
  ks * (k0 + pow(dsp1p2kd / ka, h)) / (1 + pow(dsp1p2kd / ka, h))
___r6:
 {1.0}Target = $pool
  __RATE__ * Target


#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
___r4___RATE__ = 0.1
___r6___RATE__ = 0.1
univ = 1.0


#initial values
P1 = 0.0
P2 = 0.0
Target = 0.0


#assignment rules
!F dsp1p2kd = Kd / 2.0 * (1.0 + (dsp1ksp + P2 * univ) / Kd - pow(pow(1.0 + (dsp1ksp + P2 * univ) / Kd, 2.0) - 4.0 * dsp1ksp * P2 * univ / pow(Kd, 2.0), 0.5))
!F dsp1ksp = Ksp / 2.0 * (1.0 + (s + P1 * univ) / Ksp - pow(pow(1.0 + (s + P1 * univ) / Ksp, 2.0) - 4.0 * s * P1 * univ / pow(Ksp, 2.0), 0.5))