Function dispatched to the following (C++) functions:
[1] (obe: OneBodyElementsOnGrid, mu: float, beta: float) -> float
Compute number of particles n = ∑f(β(e(k)-μ))