関数 y = exp(−x 2 )x 軸で囲まれた部分の面積 (= √ π ) がガウス積分を表す。 ガウス積分 （ガウスせきぶん、英 : Gaussian integral ）あるいはオイラー＝ポアソン積分 （オイラーポアソンせきぶん、英 : Euler–Poisson integral [ 1] ガウス関数 exp(−x 2 ) の実数全体での広義積分 ：
∫
−
∞
+
∞
e
−
x
2
d
x
=
π
{\displaystyle \int _{-\infty }^{+\infty }e^{-x^{2}}\,{\rm {d}}x={\sqrt {\pi }}}
のことである。名称は、数学・物理学者のカール・フリードリヒ・ガウス に由来する。
この積分の応用は広い。例えば、変数の微小変化に伴う正規分布 の正規化定数 の計算に用いられる。積分の上の限界を有限な値に替えることで、誤差関数 や正規分布の累積分布関数 とも深く関連する。
誤差関数を表す初等関数 はリッシュのアルゴリズム により存在しないことが証明できるが、ガウス積分の値は微分積分学 の道具立てを用いて解析的に求めることが可能である。つまり、初等関数としての不定積分
∫
e
−
x
2
d
x
{\displaystyle \textstyle \int e^{-x^{2}}\,{\rm {d}}x}
定積分
∫
−
∞
+
∞
e
−
x
2
d
x
{\displaystyle \textstyle \int _{-\infty }^{+\infty }e^{-x^{2}}\,{\rm {d}}x}
ガウス積分は物理学で非常に頻繁に現れ、またガウス積分の様々な一般化が場の量子論 に現れる。
ガウス積分を求める標準的な方法として、以下のアイデアはポアソン まで遡れる[ 2]
平面 R 2 exp{−(x 2 + y 2 )} = exp(−r 2 ) を考え、これを2通りの方法で計算する。
一つは直交座標系 に関する二重積分 として計算し、その値は求める値の平方になることを確かめる。
いま一つは極座標系 に関する二重積分（いわゆるバウムクーヘン積分 ）として計算し、その値が π 広義積分が現れることに注意して、これら2つの計算を比較して積分の値が求まる。即ち、面積要素 dA が xy -直交座標系では dA = dx dy , rθ -極座標系では dA = r dr dθ で与えられることに注意すれば、
∫
R
2
exp
(
−
(
x
2
+
y
2
)
)
d
A
=
∫
−
∞
∞
∫
−
∞
∞
e
−
(
x
2
+
y
2
)
d
x
d
y
=
(
∫
−
∞
∞
exp
(
−
t
2
)
d
t
)
2
,
{\displaystyle {\begin{aligned}\int _{\mathbb {R} ^{2}}\exp(-(x^{2}+y^{2}))\,dA=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }e^{-(x^{2}+y^{2})}\,dx\,dy=\left(\int _{-\infty }^{\infty }\exp(-t^{2})\,dt\right)^{2},\end{aligned}}}
と計算でき（フビニの定理を用いる。そのための正当性については下記を参照）、および
∫
R
2
exp
(
−
(
x
2
+
y
2
)
)
d
A
=
∫
0
2
π
∫
0
∞
exp
(
−
r
2
)
r
d
r
d
θ
=
2
π
∫
0
∞
r
exp
(
−
r
2
)
d
r
=
π
∫
−
∞
0
e
s
d
s
=
π
{\displaystyle {\begin{aligned}\int _{\mathbb {R} ^{2}}\exp(-(x^{2}+y^{2}))\,dA&=\int _{0}^{2\pi }\int _{0}^{\infty }\exp(-r^{2})\,r\,dr\,d\theta \\&=2\pi \int _{0}^{\infty }r\exp(-r^{2})\,dr=\pi \int _{-\infty }^{0}e^{s}\,ds=\pi \end{aligned}}}
と計算できる。後者では s = −r 2 ds = −2r dr となることを用いている。さてこれらの結果から
(
∫
−
∞
∞
exp
(
−
x
2
)
d
x
)
2
=
π
{\displaystyle \left(\int _{-\infty }^{\infty }\exp(-x^{2})\,dx\right)^{2}=\pi }
であり、符号を考慮して
∫
−
∞
∞
exp
(
−
x
2
)
d
x
=
π
{\displaystyle \int _{-\infty }^{\infty }\exp(-x^{2})\,dx={\sqrt {\pi }}}
を得る。
上記の考察において、広義二重積分や二つの式を等しいとおいたことに対する正当性を再考しておこう。まずは近似函数
I
(
a
)
=
∫
−
a
a
exp
(
−
x
2
)
d
x
{\displaystyle I(a)=\int _{-a}^{a}\exp(-x^{2})\,dx}
を考える。求めるガウス積分が絶対収斂 ならば、それはコーシー主値 、即ち
lim
a
→
∞
I
(
a
)
{\displaystyle \lim _{a\to \infty }I(a)}
なる極限によって求められることになる。これを見るには、
∫
−
∞
∞
|
exp
(
−
x
2
)
|
d
x
<
∫
−
∞
−
1
−
x
exp
(
−
x
2
)
d
x
+
∫
−
1
1
exp
(
−
x
2
)
d
x
+
∫
1
∞
x
exp
(
−
x
2
)
d
x
<
∞
{\displaystyle \int _{-\infty }^{\infty }|\exp(-x^{2})|\,dx<\int _{-\infty }^{-1}-x\exp(-x^{2})\,dx+\int _{-1}^{1}\exp(-x^{2})\,dx+\int _{1}^{\infty }x\exp(-x^{2})\,dx<\infty }
が成り立つという事実を確かめればよい。故に I (a )
I
(
a
)
2
=
∫
−
a
a
∫
−
a
a
exp
(
−
(
x
2
+
y
2
)
)
d
x
d
y
{\displaystyle I(a)^{2}=\int _{-a}^{a}\int _{-a}^{a}\exp(-(x^{2}+y^{2}))\,dx\,dy}
と書くことができて、フビニの定理 により、これは xy -座標平面 における面積分
∫
exp
(
−
(
x
2
+
y
2
)
)
d
A
{\displaystyle \int \exp(-(x^{2}+y^{2}))\,dA}
に等しいことが確かめられる。ただし、積分域は {(−a , a ), (a , a ), (a , −a ), (−a , −a )} を頂点集合とする正方形である。
指数函数は全実数に対して正の値を取るから、上記の積分域の内接円 上での積分は I (a )2 外接円 上での積分は I (a )2
x
=
r
cos
θ
,
y
=
r
sin
θ
,
d
x
d
y
=
r
d
r
d
θ
{\displaystyle x=r\cos \theta ,\quad y=r\sin \theta ,\quad dx\,dy=rdr\,d\theta }
なる標準的な変換でうつれば容易に計算できるから、積分を実行して
π
(
1
−
e
−
a
2
)
<
I
(
a
)
2
<
π
(
1
−
e
−
2
a
2
)
{\displaystyle \pi (1-e^{-a^{2}})<I(a)^{2}<\pi (1-e^{-2a^{2}})}
なる評価を得ることができる。a → ∞
(
∫
−
∞
∞
exp
(
−
x
2
)
d
x
)
2
=
π
{\displaystyle \left(\int _{-\infty }^{\infty }\exp(-x^{2})\,dx\right)^{2}=\pi }
が正当化できる。
ガウス積分を計算する別な方法として、以下はラプラス (1812) にまで遡れる[ 2]
y
=
x
s
,
d
y
=
x
d
s
{\displaystyle y=xs,\quad dy=x\,ds}
と置くと、y を ±∞ へ近づけるとき s の極限は x の符号で決まるから、exp(−x 2 ) が偶函数 ゆえに実数全体にわたる積分が正の実数全体にわたる積分の 2 倍となること、つまり
∫
−
∞
∞
e
−
x
2
d
x
=
2
∫
0
∞
e
−
x
2
d
x
{\displaystyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx=2\int _{0}^{\infty }e^{-x^{2}}\,dx}
であることを利用すれば計算が簡単になる。即ち、積分範囲を x ≥ 0y と s とは同じ極限を持ち、
I
2
=
4
∫
0
∞
∫
0
∞
e
−
(
x
2
+
y
2
)
d
y
d
x
{\displaystyle I^{2}=4\int _{0}^{\infty }\int _{0}^{\infty }e^{-(x^{2}+y^{2})}dy\,dx}
が成り立つ。故に
I
2
4
=
∫
0
∞
(
∫
0
∞
e
−
(
x
2
+
y
2
)
d
y
)
d
x
=
∫
0
∞
(
∫
0
∞
e
−
x
2
(
1
+
s
2
)
x
d
s
)
d
x
=
∫
0
∞
(
∫
0
∞
e
−
x
2
(
1
+
s
2
)
x
d
x
)
d
s
=
1
2
∫
0
∞
d
s
1
+
s
2
=
1
2
arctan
s
|
0
∞
=
π
4
{\displaystyle {\begin{aligned}{\frac {I^{2}}{4}}&=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-(x^{2}+y^{2})}\,dy\right)\,dx=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}(1+s^{2})}x\,ds\right)dx\\[5pt]&=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}(1+s^{2})}x\,dx\right)\,ds={\frac {1}{2}}\int _{0}^{\infty }{\frac {ds}{1+s^{2}}}\\[5pt]&=\left.{\frac {1}{2}}\arctan s\,\right|_{0}^{\infty }={\frac {\pi }{4}}\end{aligned}}}
となり、所期の I = √ π
ウォリス積分 における公式を用いて証明することができる。
x
≥
0
{\displaystyle x\geq 0}
1
−
x
2
≤
e
−
x
2
≤
1
1
+
x
2
{\displaystyle 1-x^{2}\leq e^{-x^{2}}\leq {\frac {1}{1+x^{2}}}}
微分法 により示す自然数 n に対して
∫
0
π
/
2
sin
2
n
+
1
θ
d
θ
<
∫
0
∞
e
−
n
t
2
d
t
<
∫
0
π
/
2
sin
2
n
−
2
θ
d
θ
{\displaystyle \int _{0}^{\pi /2}\sin ^{2n+1}\theta \,d\theta <\int _{0}^{\infty }e^{-nt^{2}}\,dt<\int _{0}^{\pi /2}\sin ^{2n-2}\theta \,d\theta }
が成り立つことを示す
∫
−
∞
∞
e
−
x
2
d
x
=
2
n
∫
0
∞
e
−
n
t
2
d
t
{\displaystyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx=2{\sqrt {n}}\int _{0}^{\infty }e^{-nt^{2}}\,dt}
n → ∞ のとき √ π 被積分関数が偶関数 ゆえ
∫
−
∞
∞
exp
(
−
x
2
)
d
x
=
2
∫
0
∞
exp
(
−
x
2
)
d
x
{\displaystyle \int _{-\infty }^{\infty }\exp(-x^{2})\,dx=2\int _{0}^{\infty }\exp(-x^{2})\,dx}
が成り立ち、これに変数変換 x = t 1/2 を行えばオイラー積分
2
∫
0
∞
exp
(
−
x
2
)
d
x
=
2
∫
0
∞
1
2
exp
(
−
t
)
t
−
1
/
2
d
t
=
Γ
(
1
2
)
=
π
{\displaystyle 2\int _{0}^{\infty }\exp(-x^{2})\,dx=2\int _{0}^{\infty }{\frac {1}{2}}\exp(-t)\,t^{-1/2}\,dt=\Gamma \left({\frac {1}{2}}\right)={\sqrt {\pi }}}
が得られる。ここで Γ はガンマ関数 。この式は、半整数値の階乗 が √ π
∫
0
∞
exp
(
−
a
x
b
)
d
x
=
1
b
a
−
1
/
b
Γ
(
1
b
)
{\displaystyle \int _{0}^{\infty }\exp(-ax^{b})\,dx={\frac {1}{b}}a^{-1/b}\,\Gamma \left({\frac {1}{b}}\right)}
が成り立つ。
勝手なガウス函数 の積分は
∫
−
∞
∞
exp
(
−
a
x
2
)
d
x
=
π
a
Re
{
a
}
≥
0
,
a
≠
0
{\displaystyle \int _{-\infty }^{\infty }\exp \left(-a{x^{2}}\right)dx={\sqrt {\frac {\pi }{a}}}\qquad \operatorname {Re} \{a\}\geq 0,{a}\neq 0}
(1.1 )
あるいは
∫
−
∞
∞
exp
(
−
a
(
x
−
b
)
2
)
d
x
=
π
a
Re
{
a
}
>
0
,
b
∈
C
{\displaystyle \int _{-\infty }^{\infty }\exp \left(-a{(x-b)^{2}}\right)dx={\sqrt {\frac {\pi }{a}}}\qquad \operatorname {Re} \{a\}>0,b\in \mathbb {C} }
(1.2 )
で与えられる
[ 注釈 1]
A = (αij )正定値 対称 （従って可逆な ）共変行列（二階共変テンソル ）ならば
∫
−
∞
∞
exp
(
−
1
2
x
⊺
A
x
)
d
n
x
=
∫
−
∞
∞
exp
(
−
1
2
∑
i
,
j
=
1
n
α
i
j
x
i
x
j
)
d
n
x
=
(
2
π
)
n
det
A
{\displaystyle \int _{-\infty }^{\infty }\exp \left(-{\frac {1}{2}}{\boldsymbol {x}}^{\intercal }{\boldsymbol {A}}{\boldsymbol {x}}\right)\,d^{n}{\boldsymbol {x}}=\int _{-\infty }^{\infty }\exp \left(-{\frac {1}{2}}\sum _{i,j=1}^{n}\alpha _{ij}x_{i}x_{j}\right)\,d^{n}{\boldsymbol {x}}={\sqrt {\frac {(2\pi )^{n}}{\det {\boldsymbol {A}}}}}}
(2.1 )
が成り立つ[ 注釈 2] R n
また、
∫
x
k
1
⋯
x
k
2
N
exp
(
−
1
2
x
⊺
A
x
)
d
n
x
=
(
2
π
)
n
det
A
1
2
N
N
!
∑
σ
∈
S
2
N
(
A
−
1
)
k
σ
(
1
)
k
σ
(
2
)
⋯
(
A
−
1
)
k
σ
(
2
N
−
1
)
k
σ
(
2
N
)
{\displaystyle \int x^{k_{1}}\dotsb x^{k_{2N}}\exp \left(-{\frac {1}{2}}{\boldsymbol {x}}^{\intercal }{\boldsymbol {A}}{\boldsymbol {x}}\right)d^{n}{\boldsymbol {x}}={\sqrt {\frac {(2\pi )^{n}}{\det {\boldsymbol {A}}}}}\;{\frac {1}{2^{N}N!}}\sum _{\sigma \in S_{2N}}({\boldsymbol {A}}^{-1})^{k_{\sigma (1)}k_{\sigma (2)}}\dotsb ({\boldsymbol {A}}^{-1})^{k_{\sigma (2N-1)}k_{\sigma (2N)}}}
(2.2 )
が成立する。ここで、σ は {1, ..., 2N } の置換 であり、右辺に現れる余分な因子は A −1 N 個のコピーを {1, ..., 2N } の組合せ対 (combinatorial pairing) の全体に亘って加えた和を意味する。
あるいはまた、A −1 = (βij )
∫
f
(
x
)
exp
(
−
1
2
x
⊺
A
x
)
d
n
x
=
(
2
π
)
n
det
A
exp
(
1
2
∇
⊺
A
−
1
∇
)
f
(
x
)
|
x
=
0
=
(
2
π
)
n
det
A
exp
(
1
2
∑
i
,
j
=
1
n
β
i
j
∂
2
∂
x
i
∂
x
j
)
f
(
x
)
|
x
=
0
{\displaystyle \int f({\boldsymbol {x}})\exp \left(-{\frac {1}{2}}{\boldsymbol {x}}^{\intercal }{\boldsymbol {A}}{\boldsymbol {x}}\right)d^{n}{\boldsymbol {x}}={\sqrt {\frac {(2\pi )^{n}}{\det {\boldsymbol {A}}}}}\;\left.\exp \left({\frac {1}{2}}\nabla ^{\intercal }{\boldsymbol {A}}^{-1}\nabla \right)f({\boldsymbol {x}})\right|_{{\boldsymbol {x}}={\boldsymbol {0}}}={\sqrt {\frac {(2\pi )^{n}}{\det {\boldsymbol {A}}}}}\;\left.\exp \left({\frac {1}{2}}\sum _{i,j=1}^{n}\beta _{ij}{\frac {\partial ^{2}}{\partial x^{i}\partial x^{j}}}\right)f({\boldsymbol {x}})\right|_{{\boldsymbol {x}}={\boldsymbol {0}}}}
(2.3 )
がいくつかの解析関数 f に対して成立する。f は増加具合が適当に制限されているとかあるいはほかの技術的な判定条件を満足する必要がある。これは特定の関数に対してはうまく行くがそうでないものもある。たとえば多項式ならば成立する。また微分作用素変数の指数関数 exp は冪級数 として理解され、あらたな微分作用素を定めるものである。
さらに無限次元への一般化としての汎関数積分 には厳密な定義は無く、多くの場合それは計算的でさえないが、ガウス汎関数積分を有限次元の場合の類似物として「定義」することができる。もちろん問題はあって、単純に有限次元の場合の式を無限次元の場合に適用しようとすれば (2π )∞ は無限大に発散してしまうし、汎函数行列式 (functional determinant) も一般には無限大となりうる。これらのことを考慮して比
∫
f
(
x
1
)
⋯
f
(
x
2
N
)
exp
(
−
1
2
∬
A
(
x
2
N
+
1
,
x
2
N
+
2
)
f
(
x
2
N
+
1
)
f
(
x
2
N
+
2
)
d
x
2
N
+
1
d
x
2
N
+
2
)
D
f
∫
exp
(
−
1
2
∬
A
(
x
2
N
+
1
,
x
2
N
+
2
)
f
(
x
2
N
+
1
)
f
(
x
2
N
+
2
)
d
x
2
N
+
1
d
x
2
N
+
2
)
D
f
=
1
2
N
N
!
∑
σ
∈
S
2
N
A
−
1
(
x
σ
(
1
)
,
x
σ
(
2
)
)
⋯
A
−
1
(
x
σ
(
2
N
−
1
)
,
x
σ
(
2
N
)
)
{\displaystyle {\begin{aligned}&{\frac {\displaystyle \int f(x_{1})\dotsb f(x_{2N})\exp \left(-{\dfrac {1}{2}}\displaystyle \iint {\boldsymbol {A}}(x_{2N+1},x_{2N+2})f(x_{2N+1})f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2}\right){\mathcal {D}}f}{\displaystyle \int \exp \left(-{\dfrac {1}{2}}\displaystyle \iint {\boldsymbol {A}}(x_{2N+1},x_{2N+2})f(x_{2N+1})f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2}\right){\mathcal {D}}f}}\\[5pt]&={\frac {1}{2^{N}N!}}\sum _{\sigma \in S_{2N}}{\boldsymbol {A}}^{-1}(x_{\sigma (1)},x_{\sigma (2)})\dotsb {\boldsymbol {A}}^{-1}(x_{\sigma (2N-1)},x_{\sigma (2N)})\end{aligned}}}
(2.4 )
のみを考えることにするならばガウス汎関数積分を扱うことができるという意味である。ドヴィット記法 (deWitt notation) を使えば、この等式は有限次元の場合と同じ形に書くことができる。
A = (αij )
∫
exp
(
−
1
2
∑
i
,
j
=
1
n
α
i
j
x
i
x
j
+
∑
i
=
1
n
b
i
x
i
)
d
n
x
=
(
2
π
)
n
det
A
exp
(
1
2
t
b
A
−
1
b
)
{\displaystyle \int \exp \left(-{\frac {1}{2}}\sum _{i,j=1}^{n}\alpha _{ij}x_{i}x_{j}+\sum _{i=1}^{n}b_{i}x_{i}\right)d^{n}x={\sqrt {\frac {(2\pi )^{n}}{\det {A}}}}\exp \left({\frac {1}{2}}{}^{t}{\boldsymbol {b}}A^{-1}{\boldsymbol {b}}\right)}
(3.1 )
が成り立つ。ただし、b = (bi )t 行列の転置 とする。
同様の積分として、
∫
0
∞
x
2
n
exp
(
−
x
2
a
2
)
d
x
=
π
(
2
n
−
1
)
!
!
2
n
+
1
a
2
n
+
1
=
π
(
2
n
)
!
n
!
(
a
2
)
2
n
+
1
,
∫
0
∞
x
2
n
+
1
exp
(
−
x
2
a
2
)
d
x
=
n
!
2
a
2
n
+
2
{\displaystyle {\begin{aligned}\int _{0}^{\infty }x^{2n}\exp \left(-{\frac {x^{2}}{a^{2}}}\right)dx&={\sqrt {\pi }}{\frac {(2n-1)!!}{2^{n+1}}}a^{2n+1}={\sqrt {\pi }}{\frac {\left(2n\right)!}{n!}}\left({\frac {a}{2}}\right)^{2n+1},\\\int _{0}^{\infty }x^{2n+1}\exp \left(-{\frac {x^{2}}{a^{2}}}\right)dx&={\frac {n!}{2}}a^{2n+2}\end{aligned}}}
(4.1 )
が成立する。これらを導出するには積分記号下での微分法を用いるのが簡便である：
∫
−
∞
∞
x
2
n
exp
(
−
α
x
2
)
d
x
=
(
−
1
)
n
∫
−
∞
∞
∂
n
∂
α
n
e
−
α
x
2
d
x
=
(
−
1
)
n
∂
n
∂
α
n
∫
−
∞
∞
exp
(
−
α
x
2
)
d
x
=
π
(
−
1
)
n
∂
n
∂
α
n
α
−
1
/
2
=
π
α
(
2
n
−
1
)
!
!
(
2
α
)
n
.
{\displaystyle {\begin{aligned}&\int _{-\infty }^{\infty }x^{2n}\exp(-\alpha x^{2})\,dx=\left(-1\right)^{n}\int _{-\infty }^{\infty }{\frac {\partial ^{n}}{\partial \alpha ^{n}}}e^{-\alpha x^{2}}\,dx=\left(-1\right)^{n}{\frac {\partial ^{n}}{\partial \alpha ^{n}}}\int _{-\infty }^{\infty }\exp(-\alpha x^{2})\,dx\\[8pt]&={\sqrt {\pi }}(-1)^{n}{\frac {\partial ^{n}}{\partial \alpha ^{n}}}\alpha ^{-1/2}={\sqrt {\frac {\pi }{\alpha }}}{\frac {(2n-1)!!}{\left(2\alpha \right)^{n}}}.\end{aligned}}}
(4.2 )
被積分函数の冪指数がもっと別の偶数次多項式に変わった場合も、級数解は容易に計算することができる。例えば四次多項式を冪指数とする指数函数の積分は
∫
−
∞
∞
exp
(
α
x
4
+
β
x
3
+
γ
x
2
+
δ
x
+
ϵ
)
d
x
=
1
2
exp
(
ϵ
)
∑
n
,
m
,
p
=
0
n
+
p
=
0
mod
2
∞
β
n
n
!
γ
m
m
!
δ
p
p
!
Γ
(
3
n
+
2
m
+
p
+
1
4
)
exp
(
3
n
+
2
m
+
p
+
1
4
log
(
−
α
)
)
{\displaystyle {\begin{aligned}&\int _{-\infty }^{\infty }\exp(\alpha x^{4}+\beta x^{3}+\gamma x^{2}+\delta x+\epsilon )\,dx\\[5pt]&={\frac {1}{2}}\exp(\epsilon )\sum _{n,m,p=0 \atop n+p=0{\bmod {2}}}^{\infty }{\frac {\beta ^{n}}{n!}}\;{\frac {\gamma ^{m}}{m!}}\;{\frac {\delta ^{p}}{p!}}\;{\frac {\Gamma \left({\dfrac {3n+2m+p+1}{4}}\right)}{\exp \left({\dfrac {3n+2m+p+1}{4}}\log(-\alpha )\right)}}\end{aligned}}}
(5.1 )
で表される。n + p = 0 mod 2−∞ から 0 までの積分が各項に (−1)n + p なる因子として寄与し、0 から +∞ までの積分が各項に 1/2 の因子として寄与することによる。これらの積分は場の量子論 に属する話題である。
^ 式 1.1 ξ = (x + b )/c dx = c dξ 式 1.2 平方完成 すれば 式 1.1 ^ 式 2.1 A が対角行列 であれば一変数の場合と同様にして右辺を得られる（ただし積分が収束するために A の成分はすべて正、つまり A が正定値行列であることが要求される）。非対角項がゼロでない場合、A が実対称行列であるため、積分変数 x 直交行列 O を用いて変数変換することで行列 A を対角化 できる。対角化した後の計算は対角行列の場合と同様。
Weisstein, Eric W. "Gaussian Integral" . mathworld.wolfram.com (英語). David Griffiths. Introduction to Quantum Mechanics. 2nd Edition back cover.
Abramowitz, M. and Stegun, I. A. Handbook of Mathematical Functions , Dover Publications, Inc. New York 
![{\displaystyle {\begin{aligned}{\frac {I^{2}}{4}}&=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-(x^{2}+y^{2})}\,dy\right)\,dx=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}(1+s^{2})}x\,ds\right)dx\\[5pt]&=\int _{0}^{\infty }\left(\int _{0}^{\infty }e^{-x^{2}(1+s^{2})}x\,dx\right)\,ds={\frac {1}{2}}\int _{0}^{\infty }{\frac {ds}{1+s^{2}}}\\[5pt]&=\left.{\frac {1}{2}}\arctan s\,\right|_{0}^{\infty }={\frac {\pi }{4}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fab9f8f0b30ee0d11ef4da4ee61ab59052b335a7)
![{\displaystyle {\begin{aligned}&{\frac {\displaystyle \int f(x_{1})\dotsb f(x_{2N})\exp \left(-{\dfrac {1}{2}}\displaystyle \iint {\boldsymbol {A}}(x_{2N+1},x_{2N+2})f(x_{2N+1})f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2}\right){\mathcal {D}}f}{\displaystyle \int \exp \left(-{\dfrac {1}{2}}\displaystyle \iint {\boldsymbol {A}}(x_{2N+1},x_{2N+2})f(x_{2N+1})f(x_{2N+2})\,dx_{2N+1}\,dx_{2N+2}\right){\mathcal {D}}f}}\\[5pt]&={\frac {1}{2^{N}N!}}\sum _{\sigma \in S_{2N}}{\boldsymbol {A}}^{-1}(x_{\sigma (1)},x_{\sigma (2)})\dotsb {\boldsymbol {A}}^{-1}(x_{\sigma (2N-1)},x_{\sigma (2N)})\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d04824411f58a7260ac5e860c388807ca97ef302)
![{\displaystyle {\begin{aligned}&\int _{-\infty }^{\infty }x^{2n}\exp(-\alpha x^{2})\,dx=\left(-1\right)^{n}\int _{-\infty }^{\infty }{\frac {\partial ^{n}}{\partial \alpha ^{n}}}e^{-\alpha x^{2}}\,dx=\left(-1\right)^{n}{\frac {\partial ^{n}}{\partial \alpha ^{n}}}\int _{-\infty }^{\infty }\exp(-\alpha x^{2})\,dx\\[8pt]&={\sqrt {\pi }}(-1)^{n}{\frac {\partial ^{n}}{\partial \alpha ^{n}}}\alpha ^{-1/2}={\sqrt {\frac {\pi }{\alpha }}}{\frac {(2n-1)!!}{\left(2\alpha \right)^{n}}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8374681ffdd001a8d4f83fb562fac3ce351e064f)
![{\displaystyle {\begin{aligned}&\int _{-\infty }^{\infty }\exp(\alpha x^{4}+\beta x^{3}+\gamma x^{2}+\delta x+\epsilon )\,dx\\[5pt]&={\frac {1}{2}}\exp(\epsilon )\sum _{n,m,p=0 \atop n+p=0{\bmod {2}}}^{\infty }{\frac {\beta ^{n}}{n!}}\;{\frac {\gamma ^{m}}{m!}}\;{\frac {\delta ^{p}}{p!}}\;{\frac {\Gamma \left({\dfrac {3n+2m+p+1}{4}}\right)}{\exp \left({\dfrac {3n+2m+p+1}{4}}\log(-\alpha )\right)}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c29391800df98dc1d80d4f08fe89e5e93ab30446)
