#clancy1



#reactions
reaction0000001:
 {1.0}c3 = {1.0}c2
  cell*(-(b*c2) + a*c3)
reaction0000002:
 {1.0}c2 = {1.0}c1
  cell*(-(bin*c1) + ain*c2)
reaction0000003:
 {1.0}c1 = {1.0}o
  cell*(aa*c1 - bb*o)
reaction0000004:
 {1.0}o = {1.0}i
  cell*(-(ai*i) + bi*o)
reaction0000005:
 {1.0}c1 = {1.0}i
  cell*(aa*c1 - u*i)


#parameters
F = 96.485
R = 8.314
Temp = 310.0
ain = 2.172
bin = 1.077
ki = 140.0
ko = 5.4
vhold = -40.0
vtest = 0.0
cell = 1.0


#initial values
c1 = 0.0
c2 = 0.0
c3 = 1.0
i = 0.0
o = 0.06


#assignment rules
!F aa = 0.0655*exp(0.05547153*(-36.0 + v))
!F ik = Gk*(v - vk)*o
!F a = 0.0555*exp(0.05547153*(-12.0 + v))
!F b = 0.002357*exp(-0.036588*v)
!F Gk = 0.013500000000000002*ko^0.59
!F bi = 1.030072027074169*exp(0.000942*v)*(ko^(-1.0))^0.3
!F u = (ai*bb)/bi
!F bb = 0.0029357000000000003*exp(-0.02158*v)
!F vk = (R*Temp*log(ko/ki))/F
!F ai = 1.9755*exp(-0.02352*(25.0 + v))/ko
!F v = piecewise(vtest,50.0<time<=300.0,vhold)