Докажем это используя дифференциальные формы.
Надо доказать что 1-форма

.
Рассмотрим ее дифференциал - форму

.
Так-как степень характеристического многочлена равна трем (нечетна), то в каждой точке трехмерного пространства матрица этой 2-формы имеет действительное собственное значение

.
Это действительное значение обязательно равно нулю, ведь матрица кососимметрическая:

.
То есть в каждой точке есть вектор

, такой что

.
Вдоль этого поля можно провести линии (линии ротора формы

).
Возьмем какую-нибудь поверхность которая не касается этих линий, определим на ней произвольную функцию

и продолжим ее вдоль линий на окружающее пространство,

-

не меняется вдоль поля

.
Эта функция будет искомой.
В самом деле

(это объем на скаляр, значит можно подставлять

, но он обращается в нуль одновременно для

и

).
В координатах


.
Тогда

. Значит

, и функция

.
Тогда

.
-- Пт ноя 26, 2010 17:39:31 --Но все это локально. А глобально так делать нельзя.
Можно, когда линии ротора устроены просто, а если они перепутаны - то нет.
-- Пт ноя 26, 2010 17:51:33 --Замечу еще, что поле

.