Monte Carlo Metropolis simulation of interacting anisotropic polarizable spins on a lattice