MCP๋กœ ์—ฐ๊ฒฐ โ†’

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

์ฒซ ์ธ์ˆ˜์—์„œ์˜ ๊ฐ๋งˆ ๊ฐ’
1.772454
over 5 rows (5 finite)
์ƒ์„ฑ๋œ ํ–‰ ์ˆ˜ 5
์œ ํ•œํ•œ ๊ฐ’ 5
์œ ํ•œ๊ฐ’ ์ค‘ ์ตœ์†Ÿ๊ฐ’ 0.886227
์œ ํ•œ๊ฐ’ ์ค‘ ์ตœ๋Œ“๊ฐ’ 1.772454
a Gamma(a)
0.500000 1.772454
1.000000 1.000000
1.500000 0.886227
2.000000 1.000000
2.500000 1.329340

๊ฐ๋งˆ ํ•จ์ˆ˜ ๊ณ„์‚ฐ๊ธฐ๋ž€?

\(\Gamma(a)\)๋กœ ํ‘œ๊ธฐํ•˜๋Š” ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” ์ˆ˜ํ•™์—์„œ ๊ฐ€์žฅ ์ค‘์š”ํ•œ ํŠน์ˆ˜ ํ•จ์ˆ˜ ์ค‘ ํ•˜๋‚˜์ž…๋‹ˆ๋‹ค. ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” ํŒฉํ† ๋ฆฌ์–ผ์„ ํ•ด์„์ ์œผ๋กœ ํ™•์žฅํ•œ ๊ฒƒ์œผ๋กœ, ์Œ์ด ์•„๋‹Œ ์ •์ˆ˜ \(n\)์— ๋Œ€ํ•ด \(\Gamma(n+1) = n!\)์ด ์„ฑ๋ฆฝํ•ฉ๋‹ˆ๋‹ค. ํŒฉํ† ๋ฆฌ์–ผ๊ณผ ๋‹ฌ๋ฆฌ ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” ๋ถ„์ˆ˜์™€ ์Œ์ˆ˜๋ฅผ ํฌํ•จํ•œ ๊ฑฐ์˜ ๋ชจ๋“  ์‹ค์ˆ˜(๊ทธ๋ฆฌ๊ณ  ๋ณต์†Œ์ˆ˜)์—์„œ ์ •์˜๋ฉ๋‹ˆ๋‹ค. ์ด ๊ณ„์‚ฐ๊ธฐ๋Š” ๋“ฑ๊ฐ„๊ฒฉ์œผ๋กœ ๋Š˜์–ด์„  ์ธ์ˆ˜๋“ค์— ๋Œ€ํ•ด \(\Gamma(a)\) ๊ฐ’์˜ ํ‘œ๋ฅผ ๋งŒ๋“ค๊ณ , ๊ทธ ๊ฒฐ๊ณผ๋ฅผ ๊ณก์„  ๊ทธ๋ž˜ํ”„๋กœ ๊ทธ๋ ค ์ค๋‹ˆ๋‹ค. ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” ๊ตญ๊ฐ€๋ณ„ ๊ทœ์น™์ด ์—†๋Š” ๋ณดํŽธ์ ์ธ ์ˆ˜ํ•™์œผ๋กœ, ์–ด๋””์„œ๋‚˜ ๋™์ผํ•˜๊ฒŒ ์ ์šฉ๋ฉ๋‹ˆ๋‹ค.

์‚ฌ์šฉ ๋ฐฉ๋ฒ•

