#cloutier1



#reactions
ROS_1:
 $pool = {1.0}ROS
  Neuron * V1(k1, S1, dalphasyn, alpha_syn, kalphasyn)
ROS_2:
 {1.0}ROS = $pool
  Neuron * V2_4(k2, ROS, S2_4)
aSyn_1:
 $pool = {1.0}alpha_syn
  Neuron * V3(k3, ROS, S2_4)
aSyn_2:
 {1.0}alpha_syn = $pool
  Neuron * V2_4(k4, alpha_syn, S2_4)


#parameters
S2_4 = 1.0
dalphasyn = 15.0
k1 = 17280.0
k2 = 17280.0
k3 = 0.168
k4 = 0.168
kalphasyn = 8.5
Neuron = 1.0


#initial values
ROS = 1.0
alpha_syn = 1.0


#assignment rules
!F S1 = piecewise(2.6, and(gt(time, 10.0), lt(time, 150.0)), 0.0)