\[ \boxed { \left ( 2\, \left ( y \left ( x \right ) \right ) ^{2}{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +{x}^{2} \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) +2\,y \left ( x \right ) \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) ^{3}+3\,x{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +y \left ( x \right ) =0} \]
Mathematica: cpu = 42.477894 (sec), leaf count = 45 \[ \text {DSolve}\left [\left (x^2+2 y(x)^2 y'(x)\right ) y''(x)+2 y(x) y'(x)^3+3 x y'(x)+y(x)=0,y(x),x\right ] \]
Maple: cpu = 1.529 (sec), leaf count = 54 \[ \left \{ y \left ( x \right ) ={\it ODESolStruc} \left ( {\it \_b} \left ( {\it \_a} \right ) ,[ \left \{ \left ( {\it \_b} \left ( {\it \_a } \right ) \right ) ^{2} \left ( {\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) \right ) ^{2}+{{\it \_a}}^{2}{\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) +{\it \_a }\,{\it \_b} \left ( {\it \_a} \right ) +{\it \_C1}=0 \right \} , \left \{ {\it \_a}=x,{\it \_b} \left ( {\it \_a} \right ) =y \left ( x \right ) \right \} , \left \{ x={\it \_a},y \left ( x \right ) ={\it \_b} \left ( {\it \_a} \right ) \right \} ] \right ) \right \} \]