Per determinare la capacità termica di un solido usando un metodo basato sulle frequenze di vibrazione del cristallo (vedi) si deve calcolare ∫ _{[0, k]} x⁴·e^x/(e^x−1)² dx.  Calcola in modo approssimato (usando del software) il valore di questo integrale quando k = 1.