数理解析研究所講究録
第 1692 巻 2010 年 88-100

88

ガウシアンを超えて横浜市立大学・国際総合科学部 藤井一幸 (FUJII Kazuyuki)
International College of Arts and Sciences
Yokohama City University

要約(Gaussian)とは、 ガウス関数 (分布)$e^{-(px^{2}+qx+r)}$に絡むすべての対象の略称である。例えばガウシアンビーム、 ガウシアンプロセス、 ガウシアンノイズのガウシアン

ように。現代の基礎及び応用科学の支柱の一つである。
我々の目的はガウシアンを乗り越えて、新しいパラダイムを構築することである。 これは誰かが、 どこかでなさなければならない基礎研究である。

$[I|$

はじめに

(Gaussian)とは、 ガウシアンビーム、 ガウシアンプロセス、 ガウシア
ンノイズのようにガウス関数 (分布)ガウシアン$e^{-f(x)}$;$f(x)=px^{2}+qx+r$(1)に絡むすべての対象の略称である。
ガウス関数は 2 変数に広げると$e^{-f(x,y)}$;$f(x, y)=px^{2}+qxy+ry^{2}=(x, y)(pq2$$rq2)(\begin{array}{l}xy\end{array})\equiv x^{t}Ax$(2)$(|A|>0)$ ,(3)と 2 次形式で表せる。 これの特徴は 2 重積分$\int\int_{R^{2}}e^{-(px^{2}+qxy+ry^{2})}dxdy=\frac{\pi}{\sqrt{|A|}}$$|A|=pr- \frac{q^{2}}{4}=-\frac{1}{4}(q^{2}-4pr)=-\frac{1}{4}D$で与えられる。 $|A|$ は $f(x, y)$ の係数から決まる行列式で、 $D$ は 2 次方程式 $px^{2}+qx+r=0$の判別式である。
この関数の強みは、 関数 (積分) を多次元化出来ることである。即ち、 $x$ を縦ベクト
ル、 $A$ を対応する対称行列として$f(x)=x^{t}Ax$ ;$x=(x_{1}, x_{2}, \cdots, x_{n})^{t}$,$A^{t}=A$

89と表せ、 従って公式$\int\int\cdots\int_{R^{n}}e^{-x^{t}Ax}dx_{1}dx_{2}\cdots dx_{n}=\frac{\pi^{\frac{n}{2}}}{\sqrt{|A|}}$$(|A|>0)$を得るわけである。確かに使いやすい公式で、 広範囲に使われている所以である。
では、 ガウス関数を$e^{-f(x,y)}$,$f(x, y)=ax^{3}+bx^{2}y+cxy^{2}+dy^{3}$(4)と 3 次形式に拡張したらガウシアンとどのような違いが出て来るのであろうか ? この系はもちろんはノンガウシアンで、最初に考えるガウシアンの拡張としては妥当なものである。

更にガウス関数を$e^{-f(x,y)}$

と一般の

$n$

,

$f(x, y)=a_{0}x^{n}+a_{1}x^{n-1}y+\cdots+a_{j}x^{n-j}y^{j}+\cdots+a_{n-1}xy^{n-1}+a_{n}y^{n}$(5)

次形式に拡張したらどのような違いが出て来るのであろうか ?

(4), (5) に対して (3) のような積分公式は成り立つのか ? もし成り立てば (3) の行列式
や判別式に対応するものは何か ? 更に多次元化することが出来るのか ? これは誰かが、どこかでなさなければならない基礎研究である。

[II] ノンガウシアンに向けて
この新しい系 (4) に対して積分公式 (3) に対応するものは何か ? 残念ながら$\int\int_{R^{2}}e^{-(ax^{3}+bx^{2}y+cxy^{2}+dy^{3})}dxdy=\infty$(6)と発散してしまう。 これを打開するためロシアの Morozov &Shakirov は論文 [1] で、空ではなくと同相なある曲面 $X$ 上で積分しての中の間を$R^{2}$$C^{2}$$R^{2}$$\int\int_{x^{e^{-(ax^{3}+bx^{2}y+cxy^{2}+dy^{3}}}})_{dxdy=}\frac{C}{\sqrt{}^{6}\overline{D}}$(7)

なる “公式” を与えた。 ここに $C$ は定数で、 $D$ は 3 次方程式

$ax^{3}+bx^{2}+cx+d=0$の判別式$D=b^{2}c^{2}+18abcd-4ac^{3}-4b^{3}d-27a^{2}d^{2}$である。 これはある種の繰り込みであるが、 やはり意味不明である。

(8)

90

これを打開するため Fujii は論文 [2] で変数変換 $x=t\rho,$ $y=\rho$ を行って

(6) の左辺$=$$\int\int e^{-\rho^{3}}(at^{3}+bt^{2}+ct+d)|\rho|dtd\rho=\int\{\int e^{-(at^{3}+bt^{2}+ct+d)}\rho^{3}|\rho|d\rho\}dt$$=$$\int|\sigma|e^{-\sigma^{3}}d\sigma\int\frac{1}{|\sqrt[3]{at^{3}+bt^{2}+ct+d}|\sqrt[3]{at^{3}+bt^{2}+ct+d}}dt$とした (最後に変数変換$\sigma=\vee\pi 3+lbt^{2}\mp\sim+l\rho$を使った) 。従って発散するパートは$\int_{R}|\sigma|e^{-\sigma^{3}}d\sigma$

で、 主要パートは(9)$\int_{R}\frac{1}{|\sqrt[3]{ax^{3}+bx^{2}+cx+d}|\sqrt[3]{ax^{3}+bx^{2}+cx+d}}dx$

(

に変えた)

この積分を (6) の繰り込みに採用したいのだが、 その形はバランスが悪く、 更に計算することが非常に難しい (以下の節を見よ) 。そこで “少し” ばか
り修正して

である

$t$

を

$x$

。

定義 (6) の “繰り込み” を$\ddagger\int\int_{R^{2}}e^{-(ax^{3}+bx^{2}y+cxy^{2}+dy^{3})}dxdy\ddagger=\int_{R}\frac{1}{\sqrt[3]{(ax^{3}+bx^{2}+cx+d)^{2}}}dx$(10)と定義する。

Fujii は上の論文でこの積分に関して(A) For $D<0$$\int_{R}\frac{1}{\sqrt[3]{(ax^{3}+bx^{2}+cx+d)^{2}}}dx=\frac{C_{-}}{\sqrt[6]{-D}}$,$C_{-}= \sqrt[3]{2}B(\frac{1}{2}, \frac{1}{6})$,$C_{+}=3B( \frac{1}{3}, \frac{1}{3})$

.

(11)

(B) For $D>0$$\int_{R}\frac{1}{\sqrt[3]{(ax^{3}+bx^{2}+cx+d)^{2}}}dx=\frac{C_{+}}{\sqrt[6]{D}}$(C)$C_{-}$と$c_{+}$

.

(12)

の関係
$C_{+}=\sqrt{}\overline{3}C_{-}$

なる公式を得た (証明は付録を見よ)

。

.

(13)

ここに $B$ はベータ関数$B(p, q)= \int_{0}^{1}x^{p-1}(1-x)^{q-1}dx=\int_{0}^{\infty}\frac{x^{p-1}}{(1+x)^{p+q}}dx$

(14)

91

である。
コメント

3 次方程式の判別式が積分計算に顔を出したことは (Fujii の知る限りでは)無い。 これが面白いところである。

多くの問題が出てくる。 例えば、

1 この公式からガウシアンの場合に相当する種々の結果は出てくるのか ?
$)$2 Morozov &Shakirov の結果と Fujii のそれはどのように関係しているのか ?
$)$3 多次元化は可能なのか ?
$)$1

