DSolve[Derivative[1][y][x] == (-(E^(1 + x)*x^3) - y[x] + x*y[x] + E^(1 + x)*x*y[x]^2)/((-1 + x)*x),y[x],x]