この問いはどうやら、基本的な数学関数のみでガンマ関数を実装する問題に落ち着くみたい。実装する前に、ガンマ関数についてあまり知識がなかったのでいろいろと調べてみた。日本語Wikipediaにはスターリングの近似なる項目があり、ここにある式は比較的簡単ながらそこそこの近似値が得られるようになっています。近似値の精度は(方法にもよるとはいえ)Excelのgammalnに劣り、納得できる水準ではありません。approximationという近似手法を紹介してくれた。若干他の手法に比べると劣る部分があるらしい。ついて調べてみたら、ご丁寧なことにPythonコードが載っていて、おかげさまで、すっかり実装する気がなくなってしまった。xhttp://en.wikipedia.org/wiki/Lanczos_approximation基本的に、よくわからないが、配列pの値を使ってうまーく近似しているんだと。その配列pはどうやって得られたのだろうか・・・。あの配列pの値を求めるのが、本筋なのではなかろうかと思ったが、今日はここまでとした。approximationっていう近似手法もあるらしい。紹介したすごい参考になったページにもコードが載ってた。Wikipediaに載ってたアルゴリズムの精度確認でも、ちゃんとした結果になるか、気になったので確認だけしてみた。(6+0j)あれれー?なんかおかしいぞー。この結果は、見なかったことに...はい、このエントリーのこと忘れました。いうことで、明日元気だったら、3問目に挑戦します。