ガウシアンの場合$)$$f(x, y)=px^{2}+qxy+ry^{2}=(x, y)(pq2$$rq2)(\begin{array}{l}xy\end{array})=x^{t}Ax$と行列表現が出来、 結果として$|A|=pr- \frac{q^{2}}{4}=-\frac{1}{4}(q^{2}-4pr)=-\frac{1}{4}D$より、 行列式と 2 次方程式の判別式が比例している。

しかしノンガウシアンの場合、 例えば$f(x, y)=ax^{3}+bx^{2}y+cxy^{2}+dy^{3}$を行列で表現することは出来ない。 出来るとすれば立体行列の理論が必要になる。他方(11), (12) より積分計算には 3 次方程式 $ax^{3}+bx^{2}+cx+d=0$ の判別式
$D=b^{2}c^{2}+18abcd-4ac^{3}-4b^{3}d-27a^{2}d^{2}$が出てくる。
立体行列の “行列式” と代数方程式の判別式の関係は、 昔ケーリーが研究したと聞いている (詳細は知らない)
。

2

$)$

本質的問題は、 Fujii の計算$\int_{R}\frac{1}{\sqrt[3]{}\overline{(ax^{3}+bx^{2}+cx+d)^{2}}}dx=\frac{C_{-}}{\sqrt{}^{6}\overline{-D}}$$\int_{R}\frac{1}{\sqrt[3]{}\overline{(ax^{3}+bx^{2}+cx+d)^{2}}}dx=\frac{c_{\perp}}{\sqrt{}^{6}\overline{D}}$が示すように $D=0$ を跨いで積分の定数
る$\circ$$(C_{+}$と$C_{-})$

