\[ \boxed { \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) -{x}^{2}y \left ( x \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) -x \left ( y \left ( x \right ) \right ) ^{2}=0} \]
Mathematica: cpu = 63.610077 (sec), leaf count = 34 \[ \text {DSolve}\left [x^2 (-y(x)) y'(x)+y'(x) y''(x)-x y(x)^2=0,y(x),x\right ] \]
Maple: cpu = 0.874 (sec), leaf count = 46 \[ \left \{ y \left ( x \right ) ={\it ODESolStruc} \left ( {\it \_b} \left ( {\it \_a} \right ) ,[ \left \{ - \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{2}{{\it \_a}}^{2}+ \left ( {\frac {\rm d}{ {\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) \right ) ^{2}+{ \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 \} \]