\[ \boxed { xy \left ( x \right ) \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) ^{2}+ \left ( \left ( y \left ( x \right ) \right ) ^{2}+{x}^{2} \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +xy \left ( x \right ) =0} \]
Mathematica: cpu = 0.008501 (sec), leaf count = 49 \[ \left \{\left \{y(x)\to \frac {c_1}{x}\right \},\left \{y(x)\to -\sqrt {2 c_1-x^2}\right \},\left \{y(x)\to \sqrt {2 c_1-x^2}\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 35 \[ \left \{ y \left ( x \right ) =\sqrt {-{x}^{2}+{\it \_C1}},y \left ( x \right ) ={\frac {{\it \_C1}}{x}},y \left ( x \right ) =-\sqrt {-{x}^{2} +{\it \_C1}} \right \} \]