for

$D<0$

for

$D>0$

の値が異なってしまうことであ

このことは Morozov&Shakirov の結果$\int\int_{x^{e^{-(ax^{3}+bx^{2}y+cxy^{2}+dy^{3}}}})_{dxdy=}\frac{C}{\sqrt{}^{6}\overline{D}}$92とconipatible なのであろうか ? と言う疑問である。

3$)$

更に、 ノンガウシアンの多次元化を考えないと多くの応用は期待出来ないかも知れない。 ガウシアンの場合はこれが可能だったからこそ応用範囲が広いのである。今のところ多次元化についてはサッパリわからない。

[III] おわりにもし以上のことがうまく行けば (多分ありえないが)いたい。即ち、、更なるノンガウシアンを取扱$f(x, y)=a_{0}x^{4}+a_{1}x^{3}y+a_{2}x^{2}y^{2}+a_{3}xy^{3}+a_{4}y^{4}$に対して$\int\int_{R^{2}}e^{-f(x,y)}dxdy=\int\int_{R^{2}}e^{-(a_{O}x^{4}+a_{1}x^{3}y+a_{2}x^{2}y^{2}+a_{3}xy^{3}+a_{4}y^{4})}dxdy=\frac{c_{\pm}}{1\sqrt[2]{\pm D}}$が期待されている。 ここに $D$ は 4 次方程式
$a_{0}x^{4}+a_{1}x^{3}+a_{2}x^{2}+a_{3}x+a_{4}=0$の判別式$D_{n=4}$$=$$256a_{0}^{3}a_{4}^{3}-4a_{1}^{3}a_{3}^{3}-27a_{0}^{2}a_{3}^{4}-27a_{1}^{4}a_{4}^{2}-128a_{0}^{2}a_{2}^{2}a_{4}^{2}+a_{1}^{2}a_{2}^{2}a_{3}^{2}+16a_{0}a_{2}^{4}a_{4}$$-4a_{0}a_{2}^{3}a_{3}^{2}-4a_{1}^{2}a_{2}^{3}a_{4}+144a_{0}^{2}a_{2}a_{3}^{2}a_{4}-6a_{0}a_{1}^{2}a_{3}^{2}a_{4}+144a_{0}a_{1}^{2}a_{2}a_{4}^{2}-192a_{0}^{2}a_{1}a_{3}a_{4}^{2}$$+18a_{0}a_{1}a_{2}a_{3}^{3}+18a_{1}^{3}a_{2}a_{3}a_{4}-80a_{0}a_{1}a_{2}^{2}a_{3}a_{4}$

,

である。
更に
$f(x, y)=a_{0}x^{5}+a_{1}x^{4}y+a_{2}x^{3}y^{2}+a_{3}x^{2}y^{3}+a_{4}xy^{4}+a_{5}y^{5}$ .

に対しては$\int\int_{R^{2}}e^{-(+3}a_{0}x^{5}+a_{1}x^{4}y+a_{2}x^{3}y^{2}ax^{2}y^{3}+a_{4}xy^{4}+a_{5}y^{5})dxdy=\infty$なので、 “繰り込み”$\ddagger\int\int_{R^{2}}e^{-f(x,y)}dxdy\ddagger=\int_{R}\frac{1}{\sqrt[5]{(a_{0}x^{5}+a_{1}x^{4}+a_{2}x^{3}+a_{3}x^{2}+a_{4}x+a_{5})^{2}}}dx$

をしなければならない。 このとき$\ddagger\int\int_{R^{2}}e^{-f(x,y)}dxdy\ddagger=\frac{c_{\pm}}{2\sqrt[0]{D\pm}}$は期待されるか ?

93

ここに

$D$

