#mayya



#reactions
v_1:
 {1.0}sn = {0.0684492}sc
  sexpmax * sn / (Ksexp + sn)
v_2:
 {0.0684492}sc = {1.0}sn
  simpmax*sc/(Ksimp + sc)
v_3:
 {1.0}sc = {1.0}psc
  Kcatp*sk*sc/(Kmp + sc)
v_4:
 {1.0}psc = {0.5}ps2c
  Kfd * psc**2
v_5:
 {1.0}ps2c = {2.0}psc
  Krd*ps2c
v_6:
 {0.0684492}ps2c = {1.0}ps2n
  ps2impmax*ps2c /(Kps2imp + ps2c)
v_7:
 {1.0}psn = {0.5}ps2n
  Kfd * psn**2
v_8:
 {1.0}ps2n = {2.0}psn
  Krd*ps2n
v_9:
 {1.0}psn = {1.0}sn
  Kcatde*SP*psn/(Kmde + psn)


#parameters
Kcatde = 1.0
Kcatp = 1.0
Kfd = 0.6
Kmde = 2.0
Kmp = 4.0
Kps2imp = 0.3
Kpsexp = 0.6
Kpsimp = 3.0
Krd = 0.03
Ksexp = 0.6
Ksimp = 3.0
SP = 0.05
SVnm = 1.0
Vc = 0.935
Vn = 0.064
ampl = 0.01
lag = 500.0
period = 2000.0
ps2impmax = 0.045
psexpmax = 0.06
psimpmax = 0.003
sexpmax = 0.06
simpmax = 0.003
default_compartment = 1.0


#initial values
ps2c = 0.0
ps2n = 0.0
psc = 0.0
psn = 0.0
sc = 1.0
sn = 0.0


#assignment rules
!F sk = piecewise(0,time < lag, 0, time > lag + period, ampl*sin(pi/period*(time-lag)))