불완전 감마 함수란?
불완전 감마 함수는 적분을 무한대까지 끝까지 수행하는 대신 유한한 지점에서 멈춘다는 점에서 일반적인 (완전) 감마 함수를 일반화한 것입니다. 하부 불완전 감마 함수 \(\gamma(a,x)\)는 \(t^{a-1} e^{-t}\)를 0부터 x까지 적분하고, 상부 불완전 감마 함수 \(\Gamma(a,x)\)는 x부터 무한대까지 적분합니다. 두 함수 모두 형상 모수 a와 인수 x에 의존하며, 단위가 없는 순수 실수 값을 갖습니다. 통계학(카이제곱 분포와 감마 분포의 누적분포함수 CDF), 물리학, 신뢰성 공학, 대기행렬 이론 등에서 끊임없이 등장합니다.
계산기 사용법
형상 모수 a(반드시 양수, \(a > 0\))와 인수 x(반드시 0 이상, \(x \ge 0\))를 입력하세요. 계산기는 \(\gamma(a,x)\), \(\Gamma(a,x)\), 그리고 완전 감마 함수 \(\Gamma(a)\)를 함께 반환하므로 \(\gamma(a,x) + \Gamma(a,x) = \Gamma(a)\)라는 항등식을 직접 확인할 수 있습니다. \(x = 0\)일 때 하부 함수는 0이고 상부 함수는 \(\Gamma(a)\)와 같으며, x가 커질수록 하부 함수는 \(\Gamma(a)\)에 가까워지고 상부 함수는 0에 가까워집니다.
공식과 알고리즘
정의가 되는 적분은 $$\gamma\!\left(a,\, x\right) = \int_{0}^{x} t^{\,a-1}\, e^{-t}\, dt$$ 이고, $$\Gamma\!\left(a,\, x\right) = \int_{x}^{\infty} t^{\,a-1}\, e^{-t}\, dt = \Gamma\!\left(a\right) - \gamma\!\left(a,\, x\right)$$ 으로 동일한 피적분함수를 적분한 것입니다. 이를 수치적으로 안정하게 계산하기 위해 이 도구는 정규화된 형태인 \(P(a,x) = \gamma(a,x)/\Gamma(a)\)와 \(Q(a,x) = \Gamma(a,x)/\Gamma(a)\)를 사용합니다. \(x < a+1\)일 때는 빠르게 수렴하는 멱급수로 P를 구하고, 그렇지 않으면 렌츠(Lentz) 연분수로 Q를 구합니다. 완전 감마 \(\Gamma(a)\)는 \(\ln \Gamma(a)\)에 대한 랜초스(Lanczos) 근사로 얻습니다. 이는 고전적인 Numerical Recipes의 gammp/gammq 분기 방식으로, 배정밀도(double precision)에서 약 15자리 유효숫자까지 정확합니다.
계산 예시
\(a = 1\), \(x = 2\)를 예로 들어 보겠습니다. \(t^{a-1} = t^0 = 1\)이므로 하부 함수는 \(e^{-t}\)를 0부터 2까지 적분한 값, 즉 $$1 - e^{-2} = 1 - 0.13533528 = 0.86466472$$ 가 됩니다. 상부 함수는 \(e^{-2} = 0.13533528\)이고 \(\Gamma(1) = 1\)입니다. 항등식 확인: $$0.86466472 + 0.13533528 = 1.0$$ 으로 결과가 일치함을 알 수 있습니다.
자주 묻는 질문
왜 a는 양수여야 하나요? 수렴하는 정의와 랜초스 ln-gamma 계산에는 \(a > 0\)이 필요합니다. 0 이하의 정수에서는 \(\Gamma(a)\)에 극(pole)이 존재하기 때문입니다.
x가 0이면 어떻게 되나요? \(\gamma(a,0) = 0\)이고 \(\Gamma(a,0) = \Gamma(a)\)이므로, 상부 함수는 완전 감마 함수와 같아집니다.
결과는 얼마나 정확한가요? 배정밀도 연산과 급수/연분수 분기 방식으로 유효 정의역 전반에서 약 15자리 유효숫자의 정확도를 제공합니다.
