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

import math


p = 0
k = 3
for n in range(0, k + 1):
    a = 126392*n**5 + 412708*n**4 + 531578*n**3 + 336367*n**2 + 104000*n + 12463
    zaehler = ((-1.0)**n)*a*((math.factorial(2*n+1)*math.factorial(2*n)*math.factorial(n))**3)
    nenner = math.factorial(3*n+2)*(math.factorial(4*n+3)**3)

    p = p + zaehler/nenner

print(p/24)