U=-\sum_i \vec{m}_i\cdot\vec{B}=-\vec{M}\cdot\vec{B}\sim -\frac{1}{r}