function r=beta_m(v)
r =  4*exp(-(v+65)/18);
