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