pub fn saha_hydrogen(z: f64, cosmo: &Cosmology) -> f64
Hydrogen Saha ionization fraction.
Solves X_e²N_H / (1−X_e) = (m_e k_B T / 2πℏ²)^{3/2} exp(−E_Rydberg/kT) for X_e.