To evaluate this integral, the function N ¯ ¯ ( E ) , the number of phonons with energy E , must also be known. Phonons obey Bose–Einstein statistics, and their distribution is given by the Bose–Einstein statistics formula: