The (anti)ferromagnetic and ferroelectric transitions in some multiferroic compounds seem to be strongly correlated. Even for systems that do not show spontaneous ferroelectricity such as the LiMPO 4 (M=Mn, Fe, Co, Ni) compounds, the coupling between magnetic and electric degrees of freedom is evident experimentally. Here, we present a simple numerical calculation to simulate this coupling that leads to the two transitions. We assume a magnetic sublattice consisting of classical magnetic moments coupled to a separated nonmagnetic sublattice consisting of classical electric dipoles. The coupling between them is realized through a phenomenological spin-lattice Hamiltonian, and the solution is obtained using the Monte Carlo technique. In the simplest version, the magnetic system is 2D Ising (anti)ferromagnetic lattice, with nearest neighbors interactions only, and the electric moments are permanent moments, coupled electrically. Within this approximation, the second order magnetic transition induces ferroelectricity in the electric dipoles. We show that these calculations can be extended to other magnetic systems, (x-y model and 3D Heisenberg) and to systems where the electric moments are created by strains, generated via spin-lattice coupling, so the model can be applied to model realistic systems such as the olivines mentioned above. © Published under licence by IOP Publishing Ltd.2731 Fiebig, M., (2005) J.Phys. D: Appl. Phys., 38 (8), p. 123Rivera, J.-P., (1996) Ferroelectrics, 161, p. 147Vaknin, D., Zaretsky, J.L., Rivera, J.-P., Schmid, H., (2004) Phys Rev. Lett., 92 (20), p. 207201Kosterlitz, J.M., Nelson, D.R., Fisher, M.E., (1976) Phys Rev., 13 (1), p. 412