6 lines
254 B
Text
6 lines
254 B
Text
ClearAll[Faulhaber]
|
|
bernoulliB[1] := 1/2
|
|
bernoulliB[n_] := BernoulliB[n]
|
|
Faulhaber[n_, p_] := 1/(p + 1) Sum[Binomial[p + 1, j] bernoulliB[j] n^(p + 1 - j), {j, 0, p}]
|
|
Table[Rest@CoefficientList[Faulhaber[n, t], n], {t, 0, 9}] // Grid
|
|
Faulhaber[1000, 17]
|