8
$\begingroup$

Someone asked me for a formula for the sum of the harmonic progression. So I did some calculations and gave him an approximate formula:

$\int_1^n\frac{dx}{x} = \frac{y_1 + y_2}{2} + \frac{y_2 + y_3}{2} + \cdots +\frac{y_{n-1} + y_n}{2}$ where $y_i$ is $i$th term of the HP $\ln(n) = \frac{y_1}{2} + y_2 +y_3 + \cdots +\frac{y_n}{2}$

so

$\sum_{i=1}^n y_i = \ln(n) + \frac{y_1 + y_n}{2}$

e.g. $1+1/2+\cdots+1/10 = 2.8525 $

actual result $= 2.9289$

My question is, how to correct this formula?

3 Answers 3

15

The partial sums of the harmonic series are called "harmonic numbers." The difference between the nth harmonic number and ln(n) tends to a limit as n increases, and that limit is called Euler's constant or gamma.

There's a great book about all this called Gamma: Exploring Euler's Constant.

6

A good quick approximation: $\log(n + 0.5) + \gamma$.

A yet more accurate approximation: $\log(n)+\gamma+1/2n + 1/12n^2-1/120n^4+1/252n^6-1/240n^8+1/132n^{10}-691/32760n^{12}+\cdots$, taking as many terms as is convenient.


Error analysis:

The first formula above has maximum error of $1/24n^2$. uday's first answer has error of about 0.0024 ln n. Mine is a better approximation for all $n\neq6$.

The first formula is also better than uday's revised formula for all $n\ge46$. The expected error in uday's answer is about 0.000016.

1

Hey guys.... There is no correct simple general formula for sum to n terms of the series 1+1/2+1/3+1/4+ ............. + 1/n

but the following formula will be a good approximation for sum to n-terms of the above series when n>5

S = log(n+0.5) + 0.5772 + 0.04021/(n*n + 0.8848)

Deviation from the actual value fluctuates but remains relatively low..

so i guess this may be a good approximation

  • 0
    in my formula upto 45 you may use 0.5772 as the middle term and after 45 replace 2nd and 3rd terms with γ.I guess this process will give you the optimum results for all n.2010-10-25