/*smeets generate hydrogenic collisional ionization rate coefficients */
double smeets(long int n);