y'(x)=\frac{x^2 \left(-\sqrt{x^2+y(x)^2}\right)+x y(x) \sqrt{x^2+y(x)^2}+x y(x)+y(x)}{x (x+1)}