Рассмотрена двумерная спин-псевдоспиновая система, которая моделирует плоскость CuO2 в купратах и учитывает спиновый обмен и зарядовые корреляции на узле и между узлами. Эта модель позволяет изучить конкуренцию зарядового и магнитного упорядочения в купратах. В работе мы исследовали термодинамическую устойчивость АФМ фазы для случаев сильного и слабого обмена. Для этого мы получили температурные зависимости химического потенциала и сжимаемости классическим методом Монте-Карло и сравнили результаты численных расчетов с аналитическими результатами в рамках среднего поля.