В рамках теории фаз Нееля проведено моделирование кривых намагничивания кубических магнетиков. Определены основные особенности хода кривых в различных кристаллографических направлениях. Получено хорошее согласие модельных кривых с экспериментальными данными.