function b=betah(v)
%filename: betah.m
theta = (v+27)/5;
b=4.0/(1+exp(-theta));
