begin name
mayya
end name

begin reactions
v1 {1.0}sn = {0.0684492}sc
v2 {0.0684492}sc = {1.0}sn
v3 {1.0}sc = {1.0}psc
v4 {1.0}psc = {0.5}ps2c
v5 {1.0}ps2c = {2.0}psc
v6 {0.0684492}ps2c = {1.0}ps2n
v7 {1.0}psn = {0.5}ps2n
v8 {1.0}ps2n = {2.0}psn
v9 {1.0}psn = {1.0}sn
end reactions

begin rate equations
v1 = (sexpmax*sn[t])/(Ksexp + sn[t])
v2 = (simpmax*sc[t])/(Ksimp + sc[t])
v3 = (Kcatp*sk*sc[t])/(Kmp + sc[t])
v4 = Kfd*psc[t]^2
v5 = Krd*ps2c[t]
v6 = (ps2impmax*ps2c[t])/(Kps2imp + ps2c[t])
v7 = Kfd*psn[t]^2
v8 = Krd*ps2n[t]
v9 = (Kcatde*SP*psn[t])/(Kmde + psn[t])
end rate equations

begin 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
defaultcompartment = 1.0
end parameters

begin initial conditions
ps2c[0] = ps2ci
ps2n[0] = ps2ni
psc[0] = psci
psn[0] = psni
sc[0] = sci
sn[0] = sni
end initial conditions

begin initial values
ps2ci = 0.0
ps2ni = 0.0
psci = 0.0
psni = 0.0
sci = 1.0
sni = 0.0
end initial values

begin assignment rules
sk := Piecewise[{{0, t < lag || t > lag + period}}, ampl*Sin[(Pi*(-lag + t))/period]]
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units