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