cmdVdt=−im+IeA
τm(V)dmdt=m∞(V)−m
τh(V)dhdt=h∞(V)−h
τn(V)dndt=n∞(V)−n
τa(V)dadt=a∞(V)−a
τb(V)dbdt=b∞(V)−b
im=g¯_L(V−EL)+g¯_Nam3h(V−E_Na)+g¯_Kn4(V−E_K)+g¯_Aa3b(V−EA)
αn(V)=0.02mV−1(V+45.7mV)1−exp(−0.1mV−1(V+45.7mV))
αm(V)=0.38mV−1(V+29.7mV)1−exp(−0.1mV−1(V+29.7mV))
αh(V)=0.266exp(−0.05mV−1(V+48.0mV))
βn(V)=0.25exp(−0.0125mV−1(V+55.7mV))
βm(V)=15.2exp(−0.0556mV−1(V+54.7mV))
βh(V)=3.81+exp(−0.1mV−1(V+18mV))
τn(V)=1.0msαn(V)+βn(V)
τm(V)=1.0msαm(V)+βm(V)
τh(V)=1.0msαh(V)+βh(V)
τa(V)=0.3632ms+1.158ms1.0+exp(0.0497mV−1(V+55.96mV))
τb(V)=1.24ms+2.678ms1.0+exp(0.0624mV−1(V+50.0mV))
a∞(V)=(0.0761*exp(0.0314mV−1(V+94.22mV))1+exp(0.0346mV−1(V+1.17mV))))1/3ms
b∞(V)=(11+exp(0.0688mV−1(V+53.3mV)))4ms
n∞(V)=αn(V)τn(V)
m∞(V)=αm(V)τm(V)
h∞(V)=αh(V)τh(V)
m(t=0)=0.010
n(t=0)=0.156
h(t=0)=0.966
a(t=0)=0.540
b(t=0)=0.289
V(t=0)=−68.0mV
cm=0.1μFmm2
IeA=0.35μAmm2
g¯L=0.003mSmm2
g¯Na=1.2mSmm2
g¯K=0.2mSmm2
g¯A=0.477mSmm2
EL=−17.0mV
ENa=55.0mV
EK=−72.0mV
EA=−75.0mV
M∞=11+exp(−(V+57mV)/6.2mV)
H∞=11+exp(−(V+81mV)/4mV)
τM=0.612ms+1msexp(−(V+132mV)/16.7mV)+exp((V+16.8mV)/18.2mV)
if V<−80mV:
τH=1msexp((V+467mV)/66.6mV)
else:
τH=28ms+1msexp(−(V+22mV)/10.5mV)
c∞=([Ca2+][Ca2+]+3μM)11+exp(−(V+28.3mV)/12.6mV)
τC=90.3ms−75.1ms1+exp((−V+46mV)/22.7mV)