Limit[Exp[-EulerGamma]*ExpIntegralEi[Log[n] + EulerGamma] - Sum[1/HarmonicNumber[k], {k, 1, n}], {n -> Infinity}]