function n_e_inf=n_e_inf(v);
alpha_n=0.032*(v+52)./(1-exp(-(v+52)/5));
beta_n=0.5*exp(-(v+57)/40);
n_e_inf=alpha_n./(alpha_n+beta_n);