Описано чіткий алгоритм формування системи нелінійних алгебраїчних рівнянь виходячи з мінімуму функціонала для розрахунку розподілу температури в тривимірній області з безгістерезисними нелінійними анізотропними середовищами. Сформульовано краєву задачу розрахунку тривимірного стаціонарного потенціального температурного поля. Побудовано скінченно-елементну модель розрахунку розподілу температури всередині тривимірної області, заповненої безгістерезисними нелінійними анізотропними середовищами.