\[ y'(x)=\frac {y(x) \left (x^2+x y(x)+y(x)^2+x\right )}{x^2} \] ✓ Mathematica : cpu = 0.231291 (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}}=x+c_1,y(x)\right ]\] ✓ Maple : cpu = 0.2 (sec), leaf count = 71
\[\left \{y \left (x \right ) = \frac {\sqrt {3}\, x \tan \left (\RootOf \left (2 c_{1} \sqrt {3}+2 \textit {\_Z} +2 \sqrt {3}\, x -\sqrt {3}\, \ln \left (\frac {4}{3 \left (\tan ^{2}\left (\textit {\_Z} \right )+1\right )}\right )-2 \sqrt {3}\, \ln \left (\frac {\tan \left (\textit {\_Z} \right )}{2}-\frac {\sqrt {3}}{6}\right )-\sqrt {3}\, \ln \left (3\right )\right )\right )}{2}-\frac {x}{2}\right \}\]