I want to show that there exists a constant $\gamma\in\mathbb{R}$ such that
$$ \sum_{j=1}^N \frac1{j} = \log(N)+\gamma+O(1/N). $$ I know how to prove that the Euler-Mascheroni constant exists (which I believe $\gamma$ to be), but I am having trouble with the big-$O$ notation and the subsequent bounding. I've considered
$$ \left|\left(\sum_{j=1}^N \frac1{j}\right) - \log(N)-\gamma\right|\le |K/N| $$ for some $K$, and I was approaching this by trying to show the that the left side of the inequality decays faster, but so far am stuck. Any advice for this type of problem, or analogous ones, would be appreciated. Thanks!