### Izhikevich's Leak + Instantaneous Na_p (Persistent Sodium) 1D model. ### See Dynamical Systems in Neuroscience, section 3.1.2, pp 55-56 ### Parameters # Capacitance, microFarads params c=10 # Applied current, picoAmperes params iapp=0 # Conductances, milliSiemens params gl=19,gna=74 # Reversal potentials, milliVolts params el=-67,ena=60 # Sodium activation parameters, milliVolts params vhalf_m=1.5,k_m=16 ### Functions # Sodium activation function minf(v)=1/(1+exp((vhalf_m-v)/k_m)) ### Dynamics time' = 1 v' = (iapp-gl*(v-el)-gna*minf(v)*(v-ena))/c ### Initial Conditions init time=0,v=-60 ### Settings @ TOTAL=30,DT=.05,xlo=0,xhi=10,ylo=-70,yhi=60 @ xplot=time,yplot=v #set vvst {xplot=t,yplot=v,xlo=0,xhi=30,ylo=-90,yhi=40,total=100 \ # dt=.5,meth=qualrk,iapp=.1} done