Ну пожалуйста, перешли в координаты , в которых векторное поле имеет вид

. Тут все группы симметрий выписываются сразу и тривиально. Думаете этим Ли занимался?
Не надо путать группу симметрий системы ОДУ первого порядка, и группу симметрий одного ОДУ второго или более высокого порядка. Последним соответствуют не абы какие преобразования в продолженном пространстве

, а только такие, которые получаются из симметрий пространства

по формулам продолжения, требующим сохранения равенств

,

и т.д. Т.е. преобразования производных вытекают из преобразования переменных.
Любые два невырожденных в точке векторных поля эквивалентны в окрестности этой точки. А здесь речь идет о глобально заданных векторных полях. Глобальной эквивалентности теорема о выпрямлении не даёт, так что не отвечает и на вопрос о том, что называют симметриями дифференциального уравнения.
Это не в ту степь. Речь идет именно о локальных симметриях.