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