4 lines
203 B
Text
4 lines
203 B
Text
ClearAll[Faulhaber]
|
|
Faulhaber[n_, 0] := n
|
|
Faulhaber[n_, p_] := n^(p + 1)/(p + 1) + 1/2 n^p + Sum[BernoulliB[k]/k! p!/(p - k + 1)! n^(p - k + 1), {k, 2, p}]
|
|
Table[{p, Faulhaber[n, p]}, {p, 0, 9}] // Grid
|