は 5 次方程式$a_{0}x^{5}+a_{1}x^{4}+a_{2}x^{3}+a_{3}x^{2}+a_{4}x+a_{5}=0$の判別式$D_{n=5}$$=$$3125a_{0}^{4}a_{5}^{4}-2500a_{0}^{3}a_{1}a_{4}a_{5}^{3}-3750a_{0}^{3}a_{2}a_{3}a_{5}^{3}+2000a_{0}^{3}a_{2}a_{4}^{2}a_{5}^{2}+2250a_{0}^{3}a_{3}^{2}a_{4}a_{5}^{2}$$-1600a_{0}^{3}a_{3}a_{4}^{3}a_{5}+256a_{0}^{3}a_{4}^{5}+2000a_{0}^{2}a_{1}^{2}a_{3}a_{5}^{3}-50a_{0}^{2}a_{1}^{2}a_{4}^{2}a_{5}^{2}+2250a_{0}^{2}a_{1}a_{2}^{2}a_{5}^{3}$$-2050a_{0}^{2}a_{1}a_{2}a_{3}a_{4}a_{5}^{2}+160a_{0}^{2}a_{1}a_{2}a_{4}^{3}a_{5}-900a_{0}^{2}a_{1}a_{3}^{3}a_{5}^{2}+1020a_{0}^{2}a_{1}a_{3}^{2}a_{4}^{2}a_{5}-192a_{0}^{2}a_{1}a_{3}a_{4}^{4}$$-900a_{0}^{2}a_{2}^{3}a_{4}a_{5}^{2}+825a_{0}^{2}a_{2}^{2}a_{3}^{2}a_{5}^{2}+560a_{0}^{2}a_{2}^{2}a_{3}a_{4}^{2}a_{5}-128a_{0}^{2}a_{2}^{2}a_{4}^{4}-630a_{0}^{2}a_{2}a_{3}^{3}a_{4}a_{5}$$+144a_{0}^{2}a_{2}a_{3}^{2}a_{4}^{3}+108a_{0}^{2}a_{3}^{5}a_{5}-27a_{0}^{2}a_{3}^{4}a_{4}^{2}-1600a_{0}a_{1}^{3}a_{2}a_{5}^{3}+160a_{0}a_{1}^{3}a_{3}a_{4}a_{5}^{2}$$-36a_{0}a_{1}^{3}a_{4}^{3}a_{5}+1020a_{0}a_{1}^{2}a_{2}^{2}a_{4}a_{5}^{2}+560a_{0}a_{1}^{2}a_{2}a_{3}^{2}a_{5}^{2}-746a_{0}a_{1}^{2}a_{2}a_{3}a_{4}^{2}a_{5}+144a_{0}a_{1}^{2}a_{2}a_{4}^{4}$$+24a_{0}a_{1}^{2}a_{3}^{3}a_{4}a_{5}-6a_{0}a_{1}^{2}a_{3}^{2}a_{4}^{3}-630a_{0}a_{1}a_{2}^{3}a_{3}a_{5}^{2}+24a_{0}a_{1}a_{2}^{3}a_{4}^{2}a_{5}+356a_{0}a_{1}a_{2}^{2}a_{3}^{2}a_{4}a_{5}$$-80a_{0}a_{1}a_{2}^{2}a_{3}a_{4}^{3}-72a_{0}a_{1}a_{2}a_{3}^{4}a_{5}+18a_{0}a_{1}a_{2}a_{3}^{3}a_{4}^{2}+108a_{0}a_{2}^{5}a_{5}^{2}-72a_{0}a_{2}^{4}a_{3}a_{4}a_{5}+16a_{0}a_{2}^{4}a_{4}^{3}$$+16a_{0}a_{2}^{3}a_{3}^{3}a_{5}-4a_{0}a_{2}^{3}a_{3}^{2}a_{4}^{2}+256a_{1}^{5}a_{5}^{3}-192a_{1}^{4}a_{2}a_{4}a_{5}^{2}-128a_{1}^{4}a_{3}^{2}a_{5}^{2}+144a_{1}^{4}a_{3}a_{4}^{2}a_{5}$$-27a_{1}^{4}a_{4}^{4}+144a_{1}^{3}a_{2}^{2}a_{3}a_{5}^{2}-6a_{1}^{3}a_{2}^{2}a_{4}^{2}a_{5}-80a_{1}^{3}a_{2}a_{3}^{2}a_{4}a_{5}+18a_{1}^{3}a_{2}a_{3}a_{4}^{3}+16a_{1}^{3}a_{3}^{4}a_{5}$$-4a_{1}^{3}a_{3}^{3}a_{4}^{2}-27a_{1}^{2}a_{2}^{4}a_{5}^{2}+18a_{1}^{2}a_{2}^{3}a_{3}a_{4}a_{5}-4a_{1}^{2}a_{2}^{3}a_{4}^{3}-4a_{1}^{2}a_{2}^{2}a_{3}^{3}a_{5}+a_{1}^{2}a_{2}^{2}a_{3}^{2}a_{4}^{2}$

.

である。

