ルジャンドル陪多項式（グラフ）計算機とは
この計算機は、ルジャンドル陪関数 \(P_n^m(x)\)（次数 \(n\)、位数 \(m\)）の値を指定した \(x\) の範囲で数表として求め、対応する曲線を描画します。純粋な数学であり、単位や国・地域に依存する前提は一切なく、どこでも同じ結果が得られます。ルジャンドル陪多項式は物理学や応用数学の至るところに登場し、球面調和関数、球座標におけるラプラス方程式の解、多重極展開、角運動量の量子力学などで重要な役割を果たします。
使い方
整数の次数 \(n\)（0, 1, 2, …）と整数の位数 \(m\)（\(-n \le m \le n\) を満たす値）を入力します。続いて \(x\) の開始値（\(-1\) から \(1\) の範囲）、刻み幅、行数を指定します。初期設定の \(n = 2\)、\(m = 1\)、開始値 \(= -1\)、刻み幅 \(= 0.02\)、行数 \(= 101\) では、\(x\) を \(-1\) から \(+1\)（両端を含む）まで走査します。符号規約は Type A（Wolfram 規約）または Type B（Maple 規約）から選べます。区間 \((-1, 1)\) の実数 \(x\) においては両者の絶対値は一致し、前因子の符号・位相だけが異なります。
計算式の解説
整数 \(n\) と \(0 \le m \le n\) に対して、$$P_n^m(x) = (-1)^m\,(1-x^2)^{m/2}\,\frac{d^m}{dx^m}P_n(x)$$ を用います。実際の計算には数値的に安定な漸化式を使い、まず \(P_m^m = (-1)^m(2m-1)!!(1-x^2)^{m/2}\)、\(P_{m+1}^m = x(2m+1)P_m^m\) を求め、次に $$(l-m)P_l^m = (2l-1)x\,P_{l-1}^m - (l+m-1)P_{l-2}^m$$ によって次数を上げていきます。負の \(m\) に対しては $$P_n^{-m} = (-1)^m\frac{(n-m)!}{(n+m)!}P_n^m$$ となります。この閉じた漸化式を用いることで、正の整数 \(m\) に対して \({}_2F_1\) 形式をそのまま使ったときに生じるガンマ関数の発散を回避できます。
計算例
\(n = 2\)、\(m = 1\) のとき、関数は $$P_2^1(x) = -3x\sqrt{1-x^2}$$ となります。\(x = 0\) では値は \(0\)、\(x = 0.5\) では \(-3(0.5)(0.866025) = -1.299038\)、\(x = -0.5\) では \(+1.299038\) です。曲線は \(x = -1\) で \(0\) から始まり、\(x \approx -0.577\) 付近で約 \(+1.1547\) まで上昇し、\(x = 0\) で \(0\) を横切り、\(x \approx +0.577\) 付近で約 \(-1.1547\) まで下降した後、\(x = +1\) で再び \(0\) に戻ります。
閉形式関連ルジャンドル関数 P_n^m(x)
整数次数 \(n\) と階数 \(0\le m\le n\) の関連ルジャンドル関数 \(P_n^m(x)\) は、\(P_n^m(x)=(-1)^m(1-x^2)^{m/2}\dfrac{d^m}{dx^m}P_n(x)\) から導かれます。係数 \((-1)^m\) は Condon–Shortley 位相で、Type A 規約（Wolfram に準拠）に含まれています。Type B 規約（Maple）はこれを省略するため、奇数次の \(m\) の項は符号のみ異なります。以下の表は Type A による明示的な形式を示しています。
| \(n\) | \(m\) | \(P_n^m(x)\)（Type A、符号付き） |
|---|---|---|
| 0 | 0 | \(1\) |
| 1 | 0 | \(x\) |
| 1 | 1 | \(-\sqrt{1-x^2}\) |
| 2 | 0 | \(\tfrac{1}{2}(3x^2-1)\) |
| 2 | 1 | \(-3x\sqrt{1-x^2}\) |
| 2 | 2 | \(3(1-x^2)\) |
| 3 | 0 | \(\tfrac{1}{2}(5x^3-3x)\) |
| 3 | 1 | \(-\tfrac{3}{2}(5x^2-1)\sqrt{1-x^2}\) |
| 3 | 2 | \(15x(1-x^2)\) |
| 3 | 3 | \(-15(1-x^2)^{3/2}\) |
確認として、\(x=0.5\) での P_2^1 の項は \(-3(0.5)\sqrt{1-0.25}=-1.5\sqrt{0.75}=\) -1.299038 です。\(m=0\) 列は通常のルジャンドル多項式 \(P_n(x)\) を再現します。例えば \(P_3^0(x)=\tfrac12(5x^3-3x)\) は、ルジャンドル多項式テーブル計算機で表形式にできます。
主要な用語と変数
- 次数 \(n\)
-
非負整数（
degreeN）で、基礎となるルジャンドル多項式 \(P_n(x)\) の階数を設定します。これは次数 \(n\) の多項式です。 - 階数 \(m\)
-
整数（
orderM）で、取る微分の個数を制御します。\((-1,1)\) での実数値結果の場合、通常 \(0\le m\le n\) を使用します。\(m>n\) の場合、次数 \(n\) の多項式の \(m\) 次微分は消えるため、関数は恒等的にゼロです。 - 引数 \(x\)
-
評価点（
initialXプラス \(i\cdot\)stepX）です。関数は \(-1\le x\le 1\) で実数です。物理学では \(x=\cos\theta\) です。 - Type A（Wolfram / Condon–Shortley）
-
位相係数 \((-1)^m\) を含みます。これは Wolfram の
LegendrePと標準的な量子力学の教科書で使用される規約です。 - Type B（Maple）
- \((-1)^m\) 位相を省略するため、P_n^m（Type B）\(=(-1)^m\,P_n^m\)（Type A）です。大きさは同じで、奇数次の \(m\) の項の符号のみが異なります。
- 二重階乗 \((2m-1)!!\)
- 奇数の積 \((2m-1)(2m-3)\cdots 3\cdot 1\) で、\((-1)!!=1\) です。主要係数 \(P_m^m(x)=(-1)^m(2m-1)!!\,(1-x^2)^{m/2}\) に現れます。例えば \(P_3^3\) は \(5!!=15\) を使用します。これらの値については二重階乗計算機を参照してください。
- 負の階数関係
- \(m>0\) について、\(P_n^{-m}(x)=(-1)^m\dfrac{(n-m)!}{(n+m)!}\,P_n^{m}(x)\) は、階乗を通じて正と負の階数を結びます。
表とグラフの解釈
いくつかの構造的性質により、表示された値とプロットされた曲線をサニティチェックできます：
- パリティ。 \(P_n^m(-x)=(-1)^{n+m}P_n^m(x)\)。\(n+m\) が偶数の場合、グラフは \(x=0\) について対称です。\(n+m\) が奇数の場合、反対称（したがって原点を通過）です。
- 内部のゼロ。 開区間 \((-1,1)\) では、\(P_n^m(x)\) は正確に \(n-m\) 個の単純なゼロを持ちます。例えば \(P_3^1\) は 2 つの内部ゼロを持ち、\(P_n^n\) は持ちません。
- 端点の挙動。 係数 \((1-x^2)^{m/2}\) があるため、\(m>0\) のすべての関数は \(x=\pm 1\) で消えます。\(m=0\) の場合、値は \(P_n(1)=1\) と \(P_n(-1)=(-1)^n\) です。
- 端付近の大きさ。 より高い \(m\) では、\((1-x^2)^{m/2}\) 係数は \(x\to\pm1\) のときに曲線を急激に抑制するため、最大の変動は範囲の中央に向かって発生します。
これらの関数は、球面調和関数 \(Y_n^m(\theta,\phi)\) の \(\theta\) 依存部分です。\(x=\cos\theta\) と書くと、\(Y_n^m\propto P_n^m(\cos\theta)\,e^{im\phi}\) になります。内部ゼロは緯度の節円となり、\(m>0\) の端点消失は調和関数が極で傾向するゼロに対応します。したがって同じ \(P_n^m\) の値は、選択された \(\theta\) と \(\phi\) での球面調和関数の評価に直接供給されます。
よくある質問
なぜ \(n\) と \(m\) は整数でなければならないのですか。 多項式として有限項で終わる形になるためには、\(n\) が非負整数である必要があります。また漸化式や \((n\pm m)!\) の因子を扱うには、\(m\) が \(-n \le m \le n\) を満たす整数であることが必要です。
表示されるサンプル値とは何ですか。 上部のボックスには、数表の中央の行（中央のインデックス）における \(x\) と \(P_n^m(x)\) の値が表示されます。曲線を手早く確認するための目安です。
\(m = 0\) のときはどうなりますか。 \(P_n^0(x)\) は通常のルジャンドル多項式 \(P_n(x)\) と一致します。
