function r=bh(v)
r =  1/(1+exp(-(v+35)/10));