次第に凄まじい様相を呈してきて、 コントロール出来るか否かよくわからない。
我々の研究目的はガウシアンを乗り越えて、 新しいパラダイムを構築することである。
これは誰かが、 どこかでなさなければならない基礎研究である。

幼少の伊豆の踊り子達にとって天城越えは辛かったように、我々数学者や数理物理学者
にとってもガウス越えは本当にハードなのである。

付録: 公式の証明

この付録で公式 (A), (B), (C) (各々 (11), (12), (13)) の証明を与える。

$($

A の証明
$)$

: $D<0$ の場合に (11) を証明する。

まず 3 次の多項式 $ax^{3}+bx^{2}+cx+d$ で、 $a=0$ の場合を考える$\int_{R}\frac{1}{\sqrt[3]{}\overline{(bx^{2}+cx+d)^{2}}}dx$

$0$

即ち、

94

を計算する。 この場合一 $D=b^{2}(4bd-c^{2})>0$ に注意すると$\int_{R}\frac{1}{\sqrt[3]{(bx^{2}+cx+d)^{2}}}dx$$=$$\int_{R}\frac{1}{\sqrt[3]{b^{2}}\sqrt[3]{(x^{2}+\frac{c}{b}x+\frac{d}{b})^{2}}}dx$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{R}\frac{1}{\sqrt[3]{((x+\frac{c}{2b})^{2}+\frac{d}{b}-\frac{c^{2}}{4b^{2}})^{2}}}dx$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{R}\frac{1}{\sqrt[3]{(x^{2}+\frac{4bd-c^{2}}{4b^{2}})^{2}}}dx$$T^{2} \equiv\frac{4bd-c^{2}}{4b^{2}}>0$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{R}\frac{1}{\sqrt[3]{(x^{2}+T^{2})^{2}}}dx$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{R}\frac{T}{\sqrt[3]{T^{4}}\sqrt[3]{(y^{2}+1)^{2}}}dy\Leftarrow x=Ty$$=$$\frac{2}{\sqrt[3]{b^{2}T}}\int_{0}^{\infty}\frac{1}{\sqrt[3]{(y^{2}+1)^{2}}}dy$$=$$\frac{2}{\sqrt[3]{b^{2}T}}\int_{0}^{\infty}\frac{1}{\sqrt[3]{(x+1)^{2}}}\frac{dx}{2\sqrt{x}}\Leftarrow y=\sqrt{x}$$=$$\frac{1}{\sqrt[3]{b^{2}T}}\int_{0}^{\infty}\frac{x^{-\frac{1}{2}}}{(x+1)^{\frac{2}{3}}}dx$

$=$

$=$

$=$

$=$

$\frac{B(\frac{1}{2},\frac{1}{\tau^{6}})}{\sqrt[3]{b^{2}}}\Leftarrow$(14)$\frac{B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{b^{4}T^{2}}}$$\frac{\sqrt[3]{2}B(\tilde{2}1\frac{1}{6})}{\sqrt[6]{b^{2}(4bd-c^{2})}}\Leftarrow T^{2}\equiv\frac{4bd-c^{2}}{4b^{2}}$$\frac{\sqrt[3]{2}B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{-D}}$

(15)

を得る。

条件 $D<0$ より $ax^{3}+bx^{2}+cx+d=0$ は一つの実根をもつ。 最初にそれを $x=0$ とする。 この場合 $d=0$ で$ax^{3}+bx^{2}+cx+d=x(ax^{2}+bx+c)$となる。 このとき$\int_{R}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx$$=$$\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx+\int_{-\infty}^{0}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx$

95

$=$