๋‹ค์Œ ์„ธ ๊ฐ€์ง€ ๊ฐ’์„ ์ž…๋ ฅํ•˜์„ธ์š”. a์˜ ์ดˆ๊นƒ๊ฐ’(์ฒซ ๋ฒˆ์งธ ์ธ์ˆ˜), ์ฆ๊ฐ€ํญ(step)(๊ฐ ํ–‰๋งˆ๋‹ค ๋”ํ•ด์ง€๋Š” ์ผ์ •ํ•œ ๊ฐ’), ๊ทธ๋ฆฌ๊ณ  ํ–‰ ๊ฐœ์ˆ˜์ž…๋‹ˆ๋‹ค. \(k\)๋ฒˆ์งธ ํ–‰์€ $$a_k = \text{์ดˆ๊นƒ๊ฐ’} + k \times \text{์ฆ๊ฐ€ํญ}$$ ์œผ๋กœ ์ธ์ˆ˜๋ฅผ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. ๊ณ„์‚ฐ๊ธฐ๋Š” ๊ฐ ์ธ์ˆ˜์—์„œ ๊ฐ๋งˆ ํ•จ์ˆ˜ ๊ฐ’์„ ๊ตฌํ•ด (์ธ์ˆ˜, ํ•จ์ˆซ๊ฐ’) ์Œ์„ ํ‘œ๋กœ ์ •๋ฆฌํ•˜๊ณ , ์œ ํ•œํ•œ ๊ฐ’ ์ค‘ ์ตœ์†Ÿ๊ฐ’๊ณผ ์ตœ๋Œ“๊ฐ’์„ ํ•จ๊ป˜ ๋ณด์—ฌ ์ค๋‹ˆ๋‹ค. \(a = 0, -1, -2, \dots\) ์ธ ๊ทน์ ์—์„œ๋Š” ๊ฐ’์ด ์ •์˜๋˜์ง€ ์•Š์Œ์œผ๋กœ ํ‘œ์‹œ๋ฉ๋‹ˆ๋‹ค.

๊ณต์‹ ์„ค๋ช…

๊ฐ๋งˆ ํ•จ์ˆ˜๋ฅผ ์ •์˜ํ•˜๋Š” ์ ๋ถ„์€ \(\operatorname{Re}(a) > 0\) ์ผ ๋•Œ $$\Gamma(a) = \int_{0}^{\infty} t^{\,a-1}\, e^{-t}\, dt$$ ์ž…๋‹ˆ๋‹ค. ์ˆ˜์น˜ ๊ณ„์‚ฐ์—์„œ๋Š” ์•ฝ 15์ž๋ฆฌ ์ •ํ™•๋„๋ฅผ ๋‚ด๋Š” ๋ž€์ดˆ์Šค ๊ทผ์‚ฌ(Lanczos approximation)๋ฅผ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค. \(a \le 0.5\) ์ธ ๊ฒฝ์šฐ์—๋Š” ๋ฐ˜์‚ฌ ๊ณต์‹(reflection formula) $$\Gamma(a) = \frac{\pi}{\sin(\pi \cdot a) \times \Gamma(1-a)}$$ ๋ฅผ ์ ์šฉํ•ด ์Œ์ˆ˜์™€ ์ž‘์€ ์ธ์ˆ˜๋ฅผ ์ฒ˜๋ฆฌํ•˜๊ณ  ๋ฐœ์‚ฐ์„ ํ”ผํ•ฉ๋‹ˆ๋‹ค. ๋Œ€ํ‘œ์ ์ธ ํŠน์ˆ˜๊ฐ’์œผ๋กœ๋Š” \(\Gamma(1/2) = \sqrt{\pi} = 1.77245\), \(\Gamma(1) = 1\), \(\Gamma(n+1) = n!\) ์ด ์žˆ์Šต๋‹ˆ๋‹ค.

์Œ์ด ์•„๋‹Œ ์ •์ˆ˜๊ฐ€ ์•„๋‹Œ ๊ณณ์— ๊ทน์ ์„ ๊ฐ€์ง„ ๊ฐ๋งˆ ํ•จ์ˆ˜ ๊ทธ๋ž˜ํ”„
๊ฐ๋งˆ ํ•จ์ˆ˜ ๊ณก์„ . ๊ธ‰๊ฒฉํ•œ ์ฆ๊ฐ€์™€ 0 ๋ฐ ์Œ์˜ ์ •์ˆ˜์—์„œ์˜ ์ˆ˜์ง ์ ๊ทผ์„ (๊ทน์ )์„ ๋ณด์—ฌ์ค€๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์‹œ

