The single Ion Anisotropy (SIA) term is now included in the main Hamiltonian. Also, a feature-like calculation in terms of the parameter (J) has been added where the temperature scale will be considered as KbT/J.