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