Coupled hydrogen--helium Saha equations (source code)

= Coupled hydrogen--helium Saha equations

Writing $\mathcal F=n_e/n_B$ couples the two equilibrium equations:
$$
\mathcal F\frac{X_{H^+}}{1-X_{H^+}}=K(T)e^{-I_H/(k_BT)},\qquad
\mathcal F\frac{X_{He^+}}{1-X_{He^+}}=2K(T)e^{-I_{He}/(k_BT)}.
$$