|
MATHEMATICA
|
Table[n! SeriesCoefficient[Exp[n - Exp[x] (Exp[n x] - 1)/(Exp[x] - 1)], {x, 0, n}], {n, 0, 16}]
b[n_, k_] := b[n, k] = If[n == 0, 1, -Sum[Binomial[n - 1, j - 1] Sum[i^j, {i, 1, k}] b[n - j, k], {j, 1, n}]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 16}]
|