/*fhummr evaluate Hummer's betaF(beta) function */
double fhummr(double beta); /* beta is ratio of continuum to mean line opacity,
									 * returns dest prob = beta F(beta) */