Hi, everyone!
I am new here and I hope I've found the right place to ask my question.
I've been trying to simulate the magnetization process of a ferrimagnet using Mathcad. The theory behind this is quite simple - one writes down an energy expansion, there're three terms that signify how two magnetic moments interact with each and with magnetic field. Mathematically, one needs to minimize this energy with respect to the angles the magnetic moments make with the field and some symmetry direction of the crystal lattice. The minimization should be done for each magnetic field value, e.g., B = 0,1,2,...,20 T. Then, the total magnetization, Mtotal, can be calculated using the obtaines values of the angles.
I've prepared a Mathcad document for the simulation (see attachment). The program works but gives me an incorrect unphysical result. Could you take a look and possibly exlpain to me what I've done wrong?
Thank you in advance.
Regards,
Den