ํด๋ฆฌ๊ฐ๋ง ํจ์๋?
์ฐจ์ n์ ํด๋ฆฌ๊ฐ๋ง ํจ์ \(\psi^{(n)}(a)\)๋ ๊ฐ๋ง ํจ์์ ์์ฐ๋ก๊ทธ๋ฅผ (n+1)๋ฒ ๋ฏธ๋ถํ ๊ฐ์ ๋๋ค. ์ฌ๊ธฐ์ ์ฌ์ฉํ๋ ์ธ๋ฑ์ฑ ๋ฐฉ์์์๋ \(n = -1\)์ด๋ฉด ๋ก๊ทธ๊ฐ๋ง ํจ์ \(\ln\Gamma(a)\) ์์ฒด๋ฅผ, \(n = 0\)์ด๋ฉด ๋๊ฐ๋ง ํจ์ \(\psi(a)\)๋ฅผ, \(n = 1\)์ด๋ฉด ํธ๋ฆฌ๊ฐ๋ง ํจ์๋ฅผ, \(n = 2\)์ด๋ฉด ํ ํธ๋ผ๊ฐ๋ง ํจ์๋ฅผ ๋ฐํํ๋ฉฐ ์ด๋ฐ ์์ผ๋ก ์ด์ด์ง๋๋ค. ์ด๋ฌํ ํน์ ํจ์๋ค์ ํ๋ฅ , ํต๊ณ(๊ฐ๋ง ๋ถํฌ์ ๋ฒ ํ ๋ถํฌ์ ์ต๋์ฐ๋์ถ์ ), ์ ์๋ก , ์ ๊ทผ ํด์ ์ ๋ฐ์ ๊ฑธ์ณ ๋ฑ์ฅํฉ๋๋ค. ์์ ์ํ์ ํด๋นํ๋ฏ๋ก ์ด๋์์๋ ๋์ผํ๊ฒ ์ ์ฉ๋ฉ๋๋ค.
๊ณ์ฐ๊ธฐ ์ฌ์ฉ ๋ฐฉ๋ฒ
๋๋กญ๋ค์ด์์ ๋ฏธ๋ถ ์ฐจ์ n(-1, 0, 1, 2, 3, 4 ์ค ํ๋)์ ์ ํํ๊ณ ์ธ์ a๋ฅผ ์ค์๋ก ์ ๋ ฅํ์ธ์. ๊ณ์ฐ ๋ฒํผ์ ๋๋ฅด๋ฉด \(\psi^{(n)}(a)\)๋ฅผ ์ฝ 12์๋ฆฌ ์ ํจ์ซ์๊น์ง ๊ตฌํ ์ ์์ต๋๋ค. ์ด ๊ณ์ฐ ๋ฐฉ์์ \(a > 0\)์์ ์์ ํ ์ ํจํ๋ฉฐ, a๊ฐ 0 ์ดํ์ธ ๊ฒฝ์ฐ์๋ ๋ฐ์ฌ ๊ณต์์ด ํ์ํฉ๋๋ค. ๋ํ 0 ์ดํ์ ์ ์(๊ฐ๋ง ํจ์์ ๊ทน์ )๋ ์ ์๋์ง ์์์ผ๋ก ํ์๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
\(n = -1\)์ผ ๋๋ \(\ln\Gamma\)๋ฅผ ์ํ ๋์ด์ค(Lanczos) ๊ทผ์ฌ๋ฅผ ์ฌ์ฉํด ๊ฐ์ ๊ณ์ฐํฉ๋๋ค. $$\psi^{(-1)}(a) = \ln\Gamma\!\left(a\right)$$ \(n = 0\)์ผ ๋ ๋๊ฐ๋ง๋ ์ ํ์ \(\psi(x) = \psi(x+1) - \frac{1}{x}\)๋ก ์ธ์๋ฅผ 6๋ณด๋ค ํฌ๊ฒ ๋ง๋ ๋ค, $$\psi^{(0)}(a) = \psi\!\left(a\right) = \frac{d}{da}\ln\Gamma\!\left(a\right) = \frac{\Gamma^{\prime}\!\left(a\right)}{\Gamma\!\left(a\right)}$$ ์ ๊ทผ ๊ธ์ \(\psi(x) \sim \ln x - \frac{1}{2x} - \frac{1}{12 x^2} + \cdots\) ๋ก ๊ตฌํฉ๋๋ค. \(n \geq 1\)์ผ ๋๋ ์ ํ์ \(\psi_n(x) = \psi_n(x+1) + (-1)^{n+1} \frac{n!}{x^{n+1}}\)๋ก ์ธ์๋ฅผ 10 ์ด์์ผ๋ก ์ด๋์ํจ ๋ค์, ๋ฒ ๋ฅด๋์ด ์ ๊ธฐ๋ฐ์ ์ ๊ทผ ๊ธ์๋ก ๊ณ์ฐ์ ๋ง๋ฌด๋ฆฌํฉ๋๋ค. $$\psi^{(n)}(a) = \frac{d^{n+1}}{da^{n+1}} \ln\Gamma\!\left(a\right), \quad n = \text{order}$$
๊ณ์ฐ ์์
\(n = 1\), \(a = 5\)(5์์์ ํธ๋ฆฌ๊ฐ๋ง)์ธ ๊ฒฝ์ฐ๋ฅผ ์ดํด๋ด ์๋ค. ํญ๋ฑ์ \(\psi_1(m) = \frac{\pi^2}{6} - \sum_{k=1}^{m-1} \frac{1}{k^2}\)๋ฅผ ์ฌ์ฉํ๋ฉด $$1.6449340668 - \left(1 + \frac{1}{4} + \frac{1}{9} + \frac{1}{16}\right) = 1.6449340668 - 1.4236111111 = 0.2213229557$$ ์ด ๋ฉ๋๋ค. ๊ณ์ฐ๊ธฐ๋ ์ฝ \(0.221322955737\)์ ๋ฐํํฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ a = 0์ด๊ฑฐ๋ ์์ ์ ์๋ฉด ์ ์๋์ง ์๋์? ๊ฐ๋ง ํจ์๋ 0, -1, -2, ...์์ ๊ทน์ ์ ๊ฐ์ง๋ฏ๋ก \(\ln\Gamma\)์ ๋ชจ๋ ํด๋ฆฌ๊ฐ๋ง ํจ์๊ฐ ๊ทธ ์ง์ ์์ ๋ฐ์ฐํฉ๋๋ค.
๋๊ฐ๋ง์ ํธ๋ฆฌ๊ฐ๋ง์ ์ฐจ์ด๋ ๋ฌด์์ธ๊ฐ์? ๋๊ฐ๋ง(\(n = 0\))๋ \(\ln\Gamma\)์ 1์ฐจ ๋ํจ์์ด๊ณ , ํธ๋ฆฌ๊ฐ๋ง(\(n = 1\))๋ 2์ฐจ ๋ํจ์๋ก ๋๊ฐ๋ง๋ฅผ ํ ๋ฒ ๋ ๋ฏธ๋ถํ ๊ฒ๊ณผ ๊ฐ์ต๋๋ค.
๊ฒฐ๊ณผ๋ ์ผ๋ง๋ ์ ํํ๊ฐ์? ์ด๋-์ ๊ทผ(shift-and-asymptotic) ๋ฐฉ๋ฒ์ ๋ฐฐ์ ๋ฐ๋(double-precision) ์ฐ์ฐ์ผ๋ก ์ ์ฉํ๋ฉด \(a > 0\)์์ ๋๋ต 12~15์๋ฆฌ์ ์ ํจ์ซ์๋ฅผ ์ ํํ๊ฒ ์ป์ ์ ์์ต๋๋ค.