The harmonic numbers are defined by
Harmonic numbers are sort of a discrete analog of logarithms since
As n goes to infinity, the difference between Hn and log n is Euler’s constant γ = 0.57721… 
How would you compute Hn? For small n, simply use the definition. But if n is very large, there’s a way to approximate Hn without having to do a large sum.
Since in the limit Hn – log n goes to γ, a crude approximation would be
But we could do much better by adding a couple terms to the approximation above.  That is,
The error in the approximation above is between 0 and 1/120n4.
So if you used this to compute the 1000th harmonic number, the error would be less than one part in 120,000,000,000,000. Said another way, for n = 1000 the approximation differs from the exact value in the 15th significant digit, approximately the resolution of floating point numbers (i.e. IEEE 754 double precision).
And the formula is even more accurate for larger n. If we wanted to compute the millionth harmonic number, the error in our approximation would be somewhere around the 26th decimal place.
Related: Posts on special numbers
* * *
 See Julian Havil’s excellent Gamma: Exploring Euler’s Constant. It’s popular-level book, but more sophisticated than most such books.
 There’s a sequence of increasingly accurate approximations that keep adding reciprocals of even powers of n, based on truncating an asymptotic series. See Concrete Mathematics for details.
4 thoughts on “Computing harmonic numbers”
As you talk about “log” above you mean Log_e (log base ‘e’) correct?
To put it another way, you actually man natural log throughout… or ‘ln’.
It is not the normal assumption with log, being log_10.
By “log” I always mean natural log. Here’s why.
What is the exact value of y? the article only gives some digits of precision
γ may be an irrational number, though nobody has proved this.
It’s been computed to over 600,000,000,000 digits so far.