Ну это в общем случае все равно численный расчет, просто найдите подходящую программу и все. Из готового вы возможно сможете найти только какие-нибудь предпосчитанные частный случаи типа "шар в поле точечного диполя"
Сам принцип расчета прост, сила действующая на каждую точку ферромагнетика простенько вычисляется из ее намагниченности
и градиентов магнитного поля. Сама
тоже
вроде бы тривиально вычисляется просто умножением магнитной восприимчивости ферромагнетика на величину поля. Но вот тут то и случается вычислительная засада - вы умножаете внешнее магнитное поле на восприимчивость и вроде бы получаете
. Но теперь ферромагнетик с ненулевым
и сам создает магнитное поле, складывающееся с внешним. И уже эту сумму нужно умножать на восприимчивость чтобы получить новое
, новое поле, а значит опять новую
. И вот чтобы это посчитать без бесконечной рекурсии нужны уже нетривиальные методы. В каких то частных случаях может есть аналитическое решение, но в общем - нет