์ดˆ๊นƒ๊ฐ’ = 0.5, ์ฆ๊ฐ€ํญ = 0.5, ๊ฐœ์ˆ˜ = 5 ๋กœ ๋‘๋ฉด ์ธ์ˆ˜๋Š” 0.5, 1.0, 1.5, 2.0, 2.5 ๊ฐ€ ๋ฉ๋‹ˆ๋‹ค. ๊ทธ ๊ฒฐ๊ณผ๋Š” \(\Gamma(0.5) = 1.77245\) (\(= \sqrt{\pi}\)), \(\Gamma(1.0) = 1.0\), \(\Gamma(1.5) = 0.88623\), \(\Gamma(2.0) = 1.0\), \(\Gamma(2.5) = 1.32934\) ์ž…๋‹ˆ๋‹ค. ์–‘์ˆ˜ ๊ตฌ๊ฐ„์—์„œ ๊ณก์„ ์€ \(a = 1.4616\) ๋ถ€๊ทผ์—์„œ ์ตœ์†Ÿ๊ฐ’(์•ฝ 0.8856)๊นŒ์ง€ ๋‚ด๋ ค๊ฐ”๋‹ค๊ฐ€ ๋‹ค์‹œ ์œ„๋กœ ์˜ฌ๋ผ๊ฐ‘๋‹ˆ๋‹ค.

๊ฐ๋งˆ ์ ๋ถ„์„ ๋‚˜ํƒ€๋‚ด๋Š” ๊ณก์„  t^(a-1) e^(-t) ์•„๋ž˜์˜ ๋ฉด์ 
์ ๋ถ„ ์ •์˜์˜ ๊ธฐํ•˜ํ•™์  ์˜๋ฏธ: \(\Gamma(a)\)๋Š” \(t^{a-1} \cdot e^{-t}\) ์•„๋ž˜์˜ ์ƒ‰์น ๋œ ๋ฉด์ ์ด๋‹ค.

์ž์ฃผ ๋ฌป๋Š” ์งˆ๋ฌธ

์™œ \(\Gamma(0)\)์€ ์ •์˜๋˜์ง€ ์•Š๋‚˜์š”? 0, -1, -2, โ€ฆ ๊ฐ™์€ 0 ์ดํ•˜์˜ ์ •์ˆ˜๋Š” ๊ฐ๋งˆ ํ•จ์ˆ˜๊ฐ€ ์–‘์˜ ๋ฌดํ•œ๋Œ€ ๋˜๋Š” ์Œ์˜ ๋ฌดํ•œ๋Œ€๋กœ ๋ฐœ์‚ฐํ•˜๋Š” ๋‹จ์ˆœ ๊ทน์ (simple pole)์ž…๋‹ˆ๋‹ค. ๋”ฐ๋ผ์„œ ์ด ์ง€์ ์—๋Š” ์œ ํ•œํ•œ ๊ฐ’์ด ์กด์žฌํ•˜์ง€ ์•Š์Šต๋‹ˆ๋‹ค.

a๊ฐ€ ์Œ์ˆ˜์ผ ์ˆ˜๋„ ์žˆ๋‚˜์š”? ๋„ค. ์Œ์˜ ๋น„์ •์ˆ˜ ๊ฐ’์€ ๋ชจ๋‘ ์œ ํšจํ•ฉ๋‹ˆ๋‹ค. ์—ฐ์†ํ•œ ์Œ์˜ ์ •์ˆ˜ ์‚ฌ์ด์—์„œ ๊ฐ’์˜ ๋ถ€ํ˜ธ๊ฐ€ ๋ฒˆ๊ฐˆ์•„ ๋ฐ”๋€Œ๋ฉฐ ์ ˆ๋Œ“๊ฐ’์ด ํฌ๊ฒŒ ์ปค์ง‘๋‹ˆ๋‹ค. ์˜ˆ๋ฅผ ๋“ค์–ด \(\Gamma(-0.5) = -3.5449\) ์ž…๋‹ˆ๋‹ค.

๊ฒฐ๊ณผ๋Š” ์–ผ๋งˆ๋‚˜ ์ •ํ™•ํ•œ๊ฐ€์š”? ๋ž€์ดˆ์Šค ๊ทผ์‚ฌ๋Š” ์•ฝ 15์ž๋ฆฌ์˜ ์œ ํšจ์ˆซ์ž๋ฅผ ์ œ๊ณตํ•˜๋ฏ€๋กœ, ์‚ฌ์‹ค์ƒ ๋ชจ๋“  ์‹ค๋ฌดยท๊ต์œก์šฉ ๋ชฉ์ ์— ์ถฉ๋ถ„ํ•ฉ๋‹ˆ๋‹ค.

์ตœ์ข… ์—…๋ฐ์ดํŠธ: