\[ \boxed { x{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( y \left ( x \right ) \right ) ^{2}+{x}^{2}=0} \]
Mathematica: cpu = 0.014502 (sec), leaf count = 32 \[ \left \{\left \{y(x)\to \frac {x \left (-c_1 J_1(x)-Y_1(x)\right )}{c_1 J_0(x)+Y_0(x)}\right \}\right \} \]
Maple: cpu = 0.062 (sec), leaf count = 40 \[ \left \{ y \left ( x \right ) =-{\frac {{\it \_C1}\,x{{\sl Y}_{1}\left (x \right )}}{{\it \_C1}\,{{\sl Y}_{0}\left (x\right )}+{{\sl J}_{0}\left (x \right )}}}-{\frac {{{\sl J}_{1}\left (x\right )}x}{{\it \_C1}\,{{\sl Y} _{0}\left (x\right )}+{{\sl J}_{0}\left (x\right )}}} \right \} \]