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