MCPで接続 →

計算を入力してください

公式

広告

結果

Bernoulli Number B4 (exact fraction)
-1/30
convention B1 = -1/2
インデックス n 4
厳密な分数 -1/30
Decimal (22 sig. digits) -0.03333333333333333333333

ベルヌーイ数とは

ベルヌーイ数 \(B_n\) は、整数論や解析学の至るところに登場する有理数の数列です。母関数 \(\frac{x}{e^x - 1}\) のマクローリン展開における係数として定義されます。整数のべき乗和の閉じた式、オイラー=マクローリンの公式、偶数における リーマンゼータ関数の値、さらに正接(tan)や余接(cot)関数のテイラー展開など、さまざまな場面で現れます。

最初のベルヌーイ数を正確な分数で示した表
最初の数個のベルヌーイ数を正確な分数で表示。奇数番目(B1より上)の値はゼロ。

この計算機の使い方

非負整数のインデックス \(n\)(0, 1, 2, 3, …)を入力し、表示する有効桁数を選びます。計算機は厳密な有理数を分数 \(p/q\) の形(例:\(B_{12} = -\frac{691}{2730}\))で返すとともに、丸めた小数値も表示します。有効桁数の設定は小数表示にのみ影響し、分数は常に厳密な値です。

採用している規約

本ツールでは母関数 \(\frac{x}{e^x - 1}\) に対応する規約 \(B_1 = -\frac{1}{2}\) を採用しています。もう一方の「プラス」規約では \(B_1 = +\frac{1}{2}\) とし、母関数 \(\frac{x}{1 - e^{-x}}\) を用いますが、両者の違いは \(B_1\) の符号だけです。それ以外のベルヌーイ数はどちらの規約でも同一であり、\(B_1\) より先の奇数インデックスの値はすべて厳密に 0 になります。

計算式の解説

浮動小数点による丸め誤差(たとえば \(B_2\) が 0.16666… となり、\(6 \cdot B_2\) が 0 に切り捨てられてしまうなど)を避けるため、この計算機は厳密な有理数演算を用いています。漸化式 $$B_0 = 1, \quad B_m = -\frac{1}{m+1}\sum_{k=0}^{m-1}\binom{m+1}{k}\,B_k$$ を適用します。ここで \(\binom{m+1}{k}\) は二項係数です。各 \(B_k\) は約分された分子・分母の組として保持されるため、小数に変換する前の答えは数学的に厳密です。

べき級数に展開したベルヌーイ数の母関数
べき級数の係数がベルヌーイ数Bₙを定義する母関数 x/(eˣ−1)。

計算例(n = 4)

\(B_0 = 1\)、\(B_1 = -\frac{1}{2}\)、\(B_2 = \frac{1}{6}\)、\(B_3 = 0\) から始めると、漸化式により $$B_4 = -\frac{1}{30} \approx -0.0333333\ldots$$ が得られます。既知の表と照合して確かめることもできます:\(B_6 = \frac{1}{42}\)、\(B_8 = -\frac{1}{30}\)、\(B_{10} = \frac{5}{66}\)、\(B_{12} = -\frac{691}{2730}\)。

よくある質問

なぜ奇数番目のベルヌーイ数は 0 なのですか? \(B_1\) を除き、奇数インデックスのベルヌーイ数 \(B_{2n+1}\) はすべて 0 になります。これは母関数の対称性によるものです。

なぜ偶数インデックスの値は大きくなるのですか? その大きさは急速に増大します。たとえば \(|B_{50}|\) はおよそ \(7.5\times10^{24}\) です。厳密な分数を用いれば、こうした巨大な値もオーバーフローせずに扱えます。

ここでの \(B_1\) は正ですか負ですか? 負です。この計算機は \(B_1 = -\frac{1}{2}\) を返します。

最終更新: