In this work, we consider a model which describes the CEF and magnetic interaction applied to two sublattices, where the magnetic interaction is described by a Hamiltonian who considers the Zeeman and exchange interactions in the molecular field approximation. The isothermal magnetic entropy change ΔS as a function of temperature was computationally obtained in the [001] crystallographic direction, from the equation of state of magnetization M(H, T) and entropy S (H, T) by diagonalization of the total Hamiltonian (HTot). We obtained a good agreement between the theory and the experiment was using the estimated crystalfield parameters and the exchange integral for Nd2Ni2In.