Γ関数の数値計算
Fermi関数に出て来るΓ関数ですが、計算が大変なので少しお手伝いです。
変数が実数のもの(Γ(x))
dgamma.c
変数xを渡すとΓ関数の計算結果を返す関数
変数が虚数のもの(Γ(x+iy))
cdgamma.c
これは変数 x(実数部)と y(虚数部)を渡すと|Γ(x+iy)|の2乗を返すように オリジナルのプログラムを少しだけ修正しています。(構造体が分からなくても良いように)
Header file
gammaFunc.h
一応必要、かな。
オリジナルのプログラムは
大浦拓哉さんのページ
です。