PX156 - A3b - planck's function

En=nhf

where, n∈Z, h= planck's constant

⟨E⟩=hfexp⁡(hfkBT)−1