function rate = alpha_h(v)
  rate = 0.07*exp(-v/20.0);
