/*hpfit state specific photoionization cross sections for model hydrogen atom */
double hpfit(long int iz, 
	  long int n, 
	  double e );