$\Leftarrow x=\frac{1}{y}$$\int_{\infty}^{0}\frac{1}{\sqrt[3]{\frac{1}{y^{2}}(\frac{a}{y^{2}}+\frac{b}{y}+c)^{2}}}(-\frac{dy}{y^{2}})+\int_{0}^{-\infty}\frac{1}{\sqrt[3]{\frac{1}{y^{2}}(\frac{a}{y^{2}}+\frac{b}{y}+c)^{2}}}(-\frac{dy}{y^{2}}I$$=$$\int_{0}^{\infty}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy+\int_{-\infty}^{0}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy$$=$$\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy$

$=$

$=$

$\frac{\sqrt[3]{2}B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{c^{2}(4ac-b^{2})}}\Leftarrow$(15) $(carrow b;barrow c;aarrow d)$(16)$\frac{\sqrt[3]{2}B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{-D}}$となる。
次にその実根を$\alpha$とすると$ax^{3}+bx^{2}+cx+d=(x-\alpha)(ax^{2}+kx+l)$ ;$a\alpha^{3}+b\alpha^{2}+c\alpha+d=0$と表せる。 これを展開すると$b=k-a\alpha$ ,$c=l-k\alpha$ ,(17)$d=-l\alpha$を得る。 このとき$\int_{R}\frac{1}{\sqrt[3]{(x-\alpha)^{2}(ax^{2}+kx+l)^{2}}}dx$$=$$\int_{R}\frac{1}{\sqrt[3]{y^{2}\{a(y+\alpha)^{2}+k(y+\alpha)+l\}^{2}}}dy\Leftarrow x=y+\alpha$$=$$\int_{R}\frac{1}{\sqrt[3]{y^{2}\{ay^{2}+(2a\alpha+k)y+(a\alpha^{2}+k\alpha+l)\}^{2}}}dy$$=$$\frac{\sqrt[3]{2}B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{(a\alpha^{2}+k\alpha+l)^{2}\{4a(a\alpha^{2}+k\alpha+l)-(2a\alpha+k)^{2}\}}}\Leftarrow$$=$$\frac{\sqrt[3]{2}B(\frac{1}{2l},\frac{1}{26})}{\sqrt[6]{(a\alpha^{2}+k\alpha+)(4al-k^{2})}}$

(16)
(18)

となる。

Key Lemma (17) より等式$(a\alpha^{2}+k\alpha+l)^{2}(4al-k^{2})$が成り立つ。

$=$$27a^{2}d^{2}+4ac^{3}-18abcd-b^{2}c^{2}+4b^{3}d$

$=$

$-D$

(19)

96

証明は straightforward だが、 長いので省略する。
以上のことから (18) とMain Lemma より(20)$\int_{R}\frac{1}{\sqrt[3]{(x-\alpha)^{2}(ax^{2}+kx+l)^{2}}}dx=\frac{\sqrt[3]{2}B(\frac{1}{2},\frac{1}{6})}{\sqrt[6]{-D}}$を得る。 これで (11) の証明が完成した。

(B) の証明: $D>0$ の場合に (12) を証明する。

まず $\alpha>0$ として積分$\int_{R}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx$を考える。 このとき$\int_{R}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx$$=$$\int_{-\infty}^{0}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx+\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx$$=$$\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(x+\alpha)^{2}}}dx+\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx$となる。 ここで右辺の最初の積分に対して変数変換 $xarrow-x$ を行った。

各項を計算しよう。
$\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(x+\alpha)^{2}}}dx$$=$$\frac{1}{\sqrt[3]{\alpha}}\int_{0}^{\infty}\frac{1}{\sqrt[3]{t^{2}(t+1)^{2}}}dt\Leftarrow x=\alpha t$$=$$\alpha^{-1}3\int_{0}^{\infty}\frac{t^{-2}3}{(t+1)32}dt$$=$$\alpha^{-1}3B(\frac{1}{3}, \frac{1}{3})\Leftarrow(14)$で、 他方$\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx$$\Leftarrow x=$オ$=$$\alpha^{-1}3\int_{0}^{\infty}\frac{1}{\sqrt[3]{t^{2}(t-1)^{2}}}dt$$=$$\alpha^{-1}3\{\int_{0}^{1}\frac{1}{\sqrt[3]{t^{2}(t-1)^{2}}}dt+\int_{1}^{\infty}\frac{1}{\sqrt[3]{t^{2}(t-1)^{2}}}dt\}$$=$$\alpha^{-1}3\{\int_{0}^{1}\frac{1}{\sqrt[3]{t^{2}(1-t)^{2}}}dt+\int_{1}^{\infty}\frac{1}{\sqrt[3]{t^{2}(t-1)^{2}}}dt\}$$=$$2 \alpha^{-\perp}3\int_{0}^{1}\frac{1}{\sqrt[3]{t^{2}(1-t)^{2}}}dt$$\alpha$

(21)

97

$=$

$2 \alpha^{-\frac{1}{3}}\int_{0}^{1}t^{-\frac{2}{3}}(1-t)^{-\frac{2}{3}}dt$$=$$2 \alpha^{-\frac{1}{3}}B(\frac{1}{3}, \frac{1}{3})$となる。 ここで$\int_{1}^{\infty}\frac{1}{\sqrt[3]{t^{2}(t-1)^{2}}}dt$

$=$

$=$

