* implement simple minded Y_e(betaeq) code to temporarily
  replace Jeff's code (which is more beautiful, but produces
  strange jumps in Y_e).