"""
Werte der Bernoulli Konstante
(c) Jürgen Meier
www.3d-meier.de
07.07.2024
"""

import math


    
n = 8

b = 0
    
for k in range(0, n+1):
    a = 0
    for j in range(0, k+1):    
        a = a + ((-1.0)**j) * math.factorial(k)*(j**n)/(math.factorial(j) * math.factorial(k-j))

    b = b + a/(k + 1)
print(b)