function rate = beta_n(v)
  rate = 0.125*exp(-v/80.0);
