I have been asked to calculate how many terms of the series that defines the Euler's constsant $\gamma$ to add at least, to calculate the value of $\gamma$ with error less or equal to $3 \times 10^{-3}$.

I have tried several approaches knowing that the error between the value of a series and a partial addition is defined by

$$ \left| \sum_{k=1}^{\infty}{f(k)} - \sum_{k=1}^{n}{f(k)} \right| = \left| S-S_n \right| \leq \int_{n}^{\infty}{f(t)dt} $$

and just substituting the definition of $\gamma$ as

$$ \gamma = \lim_{n \rightarrow \infty}{\left( \sum_{k=1}^{n}{\dfrac{1}{k} - \log{n}} \right)} $$

but I don't know how to calculate the value of $n$ from there... Any help or suggestion of solving this in other way will be gratefully welcome.


There is the the asymptotic estimation $$ H_n = \log n + \gamma + \frac1{2n} - \frac1{2n^2} + \frac1{120n^4} + O\left(\frac1{n^6}\right) $$ (see Introduction to Analytic Number Theory by Apostol)

EDIT: Wikipedia gives the asymptotic expansion:

$$ H_n = \log n + \gamma + \frac1{2n} - \sum_{k=1}^\infty\frac{B_{2k}}{2kn^{2k}} $$ with $B_k$ the Bernoulli numbers.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.