I want to simulate hysteresis loops using an expression similar to this one of the figure, I can already simulate using coherent rotation of the magnetization, however I have never worked with this type of Hamiltonian model ... can anyone help me write or give a hint how to implement it in Maple, Matlab, or other simulation program.
The reference of the figure: https://doi.org/10.1088/0953-8984/19/23/236225
Thanks a lot.
Oreci.