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