function b=betan(v)
%filename: betan.m
theta = (v+70)/80;
b=0.125*exp(-theta);