์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ๊ณ ์ ๋ ์ค์ ์ฐจ์ v์ ๋ํด ์ 1์ข ๋ณํ ๋ฒ ์ ํจ์ \(I_{v}(x)\)๋ฅผ ์ผ๋ จ์ x ๊ฐ ๊ตฌ๊ฐ์์ ํ๋ก ๊ณ์ฐํฉ๋๋ค. ์ฐจ์, ์์ x ๊ฐ, ์ฆ๊ฐ๋(๊ฐ๊ฒฉ), ์์ฑํ ํ ์๋ฅผ ์ ๋ ฅํ๋ฉด \(x_{i} = \text{์์๊ฐ} + i\cdot\text{๊ฐ๊ฒฉ}\) ํํ์ ๋ชฉ๋ก์ ๋ง๋ค๊ณ ๊ฐ ์ง์ ์์ \(I_{v}(x_{i})\)๋ฅผ ๊ณ์ฐํ์ฌ ํ์ ๊ทธ๋ํ๋ฅผ ํจ๊ป ์ ๊ณตํฉ๋๋ค. ์์ ์ํ์ ํน์ํจ์ ๊ณ์ฐ ๋๊ตฌ์ด๋ฏ๋ก ์ง์ญ๋ณ ๊ท์น์ด๋ ๋จ์์ ๊ตฌ์ ๋ฐ์ง ์๊ณ ์ด๋์๋ ๋์ผํ๊ฒ ์ ์ฉ๋ฉ๋๋ค.
๊ณต์
๋ณํ ๋ฒ ์ ํจ์ \(I_{v}(x)\)๋ ๋ณํ ๋ฒ ์ ๋ฐฉ์ ์ \(x^{2}y'' + xy' - (x^{2} + v^{2})y = 0\)์ ํด์ ๋๋ค. ๋ณธ ๊ณ์ฐ๊ธฐ์์๋ ๋ค์ ๋ฉฑ๊ธ์๋ก ๊ฐ์ ๊ตฌํฉ๋๋ค.
$$I_{v}(x) = \sum_{k=0}^{\infty} \frac{1}{k!\;\Gamma(v+k+1)}\left(\frac{x}{2}\right)^{v+2k}$$๊ณ์น(ํฉํ ๋ฆฌ์ผ)๊ณผ ๊ฐ๋ง ํจ์๋ฅผ ์ฌ์ฉํ๋ฏ๋ก v์๋ ์์์ ์ค์๋ฅผ ๋ฃ์ ์ ์์ต๋๋ค. ์์น์ ์์ ์ฑ์ ์ํด ๊ฐ ํญ์ \(\ln\Gamma\)์ ๋์ด์ค(Lanczos) ๊ทผ์ฌ๋ฅผ ์ด์ฉํด ๋ก๊ทธ ๊ณต๊ฐ์์ ๊ณ์ฐํ ๋ค, ํญ์ ํฌ๊ธฐ๊ฐ ๋ฌด์ํ ๋งํผ ์์์ง ๋๊น์ง ํฉ์ฐํฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ v(์: 0, 1, 2.5), x์ ์ด๊ธฐ๊ฐ, ๊ฐ ํ๋ง๋ค x์ ๋ํด์ง๋ ์ฆ๊ฐ๋, ๊ทธ๋ฆฌ๊ณ ๋ฐ๋ณต ํ์(ํ ์)๋ฅผ ์ ๋ ฅํ์ธ์. ๊ณ์ฐ ๋ฒํผ์ ๋๋ฅด๋ฉด x์ \(I_{v}(x)\)๋ก ๊ตฌ์ฑ๋ ๋ ์ด์ง๋ฆฌ ํ์, ๊ฐ์ ๊ตฌ๊ฐ์ ๋ํ ๊ทธ๋ํ๊ฐ ๋ํ๋ฉ๋๋ค.
๊ณ์ฐ ์์
v = 0, ์์๊ฐ = 0, ๊ฐ๊ฒฉ = 0.5, ํ์ = 5๋ก ์ค์ ํ๋ฉด \(x = 0, 0.5, 1, 1.5, 2\)๊ฐ ๋๊ณ ๊ฒฐ๊ณผ๋ ๋ค์๊ณผ ๊ฐ์ต๋๋ค.
\(I_{0}(0) = 1\), \(I_{0}(0.5) \approx 1.0634834\), \(I_{0}(1) \approx 1.2660658\), \(I_{0}(1.5) \approx 1.6467232\), \(I_{0}(2) \approx 2.2795853\). ์ด ๊ฐ๋ค์ ํ์ค ์ฐธ์กฐํ์ ์ผ์นํฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ฐจ์๊ฐ ์์์ด๊ฑฐ๋ ์ ์๊ฐ ์๋์ด๋ ๋๋์? ๋ฉ๋๋ค. ์์ ์ ์ ์ฐจ์์๋ ํญ๋ฑ์ \(I_{-n}(x) = I_{n}(x)\)๊ฐ ์ ์ฉ๋ฉ๋๋ค. ์ ์๊ฐ ์๋ v๋ \(x \geq 0\)์์ ์ง์๋๋ฉฐ, \(x < 0\)์ด๋ฉด์ v๊ฐ ์ ์๊ฐ ์๋ ๊ฒฝ์ฐ์๋ ๊ฐ์ด ๋ณต์์๊ฐ ๋๋ฏ๋ก NaN์ ๋ฐํํฉ๋๋ค.
\(I_{v}(x)\)๋ ์ ์ด๋ ๊ฒ ๋น ๋ฅด๊ฒ ์ปค์ง๋์? ์ง๋ํ๋ ์ผ๋ฐ ๋ฒ ์ ํจ์ \(J_{v}\)์ ๋ฌ๋ฆฌ, ๋ณํ ๋ฒ ์ ํจ์๋ x๊ฐ ์ปค์ง์๋ก ๋๋ต \(e^{x}/\sqrt{2\pi x}\)์ ๋น๋กํด ์ฆ๊ฐํฉ๋๋ค. ๊ทธ๋์ x๊ฐ ์์ฃผ ํฌ๋ฉด ๊ฐ์ด ๋ฌดํ๋๋ก ๋ฐ์ฐ(overflow)ํ ์ ์์ต๋๋ค.
\(I_{v}(0)\)์ ๊ฐ์ ๋ฌด์์ธ๊ฐ์? \(I_{0}(0) = 1\)์ด๊ณ , \(v > 0\)์ธ ๊ฒฝ์ฐ \(I_{v}(0) = 0\)์ ๋๋ค.