\[ y'(x)=\frac {y(x) \left (x^7 y(x)^2-3 x^3 y(x)-3\right )}{x \left (x^3 y(x)+1\right )} \] ✓ Mathematica : cpu = 0.134753 (sec), leaf count = 72
DSolve[Derivative[1][y][x] == (y[x]*(-3 - 3*x^3*y[x] + x^7*y[x]^2))/(x*(1 + x^3*y[x])),y[x],x]
\[\left \{\left \{y(x)\to \frac {x}{-x^4+\frac {\sqrt {x+x (-2 x+c_1)}}{\sqrt {\frac {1}{x^7}}}}\right \},\left \{y(x)\to -\frac {x}{x^4+\frac {\sqrt {x+x (-2 x+c_1)}}{\sqrt {\frac {1}{x^7}}}}\right \}\right \}\] ✓ Maple : cpu = 0.034 (sec), leaf count = 38
dsolve(diff(y(x),x) = y(x)/x*(-3*x^3*y(x)-3+y(x)^2*x^7)/(x^3*y(x)+1),y(x))
\[y \left (x \right ) = \frac {1}{x^{3} \left (\sqrt {c_{1}-2 x}-1\right )}\]