function r= alpha_m(v)
r =  .1*(v+40)/(1-exp(-(v+40)/10));
