The relationship between each harmonic numbers In Knuth's "Concrete Mathematics" in chapter about numbers below equality is given
$$H_n = \ln n + \gamma + \frac{1}{2n} - \frac{1}{12n^2} + \frac{\epsilon_n}{120n^4} $$ where $0 < \epsilon_n < 1$.
For every $i$ holds $H_{i-1} < H_i$. Is this also true that for every $i$ holds: 
$$\frac{\epsilon_i}{120i^4} < \frac{\epsilon_{i-1}}{120(i-1)^4}$$. I was trying to prove that but I can't prove the dependency between $\epsilon_i$ and $\epsilon_{i-1}$. Where can I found a prove or a contrexample?

If you persue the approximation of $H_n$ a bit further you may find
$$ H_n=\ln n+\gamma+\frac1{2n}-\frac1{12n^2}+\frac1{120n^4}-\frac1{252n^6}+\frac1{240n^8}-\frac1{132n^{10}}+\mathcal O(\frac1{n^{12}}).$$
Hence for some constant $c$ (not explicitly specified in that Wikipedia page), we have
$$ \epsilon_n=1-\frac1{252n^2}+\frac1{240n^4}-\frac1{132n^6}+\frac{\delta_n}{cn^8}$$
with $|\delta_n|<1$. Depending on the ecaxt value of $c$ this gives us that the sequence of the $\epsilon_n$ is strictly increasing either for all sufficiently large $n$ or possibly even for all $n$. Indeed, $e_{n+1}-\epsilon_n\approx\frac1{126n^3}$ already by the first term and the contribution of the last is at most $\approx \frac2{cn^8}$.
