Не могу понять, как организуется сам процесс решения. Выражение для вычисления без гистерезисной намагниченности:

Выражение для определения эффективного поля Вейса:

где

- полная намагниченность, состоящая из двух компонент:

- обратимая намагниченность

- не обратимая намагниченность
Полная намагниченность, в соответствии с моделью Джилса-Атертона, может быть вычислена исходя из выражения:

А как программно реализовать последовательность этих действий? Т.е. для численного дифференцирования последнего выражения, необходимо знать значения:

и

?