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