Visto che i dipoli magnetici sono indipendenti, la loro energia dipenderà solo dal loro orientamento rispetto al campo: supponiamo ce ne siano

controallineati e

allineati; l'energia del sistema sarà
\mu B)
, con

.
Per la distribuzione di Boltzmann, la densità di probabilità che il sistema abbia energia

è
=Ae^{-\frac{U}{kT}})
, dove

è una costante. Per determinarla, la somma sui valori discreti di

deve essere 1:
=1=A\sum e^{-\frac{U}{kT}}=A\cdot exp(\frac{N\mu B}{kT})\sum_{m=0}^{N}(e^{-\frac{2\mu B}{kT}})^m \equiv Ae^{N\alpha}\sum_{m=0}^{N}e^{-2\alpha m})
, dove

.
Conoscendo il risultato della serie geometrica, si ottiene
\alpha}})
.
L'energia media

sarà:
U(m)=A\mu B\sum_{m=0}^{N}(2m-N)e^{(N-2m)\alpha}=-N\mu B+2\mu B \frac{1-e^{-2\alpha}}{1-e^{-2(N+1)\alpha}}\sum_{m=0}^{N}m(e^{-2\alpha})^m)
Per trovare
=\sum_{m=0}^{n}mx^m)
, riscriviamola come:
^2}(nx^{n+1}-(n+1)x^n+1))
Sfruttando questo risultato, ottengo
Per trovare

, si deve avere

;

molto elevata presumo significhi

, così da poter espandere

attorno a

; in questo modo otterrei
\approx U_0-\frac{2\mu^2B^2(N+1)}{kT}+O((\frac{\mu B}{kT})^2)))
Seguirebbe che
-U(T)}{\Delta T}=C= \frac{2\mu^2B^2(N+1)}{kT^2})