|
MATHEMATICA
|
a[n_ /; 1 <= n <= 8] = 1/(n!*2^n); a[n_ /; n > 8] := (k = Quotient[n, 2]; r = Mod[n, 8]; Switch[r, 0, (1 - 2^-k)*(1 + 2^(1-k))/(k!*2)*BernoulliB[k]*Product[BernoulliB[j], {j, 2, 2k-2, 2}], 1|7, (2^k+1)/(k!*2^(2k+1))*Product[BernoulliB[j], {j, 2, 2k, 2}], 2|6, 1/((k-1)!*2^(2k+1))*EulerE[k-1]*Product[BernoulliB[j], {j, 2, 2k-2, 2}], 3|5, (2^k-1)/(k!*2^(2k+1))*Product[BernoulliB[j], {j, 2, 2k, 2}], 4, (1-2^-k)*(1-2^(1-k))/(k!*2)*BernoulliB[k]* Product[BernoulliB[j], {j, 2, 2k-2, 2}], _, Print["error n = ", n]; 0] // Abs); Table[a[n] // Denominator, {n, 1, 30}]
|