$${\rho_H = \frac{1}{2} \mu \mu_0 H^2}$$