์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ์ 1์ข ๋ณํ ๋ฒ ์ ํจ์ \(I_v(x)\)์ ์ 2์ข ๋ณํ ๋ฒ ์ ํจ์ \(K_v(x)\), ๊ทธ๋ฆฌ๊ณ ์ด๋ค์ 1์ฐจ ๋ํจ์ \(I'_v(x)\)์ \(K'_v(x)\)๋ฅผ ๊ณ์ฐํฉ๋๋ค. ๋ ํจ์๋ ๋ณํ ๋ฒ ์ ๋ฐฉ์ ์ \(x^2 y'' + x y' - (x^2 + v^2)y = 0\)์ ์๋ก ๋ ๋ฆฝ์ ์ธ ๋ ํด์ ๋๋ค. ์ํต์์์ ์ด์ ๋, ํ์ฐ, ์ ์ก์ ยท๋ํ๊ด ์ด๋ก , ํต๊ณํ ๋ฑ ๋ฌผ๋ฆฌํ๊ณผ ๊ณตํ ์ ๋ฐ์ ๊ฑธ์ณ ๋ฑ์ฅํฉ๋๋ค. ์์ ์ํ์ด๋ฏ๋ก ๊ฒฐ๊ณผ๋ ์ด๋์๋ ๋์ผํ๋ฉฐ, ์ง์ญ๋ณ ๊ท์น์ ์ ์ฉ๋์ง ์์ต๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ v(์์์ ์ค์)์ ์ธ์ x๋ฅผ ์ ๋ ฅํ์ธ์. \(I_v(x)\)๋ v๊ฐ ์ ์์ผ ๋ ๋ชจ๋ ์ค์ x์ ๋ํด ๊ณ์ฐ๋๋ฉฐ, ๊ทธ ์ธ์ ๊ฒฝ์ฐ์๋ \(x \ge 0\)์์ ๊ณ์ฐ๋ฉ๋๋ค. \(K_v(x)\)๋ \(x \to 0^+\)์์ ๋ฐ์ฐํ๊ธฐ ๋๋ฌธ์ \(x > 0\)์ด์ด์ผ ํ๋ฉฐ, \(x \le 0\)์ธ ๊ฒฝ์ฐ์๋ ์ ์๋์ง ์์(NaN)์ผ๋ก ํ์๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
\(I_v(x)\)๋ ๊ฑฐ๋ญ์ ๊ณฑ ๊ธ์๋ก๋ถํฐ ํฉ์ฐํ๋ฉฐ, ๊ฐ๋ง ํจ์์๋ ๋์ด์ค(Lanczos) ๊ทผ์ฌ๋ฅผ ์ฌ์ฉํฉ๋๋ค.
$$I_v(x) = \sum_{k=0}^{\infty} \frac{1}{k!\,\Gamma(k + v + 1)} \left(\frac{x}{2}\right)^{2k + v}$$\(K_v(x)\)๋ $$K_v(x) = \int_{0}^{\infty} e^{-x\cdot\cosh t}\cdot\cosh(vt)\,dt$$๋ฅผ ์์น ์ ๋ถํ์ฌ ๊ตฌํ๋๋ฐ, ์ด ๋ฐฉ๋ฒ์ ์ ์ ์ฐจ์์ ๋น์ ์ ์ฐจ์ ๋ชจ๋์์ ์์ ์ ์ ๋๋ค. ๋ํจ์๋ ๋์นญ ์ ํ์ $$I'_v(x) = \tfrac{1}{2}\left(I_{v-1}(x) + I_{v+1}(x)\right)$$์ $$K'_v(x) = -\tfrac{1}{2}\left(K_{v-1}(x) + K_{v+1}(x)\right)$$๋ฅผ ์ฌ์ฉํ๋ฏ๋ก x๋ก ๋๋๋ ์ฐ์ฐ์ ํผํ ์ ์์ต๋๋ค.
๊ณ์ฐ ์์ (v = 0, x = 1)
๊ธ์๋ก๋ถํฐ $$I_0(1) = 1 + 0.25 + 0.015625 + \ldots \approx 1.26606588$$์ ์ป์ต๋๋ค. ์ ๋ถ์ผ๋ก๋ \(K_0(1) \approx 0.42102444\)๋ฅผ ์ป์ต๋๋ค. \(I_{-1} = I_1\)์ด๋ฏ๋ก ๋์นญ ํํ์์ \(I'_0(1) = I_1(1) \approx 0.56515910\), \(K'_0(1) = -K_1(1) \approx -0.60190723\)์ด ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
\(K_v(x)\)๊ฐ ์ ์ ์๋์ง ์์์ผ๋ก ํ์๋๋์? \(K_v(x)\)๋ \(x > 0\)์์๋ง ์ ์๋ฉ๋๋ค. 0 ์ดํ์์๋ ๋ฐ์ฐํ๊ธฐ ๋๋ฌธ์ ๋๋ค.
๋ถ์ ์ฐจ์๋ฅผ ์ฌ์ฉํ ์ ์๋์? ๋ค. ๋ ํจ์ ๋ชจ๋ ๋น์ ์์ ์์๋ฅผ ํฌํจํ ์์์ ์ค์ ์ฐจ์๋ฅผ ๋ฐ์ต๋๋ค.
์ ํ๋๋ ์ด๋ ์ ๋์ธ๊ฐ์? ๊ฒฐ๊ณผ๋ ๋ฐฐ์ ๋ฐ๋(์ ํจ์ซ์ ์ฝ 12~15์๋ฆฌ)๋ฅผ ์ฌ์ฉํ๋ฉฐ, ์ ๋นํ ๋ฒ์์ x์ ๋ํด ํ์ค ์ฐธ์กฐํ์ ์ผ์นํฉ๋๋ค.