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