ルジャンドルの微分方程式

数学の解析学分野に於ける「ルジャンドルの微分方程式」とは次の式で表わされる線形微分方程式である。；

• $(1-x^2)y^{″}-2x{y}^{\prime }+\nu (\nu +1)y=0$

(※ただし$\nu \in \mathbb{ℝ}$とする。)

これはスツルム＝リウヴィル型微分方程式とゆー微分方程式の一例とされている。

級数解法

上記微分方程式は冪級数の項別微分を用いると解を求める事ができる。で、その無限冪級数解を

$y=\sum _{m=0}^{\mathrm{\infty }}{c}_m{x}^m$

とおくと項別微分の定理により

$y^{\prime }=\sum _{m=1}^{\mathrm{\infty }}m{c}_m{x}^{m-1},y^{″}=\sum _{m=2}^{\mathrm{\infty }}m(m-1)c_m{x}^{m-2}$

が得られる。これをルジャンドルの微分方程式に代入すれば方程式の左辺は $\sum _{m=2}^{\mathrm{\infty }}m(m-1)c_m{x}^{m-2}-\sum _{m=2}^{\mathrm{\infty }}m(m-1)c_m{x}^m-2\sum _{m=1}^{\mathrm{\infty }}m{c}_m{x}^m+\nu (\nu +1)\sum _{m=0}^{\mathrm{\infty }}{c}_m{x}^m$ と変形できる。

ここで上述の微分方程式は変数xに対して右辺が零の恒等式であるから左辺各項の$x^m$の係数はすべて零でなければならない。従って

$(m+2)(m+1)c_{m+2}-m(m-1)c_m-2m{c}_m+\nu (\nu +1)c_m=0$

が成り立つ。これを$c_{m+2}$に関して解けば漸化式

• $c_{m+2}=-\frac{(\nu -m)(\nu +m+1)}{(m+2)(m+1)}{c}_m$

が得られる。この漸化式の番号mに0か1を代入してxの冪の係数を順次求めてゆき冪級数

$y=\sum _{m=0}^{\mathrm{\infty }}{c}_m{x}^m$

に代入したら晴れてルジャンドルの微分方程式の無限級数解が得られる。

　また$\nu =n(n\in \mathbb{\mathbb{N}})$の場合基本解が、番号mが偶数の時と奇数の時の2パターンに分かれて現れ、これらの線形結合が上記微分方程式の一般解を構成する事になる。尚、この時或る一方の基本解は無限級数になる事が、もう一方の基本解は多項式になる事が知られている（詳細については専門書を参照されたい）。

ルジャンドルの多項式

ルジャンドルの微分方程式の基本解のうち多項式解の方に注目する。この多項式解の中で$x=1$の時の関数値が1になるよーな関数列をルジャンドル多項式といい$P_n(x),(n=0,1,2,3,\cdots )$で表わす。（これについても詳しい事は当該ページを御覧下さいませ☆）

このルジャンドル多項式には簡潔な表示式が存在する。以下に述べるロドリゲスの公式（またはロドリーグの公式）と呼ばれる或る関数である。

ロドリゲスの公式は次式で与えられる。；

• $P_n(x)=\frac{1}{2^{n}n!}\frac{d^n}{d{x}^n}(x^2-1{)}^n$

本稿ではこの表示式がルジャンドルの微分方程式を満たす事を示すだけにとどめ応用その他についてまでは触れない事にする。

ロドリゲスの公式の証明

まず$u=(x^2-1{)}^n$とゆー関数を考えて、これを微分したら$u^{\prime }=2nx(x^2-1{)}^{n-1}$が成り立つがこの両辺に$x^2-1$掛けたら

$(x^2-1)u^{\prime }=2nxu$

が得られる。で、この等式の両辺をライプニッツの公式

$(f\cdot g{)}^{(n)}=\sum _{k=0}^n{}_n{C}_k{f}^{(k)}{g}^{(n-k)}$

を用いてn+1回微分すれば左辺は

$(x^2-1)u^{(n+2)}+2(n+1)x{u}^{(n+1)}+n(n+1)u^{(n)}$

となり、右辺は

$2n\{x{u}^{(n+1)}+(n+1)u^{(n)}\}$

となる。ここで$y=u^{(n)}$と置けば上式は

$(1-x^2)y^{″}-2x{y}^{\prime }+n(n+1)y=0$

と変形できて上述の関数uの第n階導関数yが確かにルジャンドルの微分方程式を満たしてる事が分かる。

さて、ここで関数

$y(x)=\frac{d^n}{d{x}^n}(x^2-1{)}^n$

にライプニッツの公式を適用すると

$\begin{array}{rl}y(x)& =\frac{d^n}{d{x}^n}\{(x-1{)}^n(x+1{)}^n\}\\ & =\sum _{k=0}^n{}_n{C}_k\{(x+1{)}^n{\}}^{(n-k)}\{(x-1{)}^n{\}}^{(k)}\\ & ={}_n{C}_0\{(x-1{)}^n{\}}^{(n)}\{(x+1{)}^n{\}}^{(0)}+\cdots \end{array}$

となるが、上式の関数yはxに1を代入したら「…」以降の項がすべて零になるので

$y(1)=n!(1+1{)}^n=2^n\cdot n!$

が成り立つ。んで

$P_n(x)=Cy(x),(C=const)$

と置けばルジャンドル多項式の定義より

$P_n(1)=Cy(1)=C\cdot 2^n\cdot n!=1$

であるから「$C=1/2^{n}n!$」が得られる。ゆえに上述の等式よりロドリゲスの公式が導かれた事になる。（証明終）

関連項目

• ガウスの微分方程式
• ルジャンドル多項式