\[ y'(x)=\frac {x^2 y(x) \sqrt {x^2+y(x)^2}+x^3 \left (-\sqrt {x^2+y(x)^2}\right )+y(x)}{x} \] ✓ Mathematica : cpu = 0.121895 (sec), leaf count = 93
\[\left \{\left \{y(x)\to \frac {x \left (-2 e^{\frac {1}{3} \sqrt {2} \left (3 c_1+x^3\right )}+e^{\frac {2}{3} \sqrt {2} \left (3 c_1+x^3\right )}-1\right )}{2 e^{\frac {1}{3} \sqrt {2} \left (3 c_1+x^3\right )}+e^{\frac {2}{3} \sqrt {2} \left (3 c_1+x^3\right )}-1}\right \}\right \}\]
✓ Maple : cpu = 0.273 (sec), leaf count = 49
\[ \left \{ \ln \left ( 2\,{\frac {x \left ( \sqrt {2\, \left ( y \left ( x \right ) \right ) ^{2}+2\,{x}^{2}}+y \left ( x \right ) +x \right ) }{y \left ( x \right ) -x}} \right ) +{\frac {\sqrt {2}{x}^{3}}{3}}-\ln \left ( x \right ) -{\it \_C1}=0 \right \} \]