正の実数のログガンマ関数 ln Γ(a) を計算

使い方

1. 0より大きい実数aを入力し、有限値の宣言範囲と定義域を満たすことを確認します。
2. 求めたい量がΓ(a)ではなく自然対数ln Γ(a)であることを確認してから「計算する」を押し、正の実数のログガンマln Γ(a)の値を表示します。
3. ln Γ(1)=0や漸化式の差分と照らし、符号、桁、極や飽和への近さを読み取ります。
4. 再利用するときは入力値、表示桁、確認日、実数枝という条件も一緒に保存します。

計算式・考え方

a≥1/2ではlog Γ(a)=log√(2π)+(a-1/2)log(a+6.5)-(a+6.5)+log A(a)のLanczos形を使います。0<a<1/2ではlog Γ(a)=log Γ(a+1)-log aで安定な正側へ移します。

注意事項

• 正の実数のログガンマln Γ(a)は、正の実数で正値を取るΓ(a)の自然対数を、Γ(a)そのものへ戻さず直接評価する関数です。このページは一つまたは二つの実数入力から一つの実数値を返すスカラー計算に限定しています。グラフ、複素数、数列の一括生成は扱わず、入力した一点の値と定義域を落ち着いて確認するためのページです。
• 採用した数式は、正側のLanczos対数形とΓ(a+1)=aΓ(a)の対数log Γ(a)=log Γ(a+1)-log aを使います。表示直前だけで丸め、内部ではJavaScriptとPythonが共通に扱えるIEEE 754倍精度の有限値を保ちます。式を表計算へ転記するときも、括弧、指数、定数の位置を変えないでください。
• 入力範囲は、aを0より大きく1,000,000以下に限定し、負側のlog|Γ(a)|や複素対数は対象外にします。数学的な定義域そのものより狭い場合は、関数を別物へ縮約したのではなく、有限な倍精度結果をブラウザで再現できる契約として上限と下限を設けています。空欄、文字、NaN、無限大、範囲外の数は計算前に拒否します。
• 数値安定性については、巨大なΓ(a)をexpで生成しないため、ガンマ関数ページより広い正の入力範囲でも有限な対数値を保ちます。単純な教科書式をそのまま評価すると、近い大きな数の差、指数オーバーフロー、ゼロ除算が起こる場所があります。この実装は条件分岐を遅延評価し、選ばれなかった危険な枝を実行しません。
• 境界の読み方は、a=1ではln Γ(1)=0、a=1/2ではln Γ(1/2)=ln√πで、a→0+では正の無限大へ発散します。境界値を自動的にゼロへ置き換えたり、発散点へ便宜的な有限値を付けたりしません。定義できる特別値だけを明示的に返し、極や許容範囲外は入力エラーとして理由を示します。
• 結果は、積や階乗比を対数和へ変換するときに使え、値が負でもΓ(a)自体は正です。最後の数桁は処理系、入力の十進表現、近似式の丸めによって変わり得ます。手計算、別ライブラリ、論文の表と比べるときは、有効桁だけでなく定義、引数の並び、正規化、実数枝もそろえてください。
• 検算では、ln Γ(1)=0、ln Γ(5)=ln24、ln Γ(a+1)-ln Γ(a)=ln aを正の複数点で確認します。通常点、ゼロまたは既知の特別値、定義域外の三種類を分けて試験しています。対称性や漸化式がある関数では、一点の一致だけでなく、符号を反転した値や隣接する引数との関係も独立に確かめます。
• 参照資料は、NIST DLMF 5.2と5.5、Boost.MathのLanczos近似資料を参照し、自然対数であることを確認しました。2026年7月13日に定義、実数領域、特異点、近似式の係数を照合しました。NIST DLMFは数式の節番号、Boost.MathはLanczos近似の実装上の意味を確認するために使い、特定ライブラリの出力を固定値として埋め込んではいません。
• このページで扱わない範囲は、a≤0の符号付きガンマ、log|Γ(a)|、複素対数、任意精度の超巨大引数です。目的外の入力を無理に通すと、見かけ上は数値でも別の枝や別の規約を表すことがあります。必要な領域が範囲外なら、複素特殊関数や任意精度演算に対応した検証済みライブラリを選んでください。
• 実務への持ち込み方は、確率密度の正規化や尤度の積を対数和で扱う場合は、底10の対数ではなく自然対数であることを記録します。この計算値は学習、式の確認、テストデータ作成を助けるものです。観測モデル、材料モデル、統計分布、最適化手順の妥当性は別問題なので、前提、単位、測定誤差、採用したモデル名を計算結果と一緒に記録してください。
• 有限な実数とIEEE 754倍精度による学習・検算用の近似値です。複素枝、任意精度、物理モデル全体の妥当性は扱いません。

よくある質問