Интересует вывод
Видимо, прав
peregoudov. Формула выводится так. Сила, действующая на точечный диполь, равна (в компонентах векторов в декартовых координатах)


-

-я компонента магнитного момента, по повторяющимся индексам предполагается суммирование. Сила, действовавшая на кусочек среды, который заместила частица, равна

Если считать, что

слабо меняется на размерах частицы, то получится, что сила, действующая на кусочек жидкости равна

Поскольку жидкость в равновесии, со стороны среды на этот кусочек действует такая же сила с обратным знаком.
Если теперь вместо жидкости вставить частицу, то надо выполнить те же расчеты, только поле надо считать с учетом формы и магнитной проницаемости частицы. Однако, для "нормальных" диа- и парамагнетиков

очень мало отличается от единицы. Поэтому для них можно с хорошей точностью написать

если градиент поля мало меняется на размерах частицы.