function a=alpham(v)
%filename: alpham.m
theta=(v+45)/10;
if(theta==0)   %check for case that gives 0/0
  a=1.0;  %in that case use L'Hospital's rule
else
  a=1.0*theta/(1-exp(-theta));
end