scipy.special.bernoulli

scipy.special.bernoulli(n)[源代码]

Bernoulli数B0..Bn(含)

参数
n集成

指示要生成的Bernoulli级数中的项数。

退货
ndarray

伯努利数 [B(0), B(1), ..., B(n)]

参考文献

1

张善杰和金建明。“特殊函数的计算”,约翰·威利父子出版社,1996。https://people.sc.fsu.edu/~jburkardt/f77_src/special_functions/special_functions.html

2

“伯努利数”,维基百科,https://en.wikipedia.org/wiki/Bernoulli_number

示例

>>> from scipy.special import bernoulli, zeta
>>> bernoulli(4)
array([ 1.        , -0.5       ,  0.16666667,  0.        , -0.03333333])

维基百科的文章 ([2]) 指出了伯努利数与Zeta函数之间的关系, B_n^+ = -n * zeta(1 - n)n > 0

>>> n = np.arange(1, 5)
>>> -n * zeta(1 - n)
array([ 0.5       ,  0.16666667, -0.        , -0.03333333])

请注意,在维基百科文章中使用的符号中, bernoulli 计算 B_n^- (即,它使用的惯例是 B_1 为-1/2)。上面给出的关系是用于 B_n^+ ,因此0.5的符号与的输出不匹配 bernoulli(4)