python - ガンマ(x+1/2)/ガンマ(x)の計算

Question

かなり大きな x に対して Gamma(x+1/2)/Gamma(x) を計算する必要があります。http://docs.scipy.org/doc/scipy/reference/generated/scipy.special.gamma.htmlを使用するだけでは、分母と分子が巨大であるため失敗します。ただし、分数自体は sqrt(x) 前後です。この分数を正確に計算するにはどうすればよいですか?

Answer 1

あなたが使用することができますbeta(a, b) = gamma(a) * gamma(b) / gamma(a+b)

from scipy.special import gamma, beta
x = 10
print gamma(x+0.5)/gamma(x)
print gamma(0.5)/beta(x, 0.5)