dsolve(diff(y(x),x) = -(-1-y(x)^4+2*x^2*y(x)^2-x^4-y(x)^6+3*x^2*y(x)^4-3*x^4*y(x)^2+x^6)*x/y(x),y(x))