Monte Carlo approach to phase transitions in ferroelectromagnets