๊ฐ๋ง ํจ์๋?
๊ฐ๋ง ํจ์๋ \(\Gamma(a)\)๋ก ํ๊ธฐํ๋ฉฐ, ํฉํ ๋ฆฌ์ผ์ ์ค์(๊ทธ๋ฆฌ๊ณ ๋ณต์์) ์์ญ๊น์ง ์ฐ์์ ์ผ๋ก ํ์ฅํ ํจ์์ ๋๋ค. ์์ ์ ์ n์ ๋ํด์๋ \(\Gamma(n) = (n-1)!\) ์ด ์ฑ๋ฆฝํ๋ฏ๋ก, ์๋ฅผ ๋ค์ด \(\Gamma(5) = 4! = 24\) ๊ฐ ๋ฉ๋๋ค. ์ค์๋ถ๊ฐ ์์์ธ ์ค์ a์ ๋ํด์๋ ๋ค์๊ณผ ๊ฐ์ ์ด์์ ๋ถ์ผ๋ก ์ ์๋ฉ๋๋ค: $$\Gamma\!\left(\text{a}\right) = \int_{0}^{\infty} t^{\,\text{a} - 1}\, e^{-t}\, dt$$ ์ด ๊ณ์ฐ๊ธฐ๋ ์ ๋ ฅํ ์์์ ์ค์ a์ ๋ํด \(\Gamma(a)\) ๊ฐ์ ๋๋ ค์ค๋๋ค.
๊ณ์ฐ๊ธฐ ์ฌ์ฉ ๋ฐฉ๋ฒ
"๋ณ์ a" ์นธ์ ์ค์ ์ธ์ a๋ฅผ ์ ๋ ฅํ๊ณ , ์์์ ์๋ ๋ช ์๋ฆฌ๊น์ง ํ์ํ ์ง ์ ํํ๋ฉด ๋ฉ๋๋ค. ํผ์ ๋ถ ํจ์ \(t^{a-1}e^{-t}\)์ ์ ๋ถ ๊ตฌ๊ฐ 0๋ถํฐ ๋ฌดํ๋๊น์ง๋ ์ ์์ ๋ฐ๋ผ ๊ณ ์ ๋์ด ์์ผ๋ฏ๋ก, ์ฌ์ฉ์๋ a ๊ฐ๋ง ๋ฃ์ผ๋ฉด ๋ฉ๋๋ค. ๊ณ์ฐ๊ธฐ๋ \(\Gamma(a)\)๋ฅผ ์ถ๋ ฅํฉ๋๋ค. ๋ง์ฝ a = 0 ๋๋ ์์ ์ ์๋ฅผ ์ ๋ ฅํ๋ฉด, ๊ทธ ์ง์ ์์ ๊ฐ๋ง ํจ์๊ฐ ๊ทน(pole)์ ๊ฐ์ง๋ฏ๋ก "์ ์๋์ง ์์"์ด๋ผ๊ณ ํ์๋ฉ๋๋ค.
๊ณต์ ํ์ด
์ด ๊ณ์ฐ๊ธฐ๋ ๋งค๋ฒ ์์น์ ์ผ๋ก ์ ๋ถํ๋ ๋์ Lanczos ๊ทผ์ฌ(g = 7, 9๊ฐ์ ๊ณ์)๋ฅผ ์ฌ์ฉํฉ๋๋ค. ์ด ๋ฐฉ๋ฒ์ ์ ๋ถ๊ฐ์ ์ฝ 15๊ฐ์ ์ ํจ ์ซ์๊น์ง ์ฌํํฉ๋๋ค. \(a \le 0.5\) ์ธ ๊ฒฝ์ฐ์๋ ๋จผ์ ๋ฐ์ฌ ๊ณต์ \(\Gamma(a)\cdot\Gamma(1-a) = \pi/\sin(\pi a)\)๋ฅผ ์ ์ฉํฉ๋๋ค. ์ด ๊ณต์์ ์ธ์๋ฅผ ์์น์ ์ผ๋ก ์์ ์ ์ธ ์์ญ์ผ๋ก ์ฎ๊ฒจ ์ฃผ๋ฉฐ, ์ ์๊ฐ ์๋ ์์ ์ธ์์์ ๋ํ๋๋ ์ ํํ(๋๋ก๋ ์์์ธ) ๊ฐ๊น์ง ์ ํํ ๊ณ์ฐํด ์ค๋๋ค.
ํ์ด ์์
a = 3.5์ธ ๊ฒฝ์ฐ๋ฅผ ์ดํด๋ด ์๋ค. ์ ํ์ \(\Gamma(a) = (a-1)\cdot\Gamma(a-1)\)์ ์ฌ์ฉํ๋ฉด: $$\Gamma(3.5) = 2.5 \cdot 1.5 \cdot 0.5 \cdot \Gamma(0.5) = 1.875 \cdot \sqrt{\pi} = 1.875 \cdot 1.7724538509 \approx 3.3233509704$$ ๊ณ์ฐ๊ธฐ๋ ๋์ผํ ๊ฐ์ ๋๋ ค์ค๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
\(\Gamma(0)\)์ ์ ์ ์๋์ง ์๋์? ์ ๋ถ์ด ๋ฐ์ฐํ๋ฉฐ, ์ด ํจ์๋ 0๊ณผ ๋ชจ๋ ์์ ์ ์์์ ๋จ์๊ทน(simple pole)์ ๊ฐ์ง๋ฏ๋ก ๊ฐ์ด ๋ฌดํ๋๊ฐ ๋ฉ๋๋ค.
\(\Gamma(0.5)\)๋ ์ผ๋ง์ธ๊ฐ์? ์ ํํ \(\sqrt{\pi} \approx 1.7724538509\) ์ ๋๋ค. ๊ฐ์ฐ์ค ์ ๋ถ๊ณผ ์ฐ๊ฒฐ๋ ์ ๋ช ํ ๊ฒฐ๊ณผ์ ๋๋ค.
๊ฒฐ๊ณผ๋ ์ผ๋ง๋ ์ ํํ๊ฐ์? Lanczos ๊ทผ์ฌ๋ ์ผ๋ฐ์ ์ธ ์ธ์์ ๋ํด ์ฝ 15์๋ฆฌ๊น์ง ์ ํํ๋ฉฐ, ๊ฑฐ์ ๋ชจ๋ ์ฉ๋์์ ์ฐจ๊ณ ๋์น๋ ์ ๋ฐ๋์ ๋๋ค.