์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ์ 1์ข ๊ตฌ๋ฉด ๋ฒ ์ ํจ์ \(j_v(x)\), ์ 2์ข ๊ตฌ๋ฉด ๋ฒ ์ ํจ์ \(y_v(x)\), ๊ทธ๋ฆฌ๊ณ ๊ฐ๊ฐ์ 1์ฐจ ๋ํจ์ \(j'_v(x)\)์ \(y'_v(x)\)๋ฅผ ๊ณ์ฐํฉ๋๋ค. ์ด ํจ์๋ค์ ๊ตฌ๋ฉด ์ขํ๊ณ์์ ํ๋ ๋ฐฉ์ ์๊ณผ ํฌ๋ฆํ์ธ ๋ฐฉ์ ์์ ๋๊ฒฝ(radial) ๋ฐฉํฅ ํด์ด๋ฉฐ, ์ฐ๋ ์ด๋ก , ์ ์๊ธฐํยท์ํฅ ๋ฐฉ์ฌ, ์์์ญํ(์์ ์ ์์ ๋ถ๋ถํ) ๋ฑ ๋ฌผ๋ฆฌํ ์ ๋ฐ์์ ๋ฑ์ฅํฉ๋๋ค. ์๋กญ๊ฒ ๊ตฌํ๋ ์ด ๋ฒ์ ์ ์ ์ ์ฐจ์ \(v \ge 0\)๊ณผ ์ค์ ์ธ์ \(x > 0\)์ ์ฒ๋ฆฌํฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ \(v\)(0, 1, 2์ ๊ฐ์ ์์ด ์๋ ์ ์)์ ์ธ์ \(x\)(์์ ์ค์)๋ฅผ ์ ๋ ฅํ์ธ์. ๊ณ์ฐ ๋ฒํผ์ ๋๋ฅด๋ฉด ๋ค ๊ฐ์ง ๊ฐ์ ๋ชจ๋ ์ป์ ์ ์์ต๋๋ค. ๋ค๋ง \(y_v(x)\)์ \(y'_v(x)\)๋ \(x\)๊ฐ 0์ ๊ฐ๊น์์ง์๋ก ๋ฐ์ฐํ๋ฏ๋ก, \(x = 0\)์์๋ ๋ฌดํ๋๋ก ํ์๋ฉ๋๋ค. \(j_0(0)\)์ ๊ทนํ๊ฐ์ผ๋ก 1์ด ๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
์ด ํจ์๋ค์ ๋ฏธ๋ถ๋ฐฉ์ ์ \(x^2 w'' + 2x w' + (x^2 - v(v+1))w = 0\)์ ๋ง์กฑํฉ๋๋ค. ๋ซํ ํํ์ธ \(j_0 = \sin(x)/x\)์ \(y_0 = -\cos(x)/x\)์์ ์ถ๋ฐํ์ฌ, ๋ ๋์ ์ฐจ์๋ ์ผํญ ์ ํ์ $$f_{n+1} = \frac{2n+1}{x}\,f_n - f_{n-1}$$์ ๋ฐ๋ฆ ๋๋ค. ์ํฅ ์ ํ๋ \(y_v\)์ ๋ํด์๋ ์์ ์ ์ด์ง๋ง, \(j_v\)์ ๊ฒฝ์ฐ \(n > x\)์ผ ๋ ๋ถ์์ ํฉ๋๋ค. ๊ทธ๋์ ๋ฐ๋ฌ(Miller)์ ํํฅ ์ ํ๋ฒ์ ์ฌ์ฉํฉ๋๋ค. ์ฆ ์ถฉ๋ถํ ๋์ ์ฐจ์์์ \(f\)๋ฅผ ๊ฐ๊ฐ 0๊ณผ 1๋ก ์ค์ ํด ์์ํ ๋ค ์ฐจ์๋ฅผ ๋ฎ์ถฐ ๊ฐ๋ฉฐ ์ ํ์ํค๊ณ , 0์ฐจ ํญ์ด \(\sin(x)/x\)์ ์ผ์นํ๋๋ก ๋ชจ๋ ๊ฐ์ ๋ค์ ์ค์ผ์ผ๋งํฉ๋๋ค. ๋ํจ์๋ \(j'_v = j_{v-1} - \frac{v+1}{x}j_v\) ๊ณต์์ ์ฌ์ฉํฉ๋๋ค.
๊ณ์ฐ ์์ (v = 0, x = 2)
$$j_0(2) = \frac{\sin(2)}{2} = 0.4546487134, \qquad y_0(2) = -\frac{\cos(2)}{2} = 0.2080734183$$์ ๋๋ค. \(j'_0 = -j_1\)์ด๋ฏ๋ก \(j'_0(2) = -0.4353977750\)์ด ๋๊ณ , \(y'_0 = -y_1\)์ด๋ฏ๋ก \(y'_0(2) = 0.3506120043\)์ด ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
๋ณต์์ \(x\)๋ ์ง์ํ๋์? ์๋์. ์๋ ํ์ด์ง๋ ๋ณต์์ ์ธ์๋ฅผ ๋ฐ์ง๋ง, ์ด ๋ฒ์ ์ ๋ช ํ์ฑ๊ณผ ์๋๋ฅผ ์ํด ์ค์ \(x > 0\)์ผ๋ก ์ ํํ์ต๋๋ค.
์ \(y_v\)๋ \(x = 0\)์์ ๋ฌดํ๋์ธ๊ฐ์? ์ 2์ข ํจ์๋ ์์ ์ ๊ทน(pole)์ ๊ฐ์ง๋ฏ๋ก, \(x\)๊ฐ 0์ ๊ฐ๊น์์ง์๋ก ๊ฐ์ด ํ์์ด ์ปค์ง๋๋ค.
์ ํ๋๋ ์ด๋ ์ ๋์ธ๊ฐ์? ๋ฐฐ์ ๋ฐ๋(double precision)๋ก ๊ณ์ฐํ์ฌ ์ฝ 15์๋ฆฌ์ ์ ํจ์ซ์๋ฅผ ์ ๊ณตํ๋ฏ๋ก, ์ผ๋ฐ์ ์ธ ๊ณตํยท๋ฌผ๋ฆฌ ์์ ์๋ ์ถฉ๋ถํฉ๋๋ค.