#olivier

FIX: p s

#reactions
v_1:
 {1.0}s = {1.0}a
  piecewise(V1 * s / Ks * (1 - a / (s * Keq1)) * (1 + s / Ks + a / Ka)**(n - 1) / ((1 + s / Ks + a / Ka)**n + L0 * ((1 + p / Kp) / (1 + Cp * p / Kp))**n), gt(eq, 0), V1 * s / s05 * (1 - a / (s * Keq1)) * (s / s05 + a / a05)**(h - 1) / ((s / s05 + a / a05)**h + (1 + (p / p05)**h) / (1 + alfa * (p / p05)**h)))
v_2:
 {1.0}a = {1.0}b
  V2*(a-b/Keq2)/(a+K2a*(1+b/K2b))
v_3:
 {1.0}b = {1.0}p
  V3*(b-p/Keq3)/(b+K3b*(1+p/K3p))


#parameters
Cp = 0.026
K2a = 1.0
K2b = 1.0
K3b = 1.0
K3p = 1.0
Ka = 1000.0
Keq1 = 400.0
Keq2 = 10.0
Keq3 = 10.0
Kp = 0.1
Ks = 0.01
L0 = 10000.0
V1 = 200.0
V2 = 100.0
V3 = 100.0
a05 = 10000.0
alfa = 0.01
eq = 1.0
h = 4.0
n = 4.0
p05 = 1.0
s05 = 1.0
p = 1.0
s = 1.0
default_compartment = 1.0


#initial values
a = 10.0
b = 10.0


#assignment rules