“…Because of a limited number of exact methods available, approximate and numerical studies have their own importance. Referring to recent studies on layered magnetism in the literature, it is reported, by MFA and EFA in [23] and by MC simulations with Single cluster Wolff Algorithm in [24] that under certain range of different types of interaction strengths between lattice sites, different temperature dependencies of sublattice magnetisations cause the compensation point to appear in an Ising trilayer on square lattice. A quasi three-dimensional, spin-1/2, Ising trilayer stacking on square lattice with quenched non-magnetic impurity is studied in [44] and it is observed, the critical and compensation points, can be modulated to drift towards lower temperatures with increasing concentration of non-magnetic impurities and in [45], it was hinted that reduced residual magnetisation, another magnetocaloric quantity, possibly varies in a systematic way.…”