Solution
= Solution
In natural units, the equilibrium <phase-space distribution function> is
$$
f(\mathbf p)=\frac{1}{\exp[(E(\mathbf p)-\mu)/T]\mp1},
\qquad E(\mathbf p)=\sqrt{\mathbf p^2+m^2},
$$
where the minus sign gives the <Bose-Einstein distribution> for <bosons> and the plus sign gives the <Fermi-Dirac distribution> for <fermions>.