Pythonで与えられた根からルジャンドル級数を生成する方法

Pythonでルジャンドル級数（Legendre series）を生成するには、numpy.polynomial.legendreモジュールのlegfromroots()メソッドを使用します。このメソッドは、指定された根（roots）に対応するルジャンドル級数の係数を1次元配列として返します。

legfromroots()メソッドの特徴

• 戻り値は係数の1次元配列です。
• すべての根が実数の場合、結果は実数型の配列になります。
• 一部の根が複素数の場合、たとえ計算結果の係数がすべて実数であっても、出力は複素数型の配列になります。
• 引数rootsには、級数の根となるシーケンス（タプルやリストなど）を渡します。

手順

1. 必要なライブラリをインポートする

まず、NumPyとルジャンドル多項式モジュールをインポートします。

import numpy as np
from numpy.polynomial import legendre as L

2. ルジャンドル級数を生成する

L.legfromroots()メソッドに根を渡して、ルジャンドル級数を生成します。ここでは例として、根 (-1, 0, 1) を使用します。

print("Result...\n", L.legfromroots((-1, 0, 1)))

3. データ型を確認する

print("\nType...\n", L.legfromroots((-1, 0, 1)).dtype)

4. 配列の形状を確認する

print("\nShape...\n", L.legfromroots((-1, 0, 1)).shape)

完全なコード例

import numpy as np
from numpy.polynomial import legendre as L

# legfromroots()メソッドでルジャンドル級数を生成
print("Result...\n", L.legfromroots((-1, 0, 1)))

# データ型を取得
print("\nType...\n", L.legfromroots((-1, 0, 1)).dtype)

# 配列の形状を取得
print("\nShape...\n", L.legfromroots((-1, 0, 1)).shape)

実行結果

Result...
    [ 0.  -0.4  0.   0.4]

Type...
float64

Shape...
(4,)

解説

この例では、根 -1、0、1 を持つルジャンドル級数が生成されました。出力された係数配列 [ 0. -0.4 0. 0.4] は、次数0から3までのルジャンドル多項式の線形結合を表しています。

データ型は float64（64ビット浮動小数点数）であり、入力された根がすべて実数のため、結果も実数型になっています。また、配列の形状は (4,) で、これは4つの係数（次数0〜3に対応）が含まれていることを意味します。

なお、根の数が n 個の場合、生成されるルジャンドル級数の次数は n となり、係数配列の長さは n+1 になります。この挙動は通常の多項式における「根から多項式を構築する」操作と同じ原理ですが、係数がルジャンドル多項式基底に対して表現される点が異なります。