#bertram1



#reactions
reaction_0000001:
 $pool = {1.0}c
  cell * f * (jertot + jmemtot)
reaction_0000002:
 $pool = {1.0}cer
  -fer * sigmav * jertot * cell
reaction_000003:
 $pool = {1.0}cAMP
  cell * ETswitch * ((cAMPlow - cAMP) / taudir)


#parameters
ETswitch = 0.0
IP3 = 0.0
alpha = 4.5e-06
cAMPlow = 0.2
cm = 5300.0
dact = 0.35
dinh = 0.4
dip3 = 0.5
f = 0.01
fer = 0.01
gca = 2000.0
girk = 1000.0
gk = 3500.0
kc = 0.15
ki = 0.5
kserca = 0.4
lambda = 1.25
perl = 0.0005
sh = 70.0
sigmav = 10.0
sm = 12.0
sn = 5.0
taudir = 20000.0
tauh = 20.0
taun = 20.0
vca = 25.0
vh = -20.0
vk = -75.0
vm = -20.0
vn = -16.0
cell = 1.0


#initial values
c = 0.3
cAMP = 1.0
cer = 260.0
h = 0.0
inh = 1.0
V = -60.0
n = 0.0


#assignment rules
!F igirk = girk * h * (V - vk)
!F ica = gca * minf * (V - vca)
!F jerp = kserca * c
!F jerip3 = o * (cer - c)
!F hinfer = 1.0 / (1.0 + c / dinh)
!F ainf = 1.0 / (1.0 + dact / c)
!F ik = gk * n * (V - vk)
!F perl_inf = inh * cAMP * pow(c, 4.0) / (pow(ki, 4.0) + pow(c, 4.0))
!F jertot = jerleak + jerip3 - kserca * c
!F jerleak = perl * (cer - c)
!F jmemtot = -(alpha * gca * 1.0 / (1.0 + exp((vm - V) / sm)) * (V - vca) + kc * c)
!F o = pow(ainf, 3.0) * pow(binf, 3.0) * pow(hinfer, 3.0)
!F binf = IP3 / (IP3 + dip3)
!F hinf = 1.0 / (1.0 + exp((vh - V) / sh))
!F ninf = 1.0 / (1.0 + exp((vn - V) / sn))
!F minf = 1.0 / (1.0 + exp((vm - V) / sm))