I can get,
$$h_n \approx \ln n+\gamma$$
By, saying
$$\frac{h_{n}-h_{n-1}}{n-(n-1)}=\frac{1}{n}$$
Then letting $h_n=h(n)$ be continuous and differentiable on $(0,\infty)$. So that,
$$h'(n) \approx \frac{1}{n}$$
And thus,
$$h_n \approx \ln n+C$$
Where we will let $C$ give us the best approximation as $n \to \infty$ that is,
$$C=\lim_{n \to \infty} (h_n-\ln n)$$
What about the other terms, I'm curious to see how they are derived.
Let $x_n=H_n-(\ln n+\gamma)$ with $H_n=\sum_{k=1}^n\frac{1}{k}$ and evaluate the following limit by using the Stolz-Cesaro Theorem ($\ast/0$ case), $$\lim_{n\to \infty}\frac{x_n}{1/n}=\lim_{n\to \infty}\frac{x_{n+1}-x_n}{1/(n+1)-1/n}=\lim_{n\to \infty}\frac{1/(n+1)-\ln(1+1/n)}{1/(n+1)-1/n}=\frac{1}{2}.$$ Hence $H_n=\ln(n)+\gamma+\frac{1}{2n}+o(1/n)$.
More terms can be found in a similar way (a general formula involves the Bernoulli numbers). Another method to find the asymptotic expansion of the harmonic numbers employs Euler–Maclaurin formula.