\[ \boxed { \left ( \left ( y \left ( x \right ) \right ) ^{2}-{x}^{2} \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) ^{2}+{x}^{2}y \left ( x \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) \right ) ^{2}-4\,xy \left ( x \right ) \left ( x{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) -y \left ( x \right ) \right ) ^{3}=0} \]
Mathematica: cpu = 16.501095 (sec), leaf count = 53 \[ \text {DSolve}\left [\left (x^2 y(x) y''(x)+x^2 \left (-y'(x)^2\right )+y(x)^2\right )^2-4 x y(x) \left (x y'(x)-y(x)\right )^3=0,y(x),x\right ] \]
Maple: cpu = 0.343 (sec), leaf count = 86 \[ \left \{ y \left ( x \right ) ={\it ODESolStruc} \left ( {{\rm e}^{\int \!{\it \_b} \left ( {\it \_a} \right ) \,{\rm d}{\it \_a}+{\it \_C1}}},[ \left \{ {\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) =2\,{\frac {\sqrt {{\it \_a}\, \left ( {\it \_a}\,{\it \_b} \left ( {\it \_a} \right ) -1 \right ) }{\it \_b} \left ( {\it \_a} \right ) }{{\it \_a}}}-{\frac {1}{{{\it \_a}}^{2}} \left ( 2\,\sqrt {{ \it \_a}\, \left ( {\it \_a}\,{\it \_b} \left ( {\it \_a} \right ) -1 \right ) }+1 \right ) } \right \} , \left \{ {\it \_a}=x,{\it \_b} \left ( {\it \_a} \right ) ={\frac {{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) }{y \left ( x \right ) }} \right \} , \left \{ x={\it \_a},y \left ( x \right ) ={{\rm e}^{\int \!{\it \_b} \left ( {\it \_a} \right ) \,{\rm d}{\it \_a}+{\it \_C1}}} \right \} ] \right ) \right \} \]