from math import log

def H(n):
    """Returns an approximate value of n-th harmonic number.
       http://en.wikipedia.org/wiki/Harmonic_number
    """
gamma = 0.5772156649015329
return gamma + log(n) + 0.5 / n - 1.0 / (12 * n ** 2) + 1.0 / (120 * n ** 4)
Программа Python для расчета гармонического ряда