Модель Джилса-Атертона

[Serebro_2009]

Помогите разобраться. Хочу отобразить графически, результат моделирования используя модель Джилса-Атертона. Для простейшего случая, напряженность магнитного поля изменяется по синусоидальному закону:
$H=H_m\sin(2\pi Ft)$,
где $H_m$ - максимальная напряженность внешнего магнитного поля, $F$-частота

Далее, для реализации без гистерезисной намагниченности используется формула:
$M_{an}({H,M})=M_{s}(\coth({\frac{H+aM}{a}})-{\frac{a}{H+aM}})$

Однако, при поставление в нее значений функции описывающей напряженность внешнего поля (Н), получается...ломанный график. Подскажите, может что-то не то делаю?

Отделено от темы Моделирование частных петель гистерезиса магнитомягких мат.

[profrotter]

Тема перемещена из форума «Механика и Техника» в форум «Карантин»
Причина переноса:
Текст или формулы замещены рисунками.

(Подробно)

Forum Administration в сообщении #27356 писал(а):

I. Нарушения и санкции
1) Нарушением считается:
м) Набор любых формул без использования системы набора $\TeX$ (краткие инструкции по набору можно прочесть здесь, здесь или здесь). Использование картинок в качестве замены текста и формул (за исключением геометрических чертежей, сложных диаграмм и таблиц). Несоблюдение правил использования внешних ссылок (см. п. III-5).

Инструкции по оформлению формул здесь или здесь (или в этом видеоролике).
Оформите тему в соответствии с правилами форума и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума

[profrotter]

Тема перемещена из форума «Карантин» в форум «Механика и Техника»
Причина переноса: Вернул.

[Serebro_2009]

Не могу понять, как организуется сам процесс решения. Выражение для вычисления без гистерезисной намагниченности:
$M_{an}({H,M})=M_{s}(\coth({\frac{H_e+aM}{a}})-{\frac{a}{H_e+aM}})$
Выражение для определения эффективного поля Вейса:

$H_e=H+aM$
где $M$ - полная намагниченность, состоящая из двух компонент:
$M_{irr}$ - обратимая намагниченность
$M_{rev}$ - не обратимая намагниченность
Полная намагниченность, в соответствии с моделью Джилса-Атертона, может быть вычислена исходя из выражения:
${\frac{dM}{dH}={\frac{H_{an}-M}{K}\sigma}+\frac{C}{1+C}\frac{dM_{an}}{dH}$
А как программно реализовать последовательность этих действий? Т.е. для численного дифференцирования последнего выражения, необходимо знать значения: $M_{ank+1}$ и $M_{ank}$?