OFFSET
1,1
LINKS
Jwalin Bhatt, Table of n, a(n) for n = 1..126
Simon Plouffe, Identities inspired by Ramanujan Notebooks (part 2), April 2006.
Linas Vepštas, On Plouffe's Ramanujan identities, The Ramanujan Journal, Vol. 27 (2012), pp. 387-408; arXiv preprint, arXiv:math/0609775 [math.NT], 2006-2010.
Eric Weisstein's World of Mathematics, Riemann Zeta Function.
Wikipedia, zeta(2n+1).
FORMULA
EXAMPLE
1470*zeta(5) - 5*Pi^(5) + 3024*S_m(5) + 84*S_p(5) = 0.
18523890*zeta(9) - 625*Pi^(9) + 37122624*S_m(9) + 74844*S_p(9) = 0.
257432175*zeta(13) - 89*Pi^(13) + 514926720*S_m(13) + 62370*S_p(13) = 0.
MATHEMATICA
a[n_] := Module[{s1, s2, f},
s1 = Sum[
(-4)^(n + j) * Binomial[4*n + 2, 4*j] * BernoulliB[4*n - 4*j + 2] * BernoulliB[4*j],
{j, 0, n}
];
s2 = Sum[
(-4)^j * Binomial[4*n + 2, 2*j] * BernoulliB[4*n - 2*j + 2] * BernoulliB[2*j],
{j, 0, 2*n + 1}
];
f = 2^(4*n + 1) * (s1 + s2/2) / Factorial[4*n + 2];
2 * Denominator[f]
]
PROG
(Python)
from sympy import binomial, bernoulli, factorial
def a(n):
s1 = sum(
(-4)**(n+j) * binomial(4*n+2, 4*j) * bernoulli(4*n-4*j+2) * bernoulli(4*j)
for j in range(n+1)
)
s2 = sum(
(-4)**j * binomial(4*n+2, 2*j) * bernoulli(4*n-2*j+2) * bernoulli(2*j)
for j in range(2*n+2)
)
f = 2**(4*n+1) * (s1+s2/2) / factorial(4*n+2)
return 2*f.q
CROSSREFS
KEYWORD
nonn
AUTHOR
Jwalin Bhatt, May 05 2026
STATUS
approved
