์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ \(f(x) = b_0 + \cfrac{a_1}{b_1 + \cfrac{a_2}{b_2 + \cfrac{a_3}{b_3 + \cdots}}}\) ํํ์ ์ผ๋ฐํ ์ฐ๋ถ์(generalized continued fraction)๋ฅผ ๊ณ์ฐํฉ๋๋ค. ๊ฐ์ฅ ํฐ ํน์ง์ ์ฒซ ํญ \(b_0\), n๋ฒ์งธ ๋ถ์ \(a_n\), n๋ฒ์งธ ๋ถ๋ชจ \(b_n\)์ ๋ณ์ \(x\)์ ํญ ๋ฒํธ \(n\)์ ์์กดํ๋ ์์์ผ๋ก ์ง์ ์ ๋ ฅํ ์ ์๋ค๋ ์ ์ ๋๋ค. ์ ๋ ฅํ \(x\) ๊ฐ์ ๋์ ํ ๋ค \(n = 1, 2, 3, \ldots\) ์์๋ก ๊ฐ ํญ์ ์์ฑํ๊ณ , ์๋ ด๊ฐ \(f(x)\)์ ํจ๊ป ๋ถ๋ถ ๊ทผ์ฌ๊ฐ \(f_n(x)\)์ ํ๋ก ์ ๋ฆฌํด ๋ณด์ฌ์ค๋๋ค. ์์ํ ์์นํด์ ๋๊ตฌ์ด๋ฏ๋ก ํน์ ๊ตญ๊ฐ์ ๊ท์ ์ด๋ ๋จ์์๋ ๋ฌด๊ดํ๊ฒ ์ด๋์๋ ๋์ผํ๊ฒ ์ ์ฉ๋ฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ธ ๊ฐ์ ์์๊ณผ \(x\) ๊ฐ์ ์ ๋ ฅํ์ธ์. ๊ฐ ์์์๋ ๊ธฐํธ \(x\)์ \(n\), ์ฐ์ฐ์ + - * / ^(๊ฑฐ๋ญ์ ๊ณฑ), ํจ์ sqrt, exp, ln, log, sin, cos, tan, ๊ทธ๋ฆฌ๊ณ ์์ pi์ e๋ฅผ ์ฌ์ฉํ ์ ์์ต๋๋ค. ์ ํจ์ซ์ ์๋ฆฟ์๋ฅผ ์ ํํ๋ฉด ๊ฒฐ๊ณผ ํ์ ๋ฐฉ์์ ์กฐ์ ํ ์ ์๋๋ฐ, ์ด๋ ์๋ ด ํ์ฉ์ค์ฐจ์ ํ์ ํ์์๋ง ์ํฅ์ ์ค ๋ฟ ์ค์ ๊ณ์ฐ ์์ฒด๋ ๋ฐ๋์ง ์์ต๋๋ค. ๊ฒฐ๊ณผ ์ฐฝ์๋ \(f(x)\) ๊ฐ๊ณผ ์๋ ด์ด ์ด๋ฃจ์ด์ง ํญ ๋ฒํธ \(n\)์ด ํ์๋๋ฉฐ, ํ์์๋ ์ด๊ธฐ ๊ทผ์ฌ๊ฐ๋ค์ด ์ ์ฐจ ํ ๊ฐ์ผ๋ก ๋ชจ์ฌ๋๋ ๊ณผ์ ์ ํ์ธํ ์ ์์ต๋๋ค.
๊ณต์ ํ์ด
๊ทผ์ฌ๊ฐ์ ์ ์ง ๊ธฐ๋ณธ ์ ํ์(forward fundamental recurrence)์ผ๋ก ๊ณ์ฐ๋ฉ๋๋ค. \(A_{-1} = 1\), \(A_0 = b_0\), \(B_{-1} = 0\), \(B_0 = 1\)์์ ์ถ๋ฐํ์ฌ, ๊ฐ ๋จ๊ณ๋ง๋ค $$A_n = b_n A_{n-1} + a_n A_{n-2}, \quad B_n = b_n B_{n-1} + a_n B_{n-2}$$๋ก ๊ฐฑ์ ํ๊ณ , n๋ฒ์งธ ๊ทผ์ฌ๊ฐ์ \(f_n = A_n / B_n\)์ผ๋ก ๊ตฌํฉ๋๋ค. ์ฐ์ํ ๋ ๊ทผ์ฌ๊ฐ์ด ์์ฒญํ ์ ๋ฐ๋๊น์ง ์ผ์นํ๊ฑฐ๋, ํญ์ด 1000๊ฐ์ ๋๋ฌํ๋ฉด ๋ฐ๋ณต์ ๋ฉ์ถฅ๋๋ค.
์์ ํ์ด
๊ธฐ๋ณธ๊ฐ์ธ \(b_0 = 1\), \(a_n = x - 1\), \(b_n = 2\)๋ฅผ ๊ทธ๋๋ก ๋๊ณ \(x = 5\)๋ก ๊ณ์ฐํด ๋ด ์๋ค. ์ด๊ฒ์ \(\sqrt{x}\)์ ๋ํ ๊ณ ์ ์ ์ธ ์ฐ๋ถ์์ ๋๋ค: $$\sqrt{x} = 1 + \cfrac{x-1}{2 + \cfrac{x-1}{2 + \cdots}}.$$ \(x = 5\)์ผ ๋ ๋ชจ๋ ๋ถ์๋ 4, ๋ชจ๋ ๋ถ๋ชจ๋ 2๊ฐ ๋ฉ๋๋ค. ์ด๊ธฐ ๊ทผ์ฌ๊ฐ์ \(f_1 = 3\), \(f_2 = 2\), \(f_3 = 2.3333\ldots\), \(f_4 = 2.2\)์ด๋ฉฐ, ๋ชจ๋ \(\sqrt{5} = 2.2360679774997896\)์ผ๋ก ์๋ ดํฉ๋๋ค. ๊ณ ์ ์ \(t = 2 + \frac{4}{t}\)๋ \(t = 1 + \sqrt{5}\)๋ฅผ ๋ง์กฑํ๋ฏ๋ก \(f = 1 + \frac{4}{1+\sqrt{5}} = \sqrt{5}\)๊ฐ ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
\(a_n\)๊ณผ \(b_n\)์ด ํญ ๋ฒํธ \(n\)์ ์์กดํด๋ ๋๋์? ๋ฉ๋๋ค. ์๋ฅผ ๋ค์ด ํจ์ \(\frac{x}{e^x - 1}\)์ \(b_0 = 1 - \frac{x}{2}\), \(a_n = \frac{x^2}{4}\), \(b_n = 2n + 1\)๋ก ํํ๋๋ฉฐ, ์ฌ๊ธฐ์ \(b_n\)์ \(n\)์ด ์ปค์ง์๋ก ํจ๊ป ์ฆ๊ฐํฉ๋๋ค.
๋ถ๋ชจ๊ฐ 0์ด ๋๋ฉด ์ด๋ป๊ฒ ๋๋์? ๊ณ์ฐ๊ธฐ๋ ์์ฃผ ์์ ์ ์ค๋ก (epsilon)์ ๋์ ํด ๊ณ์ฐ์ ์ด์ด๊ฐ๋ฉฐ, ์ด๋ ์์ ๋ ์ธ (Lentz) ์๊ณ ๋ฆฌ์ฆ๊ณผ ๋์ผํ ๋ฐฉ์์ ๋๋ค. ๋๊น์ง ์๋ ดํ์ง ์์ผ๋ฉด ๋ณ๋๋ก ์๋ ค ๋๋ฆฝ๋๋ค.
์ 1000๊ฐ ํญ์์ ๋ฉ์ถ๋์? ์์ ์ ์ํ ์ํ๊ฐ์ ๋๋ค. ๊ทธ๋๊น์ง๋ ์๋ ดํ์ง ์์ผ๋ฉด ๋ง์ง๋ง ๊ทผ์ฌ๊ฐ์ ๊ฒฝ๊ณ ์ ํจ๊ป ๋ฐํํฉ๋๋ค.