この記事では一般化楕円積分にまつわるルジャンドル関係式を紹介します。
まずルジャンドル関係式とは次の恒等式のことを言うのでした。
完全楕円積分
\begin{align}
K(k)&=\int^1_0\frac{dt}{\sqrt{(1-t^2)(1-k^2t^2)}}\\
E(k)&=\int^1_0\sqrt{\frac{1-k^2t^2}{1-t^2}}dt
\end{align}
および$k'=\sqrt{1-k^2}$に対して
$$K(k)E(k')+K(k')E(k)-K(k)K(k')=\frac\pi2$$
が成り立つ。
そして芝浦工業大学の竹内慎吾教授(日本人!)はこれの一般化を考え、2016年(割と最近!)に
Legendre-type relations for generalized complete elliptic integrals
という論文で以下のように具体的な一般化を与えたようです。
$4$つのパラメーター
$$p\in(-\infty,0)\cup(1,\infty],\quad q,r\in(1,\infty),\quad k\in[0,1)$$
に対して一般化完全楕円積分$K_{p,q,r}(k),E_{p,q,r}(k)$を
\begin{align}
K_{p,q,r}(k)&=\int^1_0\frac{dt}{\sqrt[p]{1-t^q}\sqrt[r]{1-k^qt^q}}\\
E_{p,q,r}(k)&=\int^1_0\frac{\sqrt[r]{1-k^qt^q}}{\sqrt[p]{1-t^q}}dt
\end{align}
と定める。
このとき$K_{2,2,2}(k)=K(k),\ E_{2,2,2}=E(k)$が成り立つことに注意する。
一般化円周率$\pi_{p,q}$を
$$\pi_{p,q}=2\int^1_0\farc{dt}{\sqrt[p]{1-t^q}}$$
と定める。
このとき$\pi_{2,2}=\pi$が成り立つことに注意する。
写像$x\mapsto x^*$を$\frac1x+\frac1{x^*}=1$によって定め、また
$$k'=\sqrt[r]{1-k^q},\quad\frac1s=\frac1p-\farc1q$$
とおくと
$$K_{p,q,r^*}(k)E_{p,r,q}(k')+K_{p,r,q^*}(k')E_{p,q,r}(k)-K_{p,q,r^*}(k)K_{p,r,q^*}(k')=\frac{\pi_{p,q}\pi_{s,r}}{4}$$
が成り立つ(ただし$\frac1\infty=0$とする)。
実際に$p=q=r=2$のときを考えてみると$q^*=r^*=2,\;s=\infty$であり
$$\pi_{\infty,q}=2\int^1_0\frac{dt}{(1-t^q)^{\frac1\infty}}=2\int^1_0dt=2$$
に注意すると最初に挙げたルジャンドル関係式
$$K(k)E(k')+K(k')E(k)-K(k)K(k')=\frac\pi2$$
が出てくることがわかります。
ちなみにこの一般化ルジャンドル関係式は超幾何関数にまつわるElliottの恒等式
\begin{eqnarray}
&&\F{\frac12+a}{-\farc12-b}{a+c+1}{z}\F{\frac12-a}{\frac12+b}{b+c+1}{1-z}
\\&&\quad+\F{\frac12+a}{\frac12-b}{a+c+1}{z}\F{-\frac12-a}{\frac12+b}{b+c+1}{1-z}
\\&&\quad-\F{\frac12+a}{\frac12-b}{a+c+1}{z}\F{\frac12-a}{\frac12+b}{b+c+1}{1-z}
\\&=&\frac{\G(a+c+1)\G(b+c+1)}{\G(a+b+c+\frac32)\G(c+\frac12)}
\end{eqnarray}
と等価であることが知られています(このことについては
別の記事
で詳しく解説します)。
E.B.Elliottがこの恒等式を超幾何関数の観点から証明したところを、竹内教授は楕円積分の観点からの別証明を与えたという構図になっているようです。
まず$K_{p,q,r}$と$E_{p,q,r}$との間に
\begin{align}
\frac{dK_{p,q,r}}{dk}&=\frac{a'E_{p,q,r^*}-(a'-k^q)K_{p,q,r}}{k(1-k^q)}
\quad\l(a'=1+\frac q{r^*}-\frac qp\r)\\
\frac{dE_{p,q,r}}{dk}&=\frac{q(E_{p,q,r}-K_{p,q,r^*})}{rk}
\end{align}
という微分方程式が成り立つことを示す。
そして
$$K_{p,q,r}'=K_{p,q,r}(k'),\quad E_{p,q,r}'=E_{p,q,r}(k')$$
とおいたとき、上の微分方程式から
$$L=K_{p,q,r^*}E_{p,r,q}'+K_{p,r,q^*}'E_{p,q,r}-K_{p,q,r^*}K_{p,r,q^*}'$$
が
$$\frac{dL}{dk}=0$$
を満たすことを確かめる。
したがって
$$L(k)=\mathrm{Const.}$$
がわかるので$k\to0^+$において
$$K_{p,q,r^*}(0)=\farc{\pi_{p,q}}{2},\quad
E_{p,r,q}(1)=\frac{\pi_{s,r}}{2}$$
および
$$\lim_{k\to0^+}(E_{p,q,r}-K_{p,q,r^*})K_{p,r,q^*}'=0$$
と計算できることから
$$L=\lim_{k\to0^+}L(k)=\frac{\pi_{p,q}\pi_{s,r}}{4}$$
を得る。といった具合になります。
まず以下での議論を簡単にするために一般化した三角関数を導入します。
一般化正弦関数$\sin_{p,q}:[0,\farc{\pi_{p,q}}2]\to[0,1]$を関数
$$f_{p,q}(x)=\int^x_0\frac{dt}{\sqrt[p]{1-t^q}}\quad(x\in[0,1])$$
の逆関数として定義する。
また$p\neq\infty$において一般化余弦関数$\cos_{p,q}:\left(0,\farc{\pi_{p,q}}2\right)\to[0,1]$を
$$\cos_{p,q}x=(\sin_{p,q}x)'$$
と定める。
いま$\sin_{p,q}=f_{p,q}^{-1}$に注意すると逆関数の微分法から
$$\cos_{p,q}x=(f_{p,q}^{-1}(x))'=\frac{1}{f'_{p,q}(f_{p,q}^{-1}(x))}=\sqrt[p]{1-\sin_{p,q}^qx}$$
つまり
$$\cos^p_{p,q}x+\sin_{p,q}^qx=1$$
がわかります。
ただ$p=\infty$においては$\cos_{p,q}$がうまく定義されない(性質が崩れる)ため代わりの関数を
$$\cos_{*p,q}x=1-\sin^q_{p,q}x$$
と定めておきます。このとき$p=\infty$においても
\begin{align}
(\sin_{p,q}x)'&=\cos_{*p,q}^{\frac1p}x\\
(\cos_{*p,q}x)'&=-q\sin_{p,q}^{q-1}x\cos_{*p,q}^{\frac1p}x
\end{align}
が成り立ちます。
ここで変数変換$t=\sin_{p,q}\theta$を考えると$\t=f_{p,q}(t)$から
$$d\t=\frac{dt}{\sqrt[p]{1-t^q}}$$
に注意すると
\begin{align}
K_{p,q,r}(k)&=\int^{\frac12\pi_{p,q}}_0\farc{d\t}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}\\
E_{p,q,r}(k)&=\int^{\frac12\pi_{p,q}}_0\sqrt[r]{1-k^q\sin_{p,q}^q\t}d\t
\end{align}
と表せることが肝になってきます。
次に上で予告したように次の微分方程式を示します。
$$a'=1+\frac q{r^*}-\frac qp$$
とおくと
\begin{align}
\frac{dK_{p,q,r}}{dk}&=\frac{a'E_{p,q,r^*}-(a'-k^q)K_{p,q,r}}{k(1-k^q)}\\
\frac{dE_{p,q,r}}{dk}&=\frac{q(E_{p,q,r}-K_{p,q,r^*})}{rk}
\end{align}
が成り立つ。
$E_{p,q,r}$については簡単にわかる。
\begin{eqnarray}
\frac{dE_{p,q,r}}{dk}
&=&\frac qr\int^1_0\farc{-k^{q-1}t^q}{\sqrt[p]{1-t^q}(1-k^qt^q)^{1-\farc1r}}dt
\\&=&\frac q{rk}\int^1_0\farc{(1-k^qt^q)-1}{\sqrt[p]{1-t^q}(1-k^qt^q)^{1-\farc1r}}dt
\\&=&\frac q{rk}\l(\int^1_0\farc{\sqrt[r]{1-k^qt^q}}{\sqrt[p]{1-t^q}}dt
-\int^1_0\farc{dt}{\sqrt[p]{1-t^q}(1-k^qt^q)^{1-\farc1r}}\r)
\\&=&\frac q{rk}(E_{p,q,r}-K_{p,q,r^*})
\end{eqnarray}
$K_{p,q,r}$については以下の補題を用いて示す。
$\dis\farc{d}{d\t}\sqrt[r]{\frac{\cos_{*p,q}\t}{1-k^q\sin^q_{p,q}\t}} =\farc{q(k^q-1)}{r}\frac{\sin^{q-1}_{p,q}\t\cos_{*p,q}^{\farc1p-\farc1{r^*}}\t}{(1-k^q\sin^q_{p,q}\t)^{1+\frac1r}}$
$\farc1r+\frac1{r^*}=1$および$\cos_{*p,q}\t=1-\sin^q_{p,q}\t$としていたことに注意すると
\begin{eqnarray}
\farc{d}{d\t}\sqrt[r]{\frac{\cos_{*p,q}\t}{1-k^q\sin^q_{p,q}\t}}
&=&\farc1r\left(\frac{\cos_{*p,q}\t}{1-k^q\sin^q_{p,q}\t}\right)^{\farc1r-1}\cdot
\farc{d}{d\t}\left(\frac{\cos_{*p,q}\t}{1-k^q\sin^q_{p,q}\t}\right)
\\&=&\frac1r\farc{\cos_{*p,q}^{-\frac1{r^*}}\t}{(1-k^q\sin_{p,q}^q\t)^{\farc1r-1}}\cdot
\farc{(1-k^q\sin^q_{p,q}\t)(1-\sin_{p,q}^q\t)'-(1-k^q\sin^q_{p,q}\t)'\cos_{*p,q}\t}{(1-k^q\sin^q_{p,q}\t)^2}
\\&=&\frac1r\farc{\cos_{*p,q}^{-\frac1{r^*}}\t(\sin^q_{p,q}\t)'(-(1-k^q\sin_{p,q}^q\t)+k^q\cos_{*p,q}\t)}{(1-k^q\sin_{p,q}^q\t)^{1+\farc1r}}
\\&=&\frac1r\farc{\cos_{*p,q}^{-\frac1{r^*}}\t(q\sin_{p,q}^{q-1}\t\cos_{*p,q}^{\frac1p}\t)(k^q-1)}{(1-k^q\sin_{p,q}^q\t)^{1+\farc1r}}
\\&=&\farc{q(k^q-1)}{r}\frac{\sin^{q-1}_{p,q}\t\cos_{*p,q}^{\farc1p-\farc1{r^*}}\t}{(1-k^q\sin^q_{p,q}\t)^{1+\frac1r}}
\end{eqnarray}
とわかる。
いま定義より
$$\sin_{p,q}0=0,\quad
\sin_{p,q}\frac{\pi_{p,q}}{2}=1,\quad
\cos_{*p,q}\frac{\pi_{p,q}}{2}=0$$
が成り立つことに注意すると
\begin{eqnarray}
\farc{dK_{p,q,r}}{dk}
&=&\farc qr\int^{\frac12\pi_{p,q}}_0\farc{k^{q-1}\sin^q\t}{(1-k^q\sin_{p,q}^q\t)^{1+\farc1r}}d\t
\\&=&\farc{k^{q-1}}{k^q-1}\int^{\frac12\pi_{p,q}}_0\left(\farc{q(k^q-1)}{r}\frac{\sin^{q-1}_{p,q}\t\cos_{*p,q}^{\farc1p-\farc1{r^*}}\t}{(1-k^q\sin^q_{p,q}\t)^{1+\frac1r}}\right)\sin_{p,q}\t\cos_{*p,q}^{\frac1{r^*}-\frac1p}\t\;d\t
\\&=&-\frac{k^{q-1}}{k^q-1}\int^{\frac12\pi_{p,q}}_0\frac{\cos^{\frac1r}_{*p,q}\t}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}
\cdot\frac{d}{d\t}(\sin_{p,q}\t\cos_{*p,q}^{\frac1{r^*}-\frac1p}\t)d\t
\\&=&\frac{k^{q-1}}{1-k^q}\int^{\frac12\pi_{p,q}}_0\frac{\cos^{\frac1r}_{*p,q}\t}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}(\cos_{*p,q}^{\frac1{r^*}}\t-q\l(\frac1{r^*}-\frac1p\r)\sin_{p,q}^q\t\cos^{\frac1{r^*}-1}\t)d\t
\\&=&\frac{k^{q-1}}{1-k^q}\int^{\frac12\pi_{p,q}}_0\frac{\cos_{*p,q}\t-(a'-1)\sin_{p,q}^q\t}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}d\t
\\&=&\frac{1}{k(1-k^q)}\int^{\frac12\pi_{p,q}}_0\frac{k^q(1-a'\sin_{p,q}^q\t)}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}d\t
\\&=&\frac1{k(1-k^q)}\int^{\frac12\pi_{p,q}}_0\frac{k^q-a'}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}d\t+\frac1{k(1-k^a)}\int^{\frac12\pi_{p,q}}_0\frac{a'(1-k^q\sin_{p,q}^q\t)}{\sqrt[r]{1-k^q\sin_{p,q}^q\t}}d\t
\\&=&\frac{k^q-a'}{k(1-k^q)}K_{p,q,r}+\frac{a'}{k(1-k^a)}E_{p,q,r^*}
\end{eqnarray}
を得る。
$$k'=\sqrt[r]{1-k^q},\quad
K_{p,q,r}'=K_{p,q,r}(k'),\quad
E_{p,q,r}'=E_{p,q,r}(k')$$
および
$$a=1+\farc qr-\frac qp,\quad b=1+\frac rq-\frac rp$$
とおくと
\begin{align}
\farc{dK_{p,q,r^*}}{dk}&=\frac{aE_{p,q,r}-(a-k^q)K_{p,q,r^*}}{k(k')^r}&
\frac{dE_{p,q,r}}{dk}&=\frac{q(E_{p,q,r}-K_{p,q,r^*})}{rk}\\
\farc{dK'_{p,r,q^*}}{dk}&=-\frac{q(bE'_{p,r,q}-(b-(k')^r)K'_{p,r,q^*})}{rk(k')^r}&
\frac{dE_{p,r,q}}{dk}&=-\frac{k^{q-1}(E'_{p,r,q}-K'_{p,r,q^*})}{(k')^r}
\end{align}
が成り立つ。
$$\frac{dk'}{dk}=-\frac{q}{r}\frac{k^{q-1}}{(k')^{r-1}}$$
および$(x^*)^*=x$に注意するとわかる。
$$L=K_{p,q,r^*}E_{p,r,q}'+K_{p,r,q^*}'E_{p,q,r}-K_{p,q,r^*}K_{p,r,q^*}'$$
とおくと
$$\farc{dL}{dk}=0$$
が成り立つ。
定理4系に注意して計算していくと
\begin{eqnarray}
\farc{dL}{dk}
&=&\frac{aE_{p,q,r}-(a-k^q)K_{p,q,r^*}}{k(k')^r}E_{p,r,q}'-K_{p,q,r^*}\frac{k^{q-1}(E'_{p,r,q}-K'_{p,r,q^*})}{(k')^r}
\\&&\quad-\frac{q(bE'_{p,r,q}-(b-(k')^r)K'_{p,r,q^*})}{rk(k')^r}E_{p,q,r}+K'_{p,r,q^*}\frac{q(E_{p,q,r}-K_{p,q,r^*})}{rk}
\\&&\quad-\frac{aE_{p,q,r}-(a-k^q)K_{p,q,r^*}}{k(k')^r}K'_{p,r,q^*}+K_{p,q,r^*}\frac{q(bE'_{p,r,q}-(b-(k')^r)K'_{p,r,q^*})}{rk(k')^r}
\\&=&\left(-\frac{a-k^q}{k(k')^r}-\frac{k^{q-1}}{(k')^r}+\frac{qb}{rk(k')^r}\right)K_{p,q,r^*}E'_{p,r,q}
\\&&\quad+\left(\frac{q(b-(k')^r)}{rk(k')^r}+\frac{q}{rk}-\frac{a}{k(k')^r}\right)K'_{p,r,q^*}E_{p,q,r}
\\&&\quad+\left(\frac{k^{q-1}}{(k')^r}-\farc{q}{rk}+\frac{a-k^q}{k(k')^r}-\farc{q(b-(k')^r)}{rk(k')^r}\right)K_{p,q,r^*}K'_{p,r,q^*}
\\&&\quad+\left(\frac{a}{k(k')^r}-\farc{qb}{rk(k')^r}\right)E_{p,q,r}E'_{p,r,q}
\\&=&\frac{qb-ra}{rk(k')^r}(K_{p,q,r^*}E'_{p,r,q}+K'_{p,r,q^*}E_{p,q,r}-K_{p,q,r^*}K'_{p,r,q^*}-E_{p,q,r}E'_{p,r,q})
\end{eqnarray}
と整理でき、$a,b$の取り方から$qb-ra=0$が成り立つことに注意すると主張を得る。
$$\farc1v=\farc1p-\frac1r$$
とおくと
$$K_{p,q,r}(0)=\frac{\pi_{p,q}}{2},\quad E_{p,q,r}(1)=\frac{\pi_{v,q}}{2}$$
が成り立つ。
$$K_{p,q,r}(0)=\int^1_0\farc{dt}{\sqrt[p]{1-t^q}},\quad
E_{p,q,r}(1)=\int^1_0\frac{dt}{(1-t^q)^{\farc1p-\frac1r}}$$
および一般化円周率の定義から明らか。
$\dis\lim_{k\to0^+}(E_{p,q,r}-K_{p,q,r^*})K_{p,r,q^*}'=0$
\begin{eqnarray}
&&(K_{p,q,r^*}-E_{p,q,r})K_{p,r,q^*}'
\\&=&\int^{\farc12\pi_{p,q}}_0\Big(\farc{1}{(1-k^q\sin_{p,q}^q\t)^{\frac1{r^*}}}-(1-k^q\sin^q_{p,q}\t)^{\farc1r}\Big)d\t\cdot\int^{\frac12\pi_{p,r}}_0\frac{d\t}{(1-(k')^r\sin^r_{p,r}\t)^{\farc1{q^*}}}
\\&=&\int^{\farc12\pi_{p,q}}_0\farc{1-(1-k^q\sin^q\t)}{(1-k^q\sin_{p,q}^q\t)^{\frac1{r^*}}}d\t\cdot\int^{\frac12\pi_{p,r}}_0\frac{d\t}{(1-(1-k^q)\sin^r_{p,r}\t)^{\farc1{q^*}}}
\\&\leq&\int^{\farc12\pi_{p,q}}_0\farc{k^q1^q}{(1-k^q\sin_{p,q}^q\t)^{\frac1{r^*}}}d\t\cdot\int^{\frac12\pi_{p,r}}_0\frac{d\t}{(k^q(1-\sin^r_{p,r}\t)+k^q\sin^r_{p,r}\t)^{\farc{q-1}{q}}}
\\&=&k^qK_{p,q,r^*}(k)\cdot\farc{1}{k^{q-1}}\farc{\pi_{p,r}}{2}
\\&=&\frac{\pi_{p,r}}{2}kF_{p,q,r^*}(k)\to0\quad(as\;k\to0^+)
\end{eqnarray}
とわかる。
以上より
$$\frac1s=\frac1p-\frac1q$$
とおくと
$$L=\lim_{k\to0^+}L(k)=K_{p,q,r^*}(0)E_{p,r,q}(1)=\frac{\pi_{p,q}\pi_{s,r}}{4}$$
を得る。