# Excitatory cortical neurons - Golomb and Amitai
#
#
# minf^3 h --- sodium
# pinf -- persistent sodium
# n^4 -- delayed rect
# ainf^3 b - a-current
# z - M-current
number Cm=1.0
number pms=3,pns=4
number Vglu=0.0
number thetaa=-50.0,sigmaa=20.0,thetab=-80.0,sigmab=-6.0,tauBs=15.0
# 
par i0=0,ip=2,ton=50,toff=300
p gNa=24.0,gKdr=3.0,gL=0.02
p gA=1.4,gNaP=0.07,gZ=1.0
p thetaz=-39.0,sigmaz=5.0,tauZs=75.0
p phi=2.7
p thetam=-30.0,sigmam=9.5,thetah=-53.0,sigmah=-7.0
p thetan=-30.0,sigman=10.0,thetap=-40.0,sigmap=5.0
p VNa=55.0,VK=-90.0,VL=-70.0,t_tauh=-40.5,t_taun=-27.0
#
GAMMAF(V,theta,sigma)=1.0/(1.0+exp(-(V-theta)/sigma))
#
i(t)=i0+ip*heav(t-ton)*heav(toff-t)
#
V'=(-gL*(V-VL)-gNa*(Minf^pms)*h*(V-VNa)\
     -gNaP*Pinf*(V-VNa)-gKdr*(n^pns)*(V-VK)-\
      gA*Ainf^3*b*(V-VK)-gZ*z*(V-VK)+i(t))
h'=phi*(GAMMAF(V,thetah,sigmah)-h)/(1.0+7.5*GAMMAF(V,t_tauh,-6.0))
n'=phi*(GAMMAF(V,thetan,sigman)-n)/(1.0+5.0*GAMMAF(V,t_taun,-15.0))
b'=(GAMMAF(V,thetab,sigmab)-b)/tauBs
z'=(GAMMAF(V,thetaz,sigmaz)-z)/tauZs
#
Minf=GAMMAF(V,thetam,sigmam)
Pinf=GAMMAF(V,thetap,sigmap)
Ainf=GAMMAF(V,thetaa,sigmaa)
#
#
V(0)=-73.8
h(0)=0.95
n(0)=0.0122
b(0)=0.264
z(0)=0.00093
#
@ total=1000,dt=.5,meth=qualrk,tol=1e-5,xhi=1000,ylo=-85,yhi=30
@ bounds=1000
done