Дифференциальные уравнения с Mathematica

В Mathematica необходимо поместить его в следующую форму:

eqn = x*y'[x] - y[x] == Sqrt[1 + x^2]*y'[x] + Sqrt[1 + y[x]^2]
DSolve[eqn, y[x], x]

С этим я добираюсь:

{{y [x]-> InverseFunction [1/2 (ArcSinh[#1] + № 1 (-#1 + Sqrt [1 + # 1^2])) и] [1/2 (-x (x + Sqrt [1 + x^2]) - ArcSinh[x]) + C[1]] }}

с соблюдающим предупреждением:

Решите:: tdep: уравнения, кажется, включают переменные, которые будут решены для чрезвычайно неалгебраическим способом.

От справки кажется, что решатель встретился с трансцендентной функцией. Это о том, насколько я могу взять это. Надежда это помогает получить Вас движение.