\[ \boxed { {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) ={\frac {y \left ( x \right ) \left ( \left ( y \left ( x \right ) \right ) ^{2}+xy \left ( x \right ) +{x}^{2}+x \right ) }{{x}^{2}}}=0} \]
Mathematica: cpu = 0.048006 (sec), leaf count = 60 \[ \text {Solve}\left [-\frac {1}{2} \log \left (\frac {y(x)^2}{x^2}+\frac {y(x)}{x}+1\right )+\log \left (\frac {y(x)}{x}\right )-\frac {\tan ^{-1}\left (\frac {\frac {2 y(x)}{x}+1}{\sqrt {3}}\right )}{\sqrt {3}}=c_1+x,y(x)\right ] \]
Maple: cpu = 0.125 (sec), leaf count = 71 \[ \left \{ y \left ( x \right ) ={\frac {\sqrt {3}x}{2}\tan \left ( {\it RootOf} \left ( -\sqrt {3}\ln \left ( 3 \right ) -\sqrt {3}\ln \left ( { \frac {4}{3+3\, \left ( \tan \left ( {\it \_Z} \right ) \right ) ^{2}}} \right ) -2\,\sqrt {3}\ln \left ( -1/6\,\sqrt {3}+1/2\,\tan \left ( { \it \_Z} \right ) \right ) +2\,\sqrt {3}{\it \_C1}+2\,\sqrt {3}x+2\,{ \it \_Z} \right ) \right ) }-{\frac {x}{2}} \right \} \]