\[ y'(x)=\frac {x^2 \sqrt {x^2+y(x)^2}+y(x)}{x} \] ✓ Mathematica : cpu = 0.123257 (sec), leaf count = 20
\[\left \{\left \{y(x)\to x \sinh \left (\frac {1}{2} \left (x^2+2 c_1\right )\right )\right \}\right \}\] ✓ Maple : cpu = 3.691 (sec), leaf count = 30
\[\ln \left (\sqrt {y \relax (x )^{2}+x^{2}}+y \relax (x )\right )-\frac {x^{2}}{2}-\ln \relax (x )-c_{1} = 0\]