\[ y'(x)=\frac {3 x^2 y(x)-x^3-3 x y(x)^2+y(x)^3+x}{x} \] ✓ Mathematica : cpu = 0.102249 (sec), leaf count = 37
\[\left \{\left \{y(x)\to x-\frac {1}{\sqrt {-2 \log (x)+c_1}}\right \},\left \{y(x)\to x+\frac {1}{\sqrt {-2 \log (x)+c_1}}\right \}\right \}\] ✓ Maple : cpu = 0.21 (sec), leaf count = 57
\[ \left \{ y \left ( x \right ) ={ \left ( x\sqrt {-2\,\ln \left ( x \right ) +2\,{\it \_C1}}-1 \right ) {\frac {1}{\sqrt {-2\,\ln \left ( x \right ) +2\,{\it \_C1}}}}},y \left ( x \right ) ={ \left ( x\sqrt {-2\,\ln \left ( x \right ) +2\,{\it \_C1}}+1 \right ) {\frac {1}{\sqrt {-2\,\ln \left ( x \right ) +2\,{\it \_C1}}}}} \right \} \]