Twitterで見かけた積分を一般化してみた
問題の積分です
あらすじ
今週月曜日、Twitterでとても興味深い積分を見かけました。
yujitomo.hatenablog.com
意外にも、これが適当な置換積分や留数定理では全然計算できない!
埒が明かないので、一旦休息をとろうと生協へと立ち寄ることにしました。
生協の書籍コーナーをぶらぶらしていると、興味深いタイトルの本が。
- 作者: 黒川信重
- 出版社/メーカー: 岩波書店
- 発売日: 2013/11/08
- メディア: 単行本
- この商品を含むブログ (3件) を見る
気になって読んでみると、なんと三角関数が(よくわからないが)ゼータ関数に通じているではないか!
そして、この本によって偶然にもある等式に出会う。Raabe's formulaである。
ガンマ関数の対数を積分するとは、(個人的には)なんと刺激の強いことか!
すぐに研究室に駆け込んで調べてみることに。
するとRaabe's formulaの節にBinet's second expression for logΓと呼ばれる積分が!
Gamma function - Wikipedia, the free encyclopedia
これは考えるべき問題の積分とそっくりではないか!
この公式を使って計算はできないだろうかと考えました。
そして、ついに計算できました。しかもある程度の一般化もできたのです。
この問題への面白い別解と一般化を思いついたので近いうちにブログで公開します。 RT @yjjtw: 積分計算の問題 - yujitomoのブログ https://t.co/yMV6j5SFy5
— した狼@3/4(金)夜Skype人狼参加 (@shitaro2016) 2016年2月29日
なんかpdf見れないという意見があった気がするからスクショをペタペタ貼っておいた
本記事では、問題の積分を一般化した式と導出過程を簡単に紹介することにします。
出発はlog Gamma
まず、次の積分を導出します。
なぜ、この積分を考えるかというと、問題の積分の形がBinet's second expression for logΓ(以下、Binetの公式と呼ぶことにする)
に類似しているからです。
本節の積分は、このBinetの公式に対して
と置くことで得られます。
特に、
と置けば
が、
と置けば
が、
と置けば
が、それぞれ得られます。
arctan(x^3)はarctan(x)の和
arctan(x^3)は次のように表せます。
この式が本記事のキモであり、個人的に興味深いと思った、自信のある等式です。
この等式を思いついた理由は以下のとおりです:
問題の積分をBinetの公式に帰着したい。そのためにはarctan(x^3)をarctan(x)になおす必要がある。だからarctan(x^3)を(無理矢理にでも)arctan(x)で表そうとしました。
さて、問題の積分を部分積分すると
という微分が現れました。そう、arctan(x^3)を微分すると1/(x^2+1)で部分分数分解できるのです!
つまり、arctan(x^3)の微分はarctan(x)の和の微分で表せることにここで気づいたのです。実際、
なので本節の等式は成り立ちます。
より具体的に
式
のlogΓΓを具体的に書き換えたくありません?実はnを偶奇に分ければ、複素数を使わずに表現することができます。
ΓΓから、直ちに式
を使えば計算できそうだなあと予想できます。
n=2m+1の場合(m=0,1,2,...)
ΓΓは、式
より
となります。
このうち、ΓΓの部分は式
から、
が得られます。
一方、総乗は2m個の積を(x-a)(x+a)の形になるように並べ替えることで
と計算できます。
したがって、求めるΓΓは
となります。
n=2mの場合(m=1,2,...)
ΓΓは、式
より
となります。
このうち、ΓΓの部分は式
から、
が得られます。
一方、総乗は2m個の積を(x-a)(x+a)の形になるように並べ替えることで
と計算できます。
したがって、求めるΓΓは
となります。
複素数を使わない表示
というワケで、偶奇で場合分けをすれば、積分
は次のように表すことができます。
n=2m+1の場合(m=0,1,2,...)
n=2mの場合(m=1,2,...)
特にn=1,2の場合
n=1の場合の結果は、確かにyujitomoさんの答えと一致します。
また、n=1,2の場合における積分は、例えば公式集、Table of Definite and Infinite Integralsのp.94に答えだけ載っています。*1
だから本記事で示した、一般化された例の積分は多分合ってるんじゃないかなあと。
謝辞
とても奥の深い積分を我々に提供してくださったyujitomoさんに感謝します。
特に、ガンマ関数がこれほどまでに表情豊かであることを知れてとても興奮しています。
本当にありがとうございます。
それから、検算に付き合ってくれた、私と同じ研究室に所属するTypolice氏に感謝します。
*1:内部生向け。この本はチーズケーキにあります。Table of definite and infinite integrals | 東京工業大学附属図書館 蔵書検索 (OPAC)。 というか、かなりの数の公式集をあたったが、n=2まで載っている本はこの本以外に見当たらなかった。 ちなみに、この本はタイプライターで書いたようなレイアウトなので、とても見づらい。