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