/*BiDiag solve the bi-diagonal matrix for ionization balance */
void BiDiag(long int nelem, 
  int lgPrintIt);