\[ y'(x)=-x^6+3 x^4 y(x)-3 x^2 y(x)^2+y(x)^3+2 x \] ✓ Mathematica : cpu = 0.083613 (sec), leaf count = 39
DSolve[Derivative[1][y][x] == 2*x - x^6 + 3*x^4*y[x] - 3*x^2*y[x]^2 + y[x]^3,y[x],x]
\[\left \{\left \{y(x)\to x^2-\frac {1}{\sqrt {-2 x+c_1}}\right \},\left \{y(x)\to x^2+\frac {1}{\sqrt {-2 x+c_1}}\right \}\right \}\] ✓ Maple : cpu = 0.039 (sec), leaf count = 57
dsolve(diff(y(x),x) = y(x)^3-3*x^2*y(x)^2+3*y(x)*x^4-x^6+2*x,y(x))
\[y \left (x \right ) = \frac {x^{2} \sqrt {2 c_{1}-2 x}-1}{\sqrt {2 c_{1}-2 x}}\]