> Подозреваю, что решение тут гораздо проще: нужна гибкая модельная функция с ясно задаваемыми параметрами

ну так правильно, я именно это и сделал, только разжевал в уже готовую функцию, которую просто надо воткнуть в шарп. Я с шарпами не в ладах, то есть не умею от слова никак, посему оставил формулы на съедение ТС. А так в моих формулах явно задаются точки Mr, Mc, HC (мне б телепатом подрабатывать?) и есть готовая функция, (одна из трех, в зависимости от диапазона аргумента) то осталось ее вбить в шарп, скомпилить и радоваться. В точках MR и HC имеется гладкость до второго порядка, а на краю в одном случае нулевая производная, а в другом случае - нулевая вторая производная. Собственно, как и должно быть во всех кривых намагничивания.

Если немного пооптимизировать, то можно в формулах заметить одинаковые члены, которые правильно вычислить заранее, тогда формулы станут еще короче. Надеюсь, что ТС это сам осилит, если что, готов подсказать как проще.