/*rayleh compute Rayleigh scattering cross section for Lya */
double rayleh(double ener);