Plot[{HarmonicNumber[x], EulerGamma+Gamma'(x+1)/Gamma(x+1)}, {x, -1, 5}]