This is an integrate-and-fire model to simulate an action potential whose exponential term mimics the increase in the sodium current due to activation. The time evolution of the voltage can be described by the following differential equation:
${\displaystyle C_m \frac{dV}{dt} = - g_L(V-E_L) + g_L\Delta T \exp{\frac{(V-V_{thr})}{\Delta T}} + I(t)},$
where $C_m$ is the membrane capacitance, $V$ is the membrane potential, $g_L$ the leak conductance and $E_L$ the resting potential. $\Delta T$ is a parameter that indicates the slope of the sodium activation and $V_{thr}$ is the threshold potential
In [8]:
def MembranePotential(v, current):
"""
Cm * dV/dt = -gl*(v-Vrest) + I
"""
gl = 0.16 # in mV
Cm = 4.9 #
Vrest = -70. # in mV
I = current
return (gl*(v-Vrest) + I)/Cm
In [7]:
MembranePotential(-71,0)
Out[7]:
In [ ]: