Приведены результаты численного моделирования в рамках модифицированного метода Монте-Карло с алгоритмом термостата для псевдоспиновой модели редкоземельных ортоникелатов. Построены температурные фазовые диаграммы для различных степеней заполнения и при различных параметрах модели, исследовано влияние локальных корреляций на критические температуры модельного никелата. Показана возможность детектирования фазово неоднородных состояний. Результаты численного моделирования показывают хорошее качественное согласие с аналитическими результатами в методе среднего поля. Ключевые слова: редкоземельные ортоникелаты, метод Монте-Карло, приближение среднего поля.