function b=betam(v)
%filename: betam.m
theta = (v+70)/18;
b=4.0*exp(-theta);
