/*hrfit calculates state-specific recombination rates for H and H-like ions*/
void hrfit(long int iz, 
	  long int n, 
	  double t, 
	  float *r);