In [1]:

    

from scipy import *
from model.parameters import Mitral
from brian import *
psmt = Mitral()


    

In [2]:

    

tau_m = linspace(psmt.C_m / psmt.g_L, psmt.C_m / psmt.g_L * 10, 10)
R_m   = 1./psmt.g_L
E_L   = psmt.E_L
V_r   = psmt.V_r
V_t   = psmt.V_t
tau_ref = psmt.t_refract


    

In [3]:

    

currents = array(linspace(1e-10, 0.04, 10000))


    

In [4]:

    

def f_isi(I_e, tau_m):
    mask = (R_m*I_e + E_L - V_t)>0.
    tmp=zeros(I_e.shape)
    return r_[tmp[~mask],1./(tau_ref+tau_m*log((R_m*I_e[mask] + E_L - V_r)/(R_m*I_e[mask] + E_L - V_t)))]


    

In [5]:

    

for i in xrange(10):
    plot(currents, f_isi(currents, tau_m[i]), label="tau_m = "+str(tau_m[i]))
    legend()
show()


    

In [5]: