Ich würde mal folgende Substitution versuchen:
u(x) = x*v(x) bzw. v = u/x
Daraus folgt: v' = u + x*u'
Teilt man die gegebene Differentialgleichung einmal durch x, so folgt:
y' = y/x + sqrt(x²+y²)
Mit der Substitution:
u + x*u' = u + sqrt(x²+x²u²) | -u, :x
u' = sqrt(1+u²)
Das lässt sich jetzt sehr einfach mit Separation der Variablen lösen.
∫1+u2du=∫dxu=sinh−1(u)=x−x0u=sinh(x−x0)y=x∗sinh(x−x0)
Das Integral über 1/sqrt(1+u²) lässt sich leicht durch eine Substitution mit u=sinh(v) und ausnutzen des hyperbolischen Pythagoras 1=cosh²(x) - sinh²(x) lösen.