Processing math: 100%

【γ13】ガンマ関数の放物線経路によるブルゲの積分表示・ハンケル路(ガンマ関数の基礎13)

今日のテーマ

「ガンマ関数の基礎」シリーズ第13回です。前回は:

【γ12】ガンマ関数の逆数・見た目だけは簡単な積分表示

Bourguetの積分表示

Γ(z)=2azeasinπz0eax2(1+x2)z12cos[2ax+(2z1)arctanx]dx

Whittaker-Watsonのガンマ関数の例題に載っていた表示です。何に使えるのかは分かりませんが、導出すること自体は面白そうなので取り組んでみました。

古いですが有名な書物で、どんどん改訂版が出ています。前半は解析学一般、後半は特殊関数という内容で、網羅的に勉強できます。演習問題に解答がないのが昔ながらのものって感じ。2022/11/6現在、最新版は5th Editionで私も所有していますが、廉価な3rdとかでも十分かと。


A Course of Modern Analysis: fifth Edition


A Course of Modern Analysis: Third Edition

しかし自力ではなかなか導き出すことができませんでした。そこでこの公式が最初に出てきた原論文に当たりました。それは「Acta Mathematica」という雑誌の第1巻(1882年)。Bourguet(ブルゲ)という人が書いた"Sur quelques intégrales définies"と題する論文です。

フランス語?なので全然読めませんが、数式を追うだけなら何とかなります。積分路をどうするかが分かれば、あとは加法定理だとか和積の公式だとか、三角関数の計算練習みたいなものです。

ハンケルの積分表示

スタートとなるのは過去記事でも紹介した「ハンケルの積分表示」です。

Γ(z)=12isinπzCet(t)z1dt(zZ)

経路 C

経路 C を簡単に説明します。t=reiθ (0<θ<2π) とすると、まず r=,θ=+0 から原点へやってきて、原点を1周して最後は r=,θ=2π0 となるような曲線です。これを満たしていれば形が変化してもよく、例えば(1)を導出する際には下のような経路を使いました。

このあたりの話については下の記事を参照:

【γ11】ガンマ関数の積分表示導出①(ハンケルとか)

この事実を使ってブルゲの表示を導きましょう。

放物線の経路

積分路の設定

(1)の積分においてt=2acosθ+isinθ1+cosθ(π<θ<π)とします。t ではなく t で置いていることに注意してください。偏角 θ を変化させると下図のような放物線を描きます(t の複素平面)。

t 平面で描かれる放物線

x=2acosθ1+cosθ , y=2asinθ1+cosθ とおくとcosθ=x2a1x2a,sinθ=y2a1x2a2乗して辺々足すと1であることから4a24ax=y2となり、原点を焦点とする放物線であることが分かる。

値の計算

微小要素はdt=2asinθi(1+cosθ)(1+cosθ)2dθよってゴリゴリ計算していくとCet(t)z1dt=ππ(2acosθ+isinθ1+cosθ)z1e2acosθ+isinθ1+cosθ2asinθi(1+cosθ)(1+cosθ)2dθ=2zazππ(cosθ+isinθ)z1(1+cosθ)z+1e2acosθ+isinθ1+cosθ[sinθi(1+cosθ)]dθ=2zazππ[cos(z1)θ+isin(z1)θ][sinθi(1+cosθ)](1+cosθ)z+1e2acosθ+isinθ1+cosθdθ=2zazππsinzθ+sin(z1)θi(coszθ+cos(z1)θ)(1+cosθ)z+1e2acosθ+isinθ1+cosθdθ=2zazππ2sin2z12θcosθ22icos2z12θcosθ22z+1cos2z+2θ2e2acosθ+isinθ1+cosθdθ=azππsin2z12θicos2z12θcos2z+1θ2e2acosθ+isinθ1+cosθdθ=azππsin2z12θicos2z12θcos2z+1θ2e2a(111+cosθ+isinθ1+cosθ)dθ=azeaππsin2z12θicos2z12θcos2z+1θ2eatan2θ2ei2atanθ2dθ=azeaππ(sin2z12θicos2z12θ)(cos2atanθ2+isin2atanθ2)cos2z+1θ2eatan2θ2dθ=azeaππsin(2z12θ+2atanθ2)icos(2z12θ+2atanθ2)cos2z+1θ2eatan2θ2dθ分子第1項は奇関数、分子第2項は偶関数なので=2iazeaπ0cos(2z12θ+2atanθ2)cos2z+1θ2eatan2θ2dθx=tanθ2 とおくと dθ=2cos2θ2dx ですから=4iazea0cos[(2z1)arctanx+2ax](1+x2)z12eax2dxこれとΓ(z)=12isinπzCet(t)z1dtを合わせることで結論を得ます。

Γ(z)=2azeasinπz0eax2(1+x2)z12cos[2ax+(2z1)arctanx]dx

冒頭で紹介したブルゲの論文では a=1/2 としたときの同じ表示を導出していますので、ここで示した公式はその一般化といえます。三角関数の計算がごちゃごちゃしましたが、たどり着くと達成感があるものです。

0ex2cos(2x4arctanx)(1+x2)2dx=2π3e

(2)で a=1 , z=3/2 とすると導出できます。Wolframで計算すると両辺ともに約0.4347で一致しています。

ここ数回でガンマ関数や対数ガンマ関数の積分表示を見てきました。次回からはディガンマ関数の積分表示をガンガン導出していきましょう!

次回:

応援のおねがい

Please support me!

まめしば
まめしば

記事を気に入って下さった方、「応援してあげてもいいよ」という方がいらっしゃったら15円から可能なので支援していただければ幸いです。情報発信を継続していくため、サーバー維持費などに充てさせていただきます。

ご支援いただいた方は、こちらで確認できます。

Amazonギフトの場合、
Amazonギフト券- Eメールタイプ – Amazonベーシック
より、金額は空白欄に適当に(15円から)書きこんで下さい。受取人は「mamekebiamazonあっとgmail.com」です(あっとは@に置き換えてください)。贈り主は「匿名」等でOKです。全額がクリエイターに届きます。

OFUSEは登録不要で、100円から寄付できます。金額の90%がクリエイターに届きます。

codocは登録不要で、100円から寄付できます。金額の85%がクリエイターに届きます。

Remaining:
codocの場合

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA