###########################################################
#            Chay Cook 1988 b-cell model (FAST SUBSYSTEM)
###########################################################
#
#
###########################################################
#            Problem definition
###########################################################
#
#       parameter values defining model
#
p gib=250,gsb=10,gkb=1300,glb=50,vca=100,vk=-80,vl=-60
p cm=4524,tnb=9.09,vm=-22,vn=-9,vs=-22,sm=7.5,sn=10
p alph=0.000005727
#
#      parameters (Type 1a S-wave burst)
#
p f=0.002,kc=0.027,ss=10,lam=0.95,tsb=0.1
#
#      parameters (Type 1b )
#
#p f=0.00005,kc=0.027,ss=10,lam=0.17,tsb=0.1
#
#      parameters (Type III)
#
#p f=0.002,kc=0.022,ss=10,lam=0.1,tsb=0.1
#
#   Bifurcation parameter
#
p c=1
#
#       function definitions
#
a(v,c)=(vs+ss*log(c)-v)/(2*ss)
ts(v,c)=tsb/(2*cosh(a(v,c)))
tn(v)=tnb/(1+exp((v-vn)/sn))
si(v,c)=1/(1+exp(2*a(v,c)))
mi(v)=1/(1+exp((vm-v)/sm))
ni(v)=1/(1+exp((vn-v)/sn))
ii(v)=gib*mi(v)*(v-vca)
is(v,s)=gsb*s*(v-vca)
ik(v,n)=gkb*n*(v-vk)
il(v)=glb*(v-vl)
ica(v,s)=ii(v)+is(v,s)
#
#      differential equations
#
v'=-(ica(v,si(v,c))+ik(v,n)+il(v))/cm
n'=lam*(ni(v)-n)/tn(v)
#
#      initial conditions (Type 1a S-wave burst)
#
v(0)=-56.482
n(0)=0.0085927
#
#      Numerics and setup (Type 1a S-wave burst)
#
@ total=100,xhi=100,ylo=-60,yhi=-10,maxstor=1000,dt=1
#
d