function a=alphah(v)
%filename: alphah.m
theta = (v+50)/18;
a=0.128*exp(-theta);
