DSolve[y'[x] + (x + 1)/x y[x] = e^x/x, y[x], x] DSolve[{y'[x] + (x + 1)/x y[x] == e^x/x}, y[x], x]y[x] -> (E^-x C[1])/x + e^x/(x (1 + Log[e]))