この計算ツールでできること
このツールは、ラゲール陪多項式（一般化ラゲール多項式） \(L_{n}^{(\alpha)}(x)\) を、連続した x の値に対して表形式で算出します。次数 n、パラメータ \(\alpha\)、x の初期値、増分（ステップ幅）、そして生成する行数を指定すると、各 x における多項式の値が一覧で得られます。純粋な数学的計算であり、地域や国による前提は一切含まれず、どこでも同じように利用できます。
使い方
n（0 以上の整数）、\(\alpha\)（任意の実数。標準的な直交性が成り立つのは \(\alpha > -1\) の場合）、x の初期値、増分、行数を入力します。x の値は \(x_i = \text{startX} + i \times \text{stepX}\)（\(i = 0, 1, \ldots, \text{count}-1\)）として順に生成され、それぞれの \(L_{n}^{(\alpha)}(x_i)\) が計算されて表示されます。
計算式の解説
閉じた式は有限和で表され、$$L_{n}^{(\alpha)}(x) = \sum_{k=0}^{n} (-1)^{k} \binom{n+\alpha}{n-k} \frac{x^{k}}{k!}$$ となります。ここで \(\binom{n+\alpha}{n-k}\) は一般化二項係数です。ただし本ツールでは数値計算の安定性を高めるため、3 項間漸化式を用いています。すなわち \(L_{0} = 1\)、\(L_{1} = 1 + \alpha - x\)、そして $$(k+1)L_{k+1} = (2k+1+\alpha-x)L_{k} - (k+\alpha)L_{k-1}$$ です。この方法なら、n が中〜大の場合でも巨大な階乗や桁落ちを避けられます。
計算例
既定値 \(n = 3\)、\(\alpha = 1\) の場合、明示的な多項式は $$L_{3}^{(1)}(x) = 4 - 6x + 2x^{2} - \tfrac{1}{6}x^{3}$$ となります。\(x = 0\) では値は 4。\(x = 0.1\) では \(4 - 0.6 + 0.02 - 0.0001667 \approx 3.419833\)。\(x = 1\) では \(4 - 6 + 2 - 0.166667 = -0.166667\) になります。
第一関連ラゲール多項式
関連（一般化）ラゲール多項式 \(L_n^{(\alpha)}(x)\) は、\(x\) に関する次数 \(n\) の多項式であり、その係数はパラメータ \(\alpha\) に依存します。閉形式は
$$L_n^{(\alpha)}(x)=\sum_{k=0}^{n}(-1)^k\binom{n+\alpha}{n-k}\frac{x^k}{k!}.$$最初の5つは、一般的な \(\alpha\) の形で書くと次のようになります。
| \(n\) | \(L_n^{(\alpha)}(x)\) |
|---|---|
| 0 | \(1\) |
| 1 | \(-x+(\alpha+1)\) |
| 2 | \(\dfrac{x^2}{2}-(\alpha+2)x+\dfrac{(\alpha+1)(\alpha+2)}{2}\) |
| 3 | \(-\dfrac{x^3}{6}+\dfrac{(\alpha+3)x^2}{2}-\dfrac{(\alpha+2)(\alpha+3)x}{2}+\dfrac{(\alpha+1)(\alpha+2)(\alpha+3)}{6}\) |
| 4 | \(\dfrac{x^4}{24}-\dfrac{(\alpha+4)x^3}{6}+\dfrac{(\alpha+3)(\alpha+4)x^2}{4}-\dfrac{(\alpha+2)(\alpha+3)(\alpha+4)x}{6}+\dfrac{(\alpha+1)(\alpha+2)(\alpha+3)(\alpha+4)}{24}\) |
特別な場合 \(\alpha=0\)。\(\alpha=0\) を設定すると、通常のラゲール多項式 \(L_n(x)=L_n^{(0)}(x)\) が得られます。
| \(n\) | \(L_n(x)\) |
|---|---|
| 0 | \(1\) |
| 1 | \(1-x\) |
| 2 | \(1-2x+\tfrac12 x^2\) |
| 3 | \(1-3x+\tfrac32 x^2-\tfrac16 x^3\) |
| 4 | \(1-4x+3x^2-\tfrac23 x^3+\tfrac{1}{24}x^4\) |
最高次項の係数は常に \(\dfrac{(-1)^n}{n!}\) であり、\(\alpha\) に無関係です。
主要な用語および変数
- 次数 \(n\)
- 多項式の次数を与える非負整数。\(L_n^{(\alpha)}(x)\) はちょうど \(n\) 個の根を持ちます。電卓では degree フィールドです。
- パラメータ \(\alpha\)
- 実数（通常 \(\alpha>-1\)）であり、二項係数と直交重み関数をシフトさせます。フィールド alpha。\(\alpha=0\) の場合、多項式は通常のラゲール多項式に縮約されます。
- 引数 \(x\)
- 多項式が評価される点。表は \(x_i=\text{startX}+i\cdot\text{stepX}\) の範囲で変わります。直交性の自然な領域は \((0,\infty)\) です。
- 一般化二項係数
- 実数の上側指数に対して、\(\binom{n+\alpha}{n-k}=\dfrac{\Gamma(n+\alpha+1)}{\Gamma(k+\alpha+1)\,(n-k)!}\) であり、これはガンマ関数を介して \(\binom{m}{j}=m!/(j!(m-j)!)\) を非整数 \(\alpha\) に拡張します。
- 3項漸化式
- 多項式を生成する安定な方法：\((k+1)L_{k+1}^{(\alpha)}=(2k+1+\alpha-x)L_k^{(\alpha)}-(k+\alpha)L_{k-1}^{(\alpha)}\)、\(L_0^{(\alpha)}=1\) および \(L_1^{(\alpha)}=1+\alpha-x\) から始まります。
- \((0,\infty)\) での直交性
- 多項式は相互に直交しています：\(\displaystyle\int_0^\infty L_n^{(\alpha)}(x)L_m^{(\alpha)}(x)\,w(x)\,dx=\frac{\Gamma(n+\alpha+1)}{n!}\delta_{nm}\)。
- 重み関数 \(w(x)=x^{\alpha}e^{-x}\)
- 直交性が成り立つ因子。\(\alpha=0\) の場合、単純な指数重み関数 \(e^{-x}\) です。積分の収束には \(\alpha>-1\) が必要です。
表の解釈
\(L_n^{(\alpha)}(x)\) の計算表を読むことは、次の事実を理解すると簡単になります。
- 実根の個数 \(\alpha>-1\) の場合、\(L_n^{(\alpha)}(x)\) はちょうど \(n\) 個の単純な実ゼロを持ち、すべて開区間 \((0,\infty)\) に含まれます。表の列が \(n\) 回ゼロを通過する場合、すべてのゼロを見つけたことになります。
- 符号の変化 すべてのゼロは単純であるため、多項式は各ゼロで符号が変わります。2つの連続する根の間で、値は一定の符号を保つため、隣接する行の間で符号が反転することは根を括る開始区間となり、二分法や ニュートン法の根探索に役立ちます。
- 原点での値 すべての関連ラゲール多項式は \(L_n^{(\alpha)}(0)=\binom{n+\alpha}{n}=\dfrac{\Gamma(n+\alpha+1)}{n!\,\Gamma(\alpha+1)}\) を満たします。例えば、\(n=4,\ \alpha=0\) の場合、\(x=0\) の最初の行は 1 であり、\(n=4,\ \alpha=2\) の場合は \(\binom{6}{4}=\) 15 です。
- 量子力学 水素原子の波動関数の動径部分は \(L_{n-\ell-1}^{(2\ell+1)}\!\left(2r/(na_0)\right)\) から構成されます。多項式のノードは軌道の動径ノードに対応します。
- ガウス–ラゲール求積法 表に記載されたゼロはまさに \(\int_0^\infty f(x)\,x^{\alpha}e^{-x}\,dx\) を近似するために使用される横座標であり、重みは同じ多項式から導出されます。
これは一般的な数学的参考情報です。重大な応用で依拠する任意の値を検証してください。
よくある質問
n = 0 のときはどうなりますか？ すべての x、すべての \(\alpha\) に対して \(L_{0}^{(\alpha)}(x) = 1\) になります。
\(\alpha\) に負の値や非整数を指定できますか？ はい。和の式も漸化式も任意の実数 \(\alpha\) で機能します。区間 \((0, \infty)\) における古典的な直交性が成り立つのは \(\alpha > -1\) の場合です。
増分を 0 や負の値にできますか？ できます。負の増分では x が小さくなる方向に進みます。増分を 0 にすると同じ x が繰り返され、x が一定の縮退した表になります。