$\int_{1}^{0}\frac{}{\sqrt[3]{\frac{1}{s^{2}}\frac{1(1-s)^{2}}{s^{2}}}}(-\frac{ds}{s^{2}})\Leftarrow t=\frac{1}{s}$$\int_{0}^{1}\frac{1}{\sqrt[3]{s^{2}(1-s)^{2}}}ds=\int_{0}^{1}\frac{1}{\sqrt[3]{t^{2}(1-t)^{2}}}dt$を使った。 (21) より$\int_{R}\frac{1}{\sqrt[3]{x^{2}(x-\alpha)^{2}}}dx=3\alpha^{-\frac{1}{3}}B(\begin{array}{l}11\overline{3}’\overline{3}\end{array})=\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[3]{\alpha}}$

(22)

を得る。
次に 3 次の多項式 $ax^{3}+bx^{2}+cx+d$ で、 $a=0$ の場合を考える。 このとき $D=$$b^{2}(c^{2}-4bd)>0$であるから$\int_{R}\frac{1}{\sqrt[3]{(bx^{2}+cx+d)^{2}}}dx$$=$$\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{\{b(x+\frac{c}{2b})^{2}-\frac{c^{2}-4bd}{4b}\}^{2}}}dx$$=$$\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{(bx^{2}-\frac{c^{2}-4bd}{4b})^{2}}}dx$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{(x^{2}-\alpha^{2})^{2}}}dx\Leftarrow\alpha^{2}=\frac{c^{2}-4bd}{4b^{2}}>0$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{(x-\alpha)^{2}(x+\alpha)^{2}}}dx$$=$$\frac{1}{\sqrt[3]{b^{2}}}\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{y^{2}(y-2\alpha)^{2}}}dy\Leftarrow y=x+\alpha$$=$$\frac{1}{\sqrt[3]{b^{2}}}\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[3]{2\alpha}}\Leftarrow(22)$$=$$\frac{3B(1\tilde{3}\frac{1}{32})}{\sqrt[3]{2\alpha b}}$$3B( \frac{1}{3}, \frac{1}{3})$$=$$\overline{\sqrt[6]{4\alpha^{2}b^{4}}}$$=$$\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[6]{b^{2}(c^{2}-4bd)}}\Leftarrow\alpha^{2}=\frac{c^{2}-4bd}{4b^{2}}$$=$$\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[6]{D}}$

(23)

98

を得る。
条件 $D>0$ より $ax^{3}+bx^{2}+cx+d=0$ は一つの実根をもつ。最初にそれを $x=0$ とする。 この場合 $d=0$ で$ax^{3}+bx^{2}+cx+d=x(ax^{2}+bx+c)$となる。 このとき $D=c^{2}(b^{2}-4ac)$ で$\int_{R}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx$$=$$\int_{0}^{\infty}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx+\int_{-\infty}^{0}\frac{1}{\sqrt[3]{x^{2}(ax^{2}+bx+c)^{2}}}dx$$=$$\Leftarrow x=\frac{1}{y}$$\int_{\infty}^{0}\frac{1}{\sqrt[3]{\frac{1}{y^{2}}(\frac{a}{y^{2}}+\frac{b}{y}+c)^{2}}}(-\frac{dy}{y^{2}})+\int_{0}^{-\infty}\frac{1}{\sqrt[3]{\frac{1}{y^{2}}(\frac{a}{y^{2}}+\frac{b}{y}+c)^{2}}}(-\frac{dy}{y^{2}})$$=$$\int_{0}^{\infty}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy+\int_{-\infty}^{0}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy$$=$$\int_{-\infty}^{\infty}\frac{1}{\sqrt[3]{(cy^{2}+by+a)^{2}}}dy$$=$$\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[6]{c^{2}(b^{2}-4ac)}}\Leftarrow$$=$(23) $(carrow b;barrow c;aarrow d)$(24)$\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[6]{D}}$となる。
次にその実根を$\alpha$とすると$ax^{3}+bx^{2}+cx+d=(x-\alpha)(ax^{2}+kx+l)$ ;$a\alpha^{3}+b\alpha^{2}+c\alpha+d=0$と表せ、 これを展開して$b=k-a\alpha$ ,$c=l-k\alpha$ ,$d=-l\alpha$を得る。 このとき$\int_{R}\frac{1}{\sqrt[3]{(x-\alpha)^{2}(ax^{2}+kx+l)^{2}}}dx$$=$$\int_{R}\frac{1}{\sqrt[3]{y^{2}\{a(y+\alpha)^{2}+k(y+\alpha)+l\}^{2}}}dy\Leftarrow x=y+\alpha$$=$$\int_{R}\frac{1}{\sqrt[3]{y^{2}\{ay^{2}+(2a\alpha+k)y+(a\alpha^{2}+k\alpha+l)\}^{2}}}dy$$=$$\frac{3B(\frac{1}{3},\frac{1}{3)})}{\sqrt[6]{(a\alpha^{2}+k\alpha+l)^{2}\{(2a\alpha+k^{2}-4a(a\alpha^{2}+k\alpha+l)\}}}\Leftarrow$

(24)

99

$=$

$\frac{3B(\frac{1}{3},\frac{1}{)3})}{\sqrt[6]{(a\alpha^{2}+k\alpha+l^{2}(k^{2}-4al)}}$$=$$\frac{3B(\frac{1}{3},\frac{1}{3})}{\sqrt[6]{D}}\Leftarrow$

(19)

(25)

となる。 これで (12) の証明が完成した。

(C) の証明: (13) の証明をする。

少し準備をする。 ガンマ関数$\Gamma(x)=\int_{0}^{\infty}e^{-t}t^{x-1}dt$$(x>0)$(26)$(x, y>0)$(27)$(0<x<1)$(28)に対して公式$B(x, y)$$=$$\Gamma(x)\Gamma(1.-x)$$=$$\Gamma(\frac{x}{2})\Gamma(\frac{x+1}{2})$$=$$\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}$$\frac{\pi}{\sin(\pi x)}$$\frac{\sqrt{}\overline{\pi}}{2^{x-1}}\Gamma(x)=2^{1-x}\Gamma(\frac{1}{2})\Gamma(x)$(29)

