다항식만으로 구현하는 임의 차수 수렴 뿌리 찾기

초록

본 논문은 양의 실수 a와 자연수 M에 대해 t⁽ᴹ⁾‑a=0 의 실근인 a¹ᐟᴹ을, 적절히 선택한 정수 P에 따라 차수(P+1)의 수렴을 보장하는 다항식 고정점 반복법으로 계산하는 방법을 제시한다. 고정점 함수 F(x)를 적분과 이항정리를 이용해 전개하고, 그 도함수들의 특성을 분석해 Banach 고정점 정리를 적용함으로써 수렴성과 수렴 차수를 증명한다. 또한 P=1‒4에 대한 구체적 다항식 형태와 대규모 계산 예시(√2의 백만 자리) 등을 통해 실용성을 확인한다.

상세 분석

이 논문은 전통적인 뉴턴 방법이 2차 수렴을 보이는 데 반해, 사용자가 원하는 임의의 정수 P에 대해 (P+1) 차수의 수렴을 달성할 수 있는 새로운 고정점 반복법을 제안한다. 핵심 아이디어는 f(t)=1−tᴹ/a 라는 단순 다항식의 P제곱을 적분하여 얻은 함수

 F(x)=∏{ℓ=1}^{P}!\left(1+\frac{1}{ℓM}\right)\int{0}^{x}!\bigl(1-\frac{t^{M}}{a}\bigr)^{P}dt

를 고정점 함수로 삼는 것이다. 베타 함수와 감마 함수를 이용해 적분을 정확히 계산하면, F(x)는

 F(x)=\sum_{k=0}^{P}(-1)^{k}\binom{P}{k}\frac{x^{kM+1}}{a^{k}(kM+1)}\prod_{ℓ=1}^{P}!\left(1+\frac{1}{ℓM}\right)

라는 형태의 다항식이 된다. 여기서 중요한 점은 계수들이 모두 사전에 계산 가능하며, 차수가 P·M+1 로 제한된다는 것이다.

정리 1에서는 (a) a¹ᐟᴹ이 F의 고정점임을 적분 변수를 u=tᴹ/a 로 바꾸어 베타 적분 형태로 변형하고, 감마 함수의 재귀 관계를 이용해 ∏{ℓ=1}^{P}(1+ℓ⁻¹M)·a¹ᐟᴹ·a⁻¹ᐟᴹ= a¹ᐟᴹ임을 보인다. (b)에서는 1차 도함수 F’(x)=∏{ℓ=1}^{P}(1+ℓ⁻¹M)(1−xᴹ/a)^{P} 를 구하고, x=a¹ᐟᴹ에서 0이 되는 것을 확인한다. (c)에서는 고차 도함수들을 전개하여 k≤P에서는 모두 0, k=P+1에서는

 F^{(P+1)}(a¹ᐟᴹ)=(-1)^{P}·M·a¹ᐟᴹ·∏_{ℓ=1}^{P}(1+ℓM)

임을 증명한다.

Banach 고정점 정리를 적용하면, F’(a¹ᐟᴹ)=0이므로 근처에서 Lipschitz 상수 L<1을 만족하는 구간이 존재하고, 따라서 초기값 x₀가 충분히 근접하면 x_{n+1}=F(x_n) 은 전역적으로 a¹ᐟᴹ으로 수렴한다.

수렴 차수는 Taylor 전개를 이용해 F(x)=a¹ᐟᴹ+F^{(P+1)}(a¹ᐟᴹ)/(P+1)!·(x−a¹ᐟᴹ)^{P+1}+o((x−a¹ᐟᴹ)^{P+1}) 로 나타내어, 오류가 (x_n−a¹ᐟᴹ)^{P+1}에 비례함을 보인다. 따라서 정확히 (P+1) 차수의 수렴을 얻는다.

실제 구현에서는 계수 c_k=∏_{ℓ=1}^{P}(1+ℓ⁻¹M)(-1)^{k}a^{-k}\binom{P}{k}/(kM+1) 를 미리 계산하고, Horner 방식 등 효율적인 다항식 평가 알고리즘을 사용한다. 논문은 P=1,2,3,4에 대한 구체적인 다항식 형태를 제시하고, 3√10, √2(백만 자리) 등 두 가지 사례를 통해 이론적 수렴 차수와 실제 연산 시간·정밀도 사이의 일치를 실험적으로 확인한다. 전체적으로, 이 방법은 고차 수렴을 필요로 하는 고정밀 계산에 적합하며, 복소수 확장이나 다른 비선형 방정식에도 적용 가능성을 시사한다.