function b=betan(v)
%filename: betan.m
theta = (v+57)/40;
b=0.5*exp(-theta);