が成り立つ。

最後の式はルジャンドルの関係式 と呼ばれている。 この式で特に $x=2/3$ を代入す
ると$\Gamma(\frac{1}{3})\Gamma(\frac{5}{6})=\sqrt[3]{2}\Gamma(\frac{1}{2})\Gamma(\frac{2}{3})$である。両辺に $\Gamma(1/6)$ を掛け、計算を続けると$\Gamma(\frac{1}{3})\Gamma(\frac{5}{6})\Gamma(\frac{1}{6})=\sqrt[3]{2}r(\frac{1}{2})\Gamma(\frac{2}{3})\Gamma(\frac{1}{6})$ニ$\Gamma(\frac{1}{3})\frac{\pi}{\sin(\begin{array}{l}x6\end{array})}=\sqrt[3]{2}\Gamma(\frac{1}{2})\Gamma(\frac{1}{6})\Gamma(\frac{2}{3})$ニ$2 \pi\Gamma(\frac{1}{3})=\sqrt[3]{2}\Gamma(\frac{1}{2})\Gamma(\frac{1}{6})\Gamma(\frac{2}{3})$$\Leftrightarrow$$2\pi_{\Gamma()}\Gamma=\lrcorner^{1}\mathbb{A}_{2}23\sqrt[3]{2}^{\Gamma\Gamma}m^{11}$$\Leftarrow(28)$$\Gamma(\begin{array}{l}23\end{array})$$\Leftrightarrow$$2 \pi\frac{\Gamma(\begin{array}{l}1q\end{array})\Gamma(\begin{array}{l}1q\end{array})}{\Gamma(\begin{array}{l}23\end{array})\Gamma(\begin{array}{l}23\end{array})\Gamma(\begin{array}{l}13\end{array})}=\sqrt[3]{2}B(\frac{1}{2}, \frac{1}{6})\Leftarrow(27)$$\Leftrightarrow$$2 \pi^{\frac{\Gamma(\begin{array}{l}1q\end{array})}{\Gamma(\begin{array}{l}23\end{array})x3}}=\sqrt[3]{2}B(\frac{1}{2}, \frac{1}{6})\Leftarrow(28)$$\Leftrightarrow$$\sqrt{3}B(\frac{1}{3}, \frac{1}{3})=\sqrt[3]{2}B(\begin{array}{ll}\underline{1} \underline{1}2’ 6 \end{array}) \Leftarrow(27)$100を得る。 これで (13) の証明が完成した。

謝辞 計算をチェックしてもらった舟橋久仁夫氏と鈴木達夫氏に感謝します。

参考文献[1] A. Morozov and Sh. Shakirov: Introduction to Integral Discriminants,
arXiv: 0903.2595 [math-ph].
[2] K. Fujii: Beyond Gaussian: A Comment, arXiv: 0905.1363 [math-ph].

[3 藤井一幸 :3 次方程式の一般論, レクチャーノート (横浜市大) .
[4 沼澤秀亮 :4 次方程式の一般論, 横浜市立大学卒業論文 (横浜市大) .
[5 佐竹一郎 : 線型代数学, 裳華房.

[6] E. アルティン: ガンマ関数入門, 上野建爾 [訳解説], 日本評論社.