euler mascheroni constant

import math
def eulermascheroni(n):
    parts = range(1, n+1)
    series = []
    for i in parts:
        series.append(1/i)
    sigma = sum(series)
    return sigma - math.log(n)

#try for n = 10000000