Вернусь к результатам первоначальной задачи.
Уравнения движения не представляют труда для интегрирования, однако уже здесь проявляются некоторые особенности негологомных систем.
Имеем два перых интеграла, линейных по скоростям.

(рассматривать пока не будем).
Вот другой:

.
На его примере покажем, как из общей конструкции получается линейный интеграл.
Введем три линейно независимые 1-формы

напомню,что

- связь и три дуальных им векторных поля

. Легко видеть, что

Коммутаторы
![$[X_1,X_2]=X_3,[X_1,X_3]=0,[X_2,X_3]=0$ $[X_1,X_2]=X_3,[X_1,X_3]=0,[X_2,X_3]=0$](https://dxdy-03.korotkov.co.uk/f/a/7/e/a7eeb92636a838f7f069a42b8bc0ad8b82.png)
. Здесь обнаруживается, что мы имеем дело с классическими объектами - контактной структурой

на

и нильпотентной алгеброй Ли, натянутой на

- алгеброй Гейзенберга.
Рассмотрим квадратичную форму, соответствующую кинетической энергии

и векторное поле

. Пусть

- производная Ли. Тогда в нашем случае

, где

- 1-форма.
Последнее равенство влечет существование линейного по скоростям интеграла. Поскольку

, то первый интеграл имеет вид

, что и есть искомый первый интеграл. Описанная выше конструкция есть аналог теоремы Нётер в голономном случае.