#wu1



#reactions
reaction_0000005:
 $pool = {1.0}c
  cell * fcyt * (jmem + jer)
reaction_0000006:
 $pool = {1.0}cer
  -fer * vcytver * jer * cell


#parameters
alpha = 4.5e-06
cm = 5300.0
dact = 0.35
dinact = 0.4
dip3 = 0.5
fcyt = 0.01
fer = 0.01
gca = 1000.0
gir = 5.0
gk = 1400.0
gkatp = 500.0
gkca = 900.0
ip3 = 0.0
kd = 0.3
kpmca = 0.2
kserca = 0.4
pleak = 0.0005
r = 0.14
sa = 0.1
sm = 12.0
sn = 5.0
taua = 300000.0
taun = 16.0
vca = 50.0
vcytver = 5.0
vir = -75.0
vk = -75.0
vm = -20.0
vn = -16.0
cell = 1.0


#initial values
c = 0.1
cer = 100.0
V = -60.0
a = 0.46
n = 0.01
rIR = 0.282
nIR = 0.008


#assignment rules
!F nIRinf = 1.0 / (1.0 + betaIRn / alphaIRn)
!F betaIRn = 0.00035 * exp(0.07 * (V + 25.0))
!F alphaIRn = 0.09 / (1.0 + exp(0.11 * (V + 100.0)))
!F oinf = c / (dact + c) * (ip3 / (dip3 + ip3)) * (dinact / (dinact + c))
!F jleak = pleak * (cer - c)
!F jserca = kserca * c
!F jmem = -(alpha * ica + kpmca * c)
!F ikca = gkca * w * (V - vk)
!F ainf = 1.0 / (1.0 + exp((r - c) / sa))
!F ninf = 1.0 / (1.0 + exp((vn - V) / sn))
!F ica = gca * minf * (V - vca)
!F ikatp = gkatp * a * (V - vk)
!F w = pow(c, 5.0) / (pow(c, 5.0) + pow(kd, 5.0))
!F ik = gk * n * (V - vk)
!F minf = 1.0 / (1.0 + exp((vm - V) / sm))
!F tauIRr = 1.0 / (alphaIRr + betaIRr)
!F rIRinf = 1.0 / (1.0 + betaIRr / alphaIRr)
!F betaIRr = 0.15 / (1.0 + exp(-0.05 * (V + 120.0)))
!F alphaIRr = 30.0 / (1.0 + exp(0.04 * (V + 230.0)))
!F tauIRn = 1.0 / (alphaIRn + betaIRn)
!F iir = gir * nIR * rIR * (V - vir)
!F jer = jleak + jip3 - jserca
!F jip3 = oinf * (cer - c)