Описания формулы я не нашел, пришлось взять что попроще.

В статье по ссылке исчерпывающее описание формулы.

именно так для теста и сделал

Это было ошибкой.

ак мне нужен вначале нужен график, который слева.

Я не знаю, что такое "график слева". Есть восходящая и нисходящая ветви. У Вас же часть от восходящей, а часть от нисходящей. В итоге - сумбур вместо музыки. ОБЕ ветви определены на одном и том же интервале значений напряженности магнитного поля, и вычислять значение намагниченности для обоих ветвей надо для каждой точки внутри -H_max...H_max. Принцип подходит прекрасно. Или мне надо обязательно запостить скришноты генерации графика в R или Origin? Проблему при использовании данной модельной функции я уже указал.

Но один скриншот ничего не даст, как и таблица. Устройство то для того и сделано чтобы снимать реальные данные.

Да, БЕЗУСЛОВНО, от цикла к циклу перемагничивания одни параметры кривой могут меняться. А другие - нет. Или меняться крайне незначительно по сравнению с первыми. Но для начала было бы неплохо выяснить, можно ли в принципе описать данные тем или иным образом. Хотя бы на одном примере.