Для интереса решил найти первые интегралы системы (см. оффтоп).
(Оффтоп)
Перейдём к переменным

по формулам

В новых переменных система примет вид:

Решение:

Индекс

обозначает начальное значение переменной (при

).
Трактуя

как декартовы координаты, перейдём стандартным образом к цилиндрическим

:

Отсюда следует

, что вместе с

даёт решение в цилиндрических координатах.
Перепишем решение в виде

Будем считать, что в правых частях множители перед

не равны нулю (случай равенства нулю надо рассмотреть отдельно). Тогда отношение любых двух левых частей — константа, так что получаем три первых интеграла. Конечно, они не являются функционально независимыми. Им соответствуют три семейства поверхностей (

— произвольная ненулевая константа):

— конусы

— «логарифмически-спиральные цилиндры»

— «логарифмические винтовые поверхности»
Совершая обратные преобразования координат, можно выразить первые интегралы через

.