Частные петли гистерезиса

[salevan]

Пытаюсь сделать модель ТТ в программной среде Matlab, с помощью модели Дж. Чана.
Предельная петля и симметричные частные петли нормально сходятся с экспериментом, а вот с несимметричными - беда.
Принцип построения несимметричной петли я более менее понял, но я 4 года к матлабу не подходил вообще, и возникла проблема в программировании этих сдвигов.
Можете помочь пожалуйста с самим кодом или ещё раз попонятнее объяснить метод построения несимметричных петлей.

Заранее спасибо!