"""
Berechnung von Zeta(3)
(c) Jürgen Meier
www.3d-meier.de
05.07.2024
"""

import math


p = 0
k = 5
for n in range(1, k + 1):
    p = p + 4.0/((n**3)*(9*(n**8) + 18*(n**6) + 21*(n**4) + 4))

print(9.0/8 + p)