function a=alphah(v)
%filename: alphah.m
theta = (v+70)/20;
a=0.07*exp(-theta);