■求積の多様性を考える(その24)

 ケルビン卿の銘言に「数学者とは

  ∫(-∞,∞)exp(-x^2)dx=√π

を1+1=2のように自明だと思っている人である」というのがある.

 ガウス積分をn次元に拡張し,

  I=∫(-∞,∞)exp(-x1^2+x2^2+・・・+xn^2)dx1dx2・・・dxn

を考えると∫(-∞,∞)exp(-x^2)dx=√πのn重積分より,直ちに

  I=π^(n/2)

を得ることができます.

 n次元ガウス積分を別の方法,直交座標でなく,極座標で求めてみましょう.

===================================

 球に相当するn次元の図形を超球と呼びます.n次元単位超球{x1^2+x2^2+・・・+xn^2≦1}の体積をVnとすると,V1=2(直径),V2=π(面積),V3=4π/3(体積)はご存知でしょう.また,単位超球の表面積Sn-1はnVn,半径rのn次元球の体積はvnr^n,表面積はnVnr^n-1となります.

 ガウス積分の被積分関数を原点を中心とする半径rの球面上で積分し,次にr=0からr=∞まで積分すると,半径rの球面上で被積分関数は一定値exp(-r2)をとり,表面積はnVnr^n-1ですから,

  I=∫(0,∞)exp(-r^2)nVnr^n-1dr

=nVn∫(0,∞)r^(n-1)exp(-r^2)dr

 z=r^2と変数変換するとdz=2rdrより

  I=nVn/2∫(0,∞)z^(n/2-1)exp(-z)dz

=Vnn/2Γ(n/2) n/2Γ(n/2)=Γ(n/2+1)

=VnΓ(n/2+1)

 したがって,

  Vn=π^(n/2)/Γ(n/2+1)

を得ることができます.この結果は,形式的にVn=π^(n/2)/(n/2)!と書くことができます. Γ(m+1)=m!

 これより,半径rのn次元超球の超体積は

  Vnr^n=(πr^2)^(n/2)/Γ(n/2+1)となります.

===================================

[まとめ]ここでは球殻分割法を用いましたが,関数によっては円筒分割法が有効のことがあります.

===================================