Ve -e_const*V "Electron potential energy" eta beta*(Ef-Ve) "Energy factor" n_TF N0TF*schr.FD_half(eta) "Electron density, Thomas-Fermi approx" rho_TF -e_const*n_TF "Space charge density, Thomas-Fermi approx"