積分からBeta関数を計算

使い方

1. ページ先頭の計算フォームで関数式と数値条件を確認し、暗黙の掛け算、未対応名、定義域外がないよう入力します。
2. 最大ステップ・最大段階・次数と、必要な基礎則を選び、計算ボタンでブラウザ内の収束表を更新します。
3. 複数行の推定値を既知の恒等式または別の積分器と比較し、隣接値の一致だけで精度を断定しません。
4. 列見出し付きTSVをコピーし、関数式、区間、パラメータ、方式、最大12有効桁表示、確認日を一緒に保存します。

計算式・考え方

B(a,b)=∫_0^1x^(a−1)(1−x)^(b−1)dx=Γ(a)Γ(b)/Γ(a+b)（a,b>0）です。選択区間では同じ核の部分積分を計算します。

注意事項

• 既定値a=3、b=5、区間[0,1]、次数20では、19次が約0.009523809523809521、20次が約0.009523809523809528です。核は次数6の多項式x^2(1−x)^4なので、十分なGauss-Legendre次数では丸め誤差の範囲で厳密に積分でき、閉形式B(3,5)=1/105と照合できます。二行の差が小さい理由を、一般関数の精度保証へそのまま拡張しないでください。
• aとbは0.05以上100以下へ制限します。数学上は両方が正ならEuler積分が収束しますが、1未満では対応する端点で核が発散し、0へ近いほど端点特異性が強くなります。Gauss-Legendre分点は端点そのものを評価しないものの、重みが特異性へ特化していないため収束が遅い場合があります。必要ならGauss-Jacobi、Tanh–Sinh、変数置換、多倍長積分と比較します。
• 完全なBeta関数として解釈できるのは区間[0,1]です。区間を狭めれば同じ核の部分積分となり、正規化すれば不完全Beta関数に関連します。逆向き区間では向き付き積分として符号が反転します。0より外側では非整数指数が実数にならない場合があるため、この定義は端点入力を0から1へ制限し、複素数への分岐を扱いません。
• 対称性B(a,b)=B(b,a)は独立した検算に使えます。aとbを交換し、同時にxを1−xへ置換すると完全区間の積分値は変わりません。またGamma関数との関係は、積分器とは異なるアルゴリズムで値を確かめる手段です。大きいパラメータでは核が非常に狭くなるため、Gamma比は対数領域で計算し、次数を増やしてピークを捉えたか確認してください。
• 数値積分は、有限個の標本点で被積分関数を評価し、その重み付き和から積分を近似します。解析的な原始関数が得られない場合にも役立ちますが、有限個の値だけから積分の存在や表示桁すべての正しさを自動証明する機能ではありません。結果表では段階または次数を変えた複数の推定値を比較し、既知の恒等式、別方式、多倍長計算、概算面積など独立した根拠を併用してください。
• 関数式には変数x、定数piとe、四則演算、べき乗^、丸括弧、abs、sqrt、exp、ln、log、log10、sin、cos、tan、asin、acos、atan、sinh、cosh、tanh、floor、ceilを使えます。2xのような暗黙の掛け算は使えないため2*xと書きます。式は200文字、128トークン、演算子と括弧の深さ32段までに制限し、未知の名前や不完全な式は計算開始前に拒否します。
• Gamma・Betaページのaとbは、画面の数値入力から予約済みの式パラメータへ束縛します。利用者が式の中へ任意の変数名やコードを追加する方式ではありません。字句解析と構文解析で許可した演算だけを逆ポーランド記法へ変換して評価し、eval、new Function、プロパティ参照、配列、代入、セミコロン、文字列を実行しません。パラメータ入力も有限値と所定範囲を検証してから式へ渡します。
• 計算は利用者のブラウザ内で完結します。外部API、CDN、fetch、XMLHttpRequest、WebSocket、サーバー側計算を必要とせず、入力した関数式、区間、パラメータ、結果を送信または保存しません。サーバーの計算資源を消費せずに条件を変更できますが、端末を閉じれば入力は残りません。必要な結果は式、条件、方式、確認日を添えて列見出し付きTSVへ保存してください。
• 内部の実数演算はIEEE 754 binary64で、画面とTSVは最大12有効桁に整形します。参照資料や多倍長計算が20桁以上を示しても、同じ桁数を転載してbinary64以上の精度があるようには見せません。大きさが異なる項には安定した総和を使いますが、桁落ち、正負相殺、オーバーフロー、アンダーフローを完全には排除できません。高精度定数表の作成には任意精度実装が必要です。
• 収束表の隣接値が近いことは、二つの有限近似が一致したという観察です。真値との差、厳密な誤差上界、正しい桁数の保証と同じではありません。標本点の間に狭いピークがある場合、複数段階が同じ構造を見落として見かけ上安定することがあります。段階や次数を変えるだけでなく、区間分割、変数変換、関数のグラフ、解析的な上限下限との照合も行ってください。
• 積分区間内の極、不連続、平方根や対数の定義域境界は事前に確認します。分点が極へ偶然一致しないだけでは通常積分が存在する根拠にならず、Cauchy主値を自動的に選ぶこともありません。既知の内部特異点では区間を分け、端点特異性には問題に合う変数変換やTanh–Sinhを検討します。関数評価や重み付き和が有限値にならなければ、値らしい数を補間せずエラーにします。
• 滑らかな関数でも、振動周波数が高い、変化の尺度が区間幅より極端に小さい、非常に大きな正負の寄与が相殺する、といった場合は通常の収束傾向が現れません。最高設定へ一度に上げるのではなく、低い設定から履歴を見て、結果が安定する領域と丸め誤差が支配する領域を区別します。正値関数の積分が負になるなど基本的な性質に反する場合は採用しません。
• 単位を持つ量へ利用するときは、積分変数と被積分関数の単位の積が結果の単位になります。Gamma関数とBeta関数の定義式は無次元変数を前提にしているため、物理量の密度や確率分布へ転用するときは尺度変換と正規化定数を別途確認してください。数値だけが既知値に近くても、変数、端点、指数パラメータの意味が違えば同じ問題を解いたことにはなりません。
• TSVコピーは表の列見出しと表示順を含みます。再現性のため、Rombergでは基礎則と最大ステップ、Gammaではa・有限下端・最大段階、Betaではa・b・区間・次数を必ず残します。binary64・最大12有効桁であること、隣接値の差が誤差保証ではないことも併記します。結果の一列だけを資料へ移すと、どの公式と入力で得た値か後から検証できません。
• この計算機は数値解析の学習、プログラム実装の照合、一般的な試算を目的とします。医療投与量、安全限界、法令判定、金融取引など高い保証が必要な判断を単独で決めるものではありません。重要な用途では式とデータの出典、許容誤差、独立実装との一致をレビューし、必要に応じて専門家と証明付き・多倍長の計算環境を利用してください。
• 求積公式、Gaussian求積、外挿法の一般的な定義と注意点はNIST Digital Library of Mathematical Functions（DLMF）§3.5を2026年7月13日に照合しました。GammaのEuler積分は同資料§5.2、BetaのEuler積分は§5.12と照合しています。本文では各定義の実数正パラメータ域だけを扱い、解析接続、複素積分、極の回避、正則化を実装範囲に含めません。
• 半無限区間に用いる二重指数型変換については、高橋秀俊・森正武のTakahasi–Mori公式と京都大学数理解析研究所（RIMS）で公開された関連資料を2026年7月13日に確認しました。変換後の有限台形和は十分に減衰する関数で有効ですが、方式名だけで収束を保証しません。遅いべき減衰、発散積分、無限に続く振動には問題専用の尾部処理が必要です。
• binary64・最大12有効桁の有限近似です。収束差や最終対角値は厳密な誤差保証ではありません。

よくある質問