Original: n_i = frac {g_i} {e ^ {(varepsilon_i - mu) \/O/ kT} -1}
Simplified: n (varepsilon) = frac {1} {e ^ {(varepsilon - mu) \/O/ kT} -1}