DG representations of the equilibrium Fermi–Dirac distribution with the bound-enforcing limiter applied for electron neutrinos in conditions from a collapsed stellar core around 100 ms after core bounce. The left panel shows results for T = 15 MeV and ﹩{\mu }_{{\nu }_{{{\rm{e}}}}}=113﹩ MeV, characteristic of conditions for electron neutrinos at high density (ρ ∼ 10¹⁴ g cm⁻³ and Y_e ∼ 0.27). The middle panel shows results for T = 7.6 MeV and ﹩{\mu }_{{\nu }_{{{\rm{e}}}}}=3.5﹩ MeV, characteristic of conditions at lower density (ρ ∼ 10¹² g cm⁻³ and Y_e ∼ 0.13). The right panel shows the rms error of the approximation against the exact Fermi–Dirac distribution for the low-density case in the middle panel vs. the number of energy elements ﹩{\mathfrak{N}}﹩ for various polynomial degrees k, using the energy grids specified in Table 2.