A fixed point iteration method for the arctangent with any odd order of convergence based on sine and cosine

तकनीकी सारांश: किसी भी विषम क्रम के अभिसरण (Order of Convergence) के साथ आर्कटैन्जेंट (Arctangent) के लिए एक फिक्स्ड पॉइंट इटरेशन विधि

समस्या विवरण
यह शोध पत्र $t \in \mathbb{R}^+$ के लिए आर्कटैन्जेंट फलन, $\arctan(t)$, की गणना करने की कम्प्यूटेशनल चुनौती को संबोधित करता है। जबकि $\arctan(u)$ के लिए मानक पावर सीरीज़ विस्तार $|u| < 1$ के लिए अभिसरित होता है, उच्च परिशुद्धता (precision) प्राप्त करने के लिए अक्सर कई पदों या विशिष्ट रूपांतरणों की आवश्यकता होती है। लेखक एक फिक्स्ड-पॉइंट इटरेशन विधि की खोज करते हैं जो किसी भी विषम क्रम के अभिसरण को प्राप्त कर सकती है, जिससे परिणाम के सही अंकों के प्राप्त होने की दर में तेजी आती है।

कार्यप्रणाली (Methodology)
प्रस्तावित विधि का मूल आधार आर्कटैन्जेंट के ट्रंकेटेड पावर सीरीज़ से व्युत्पन्न एक फिक्स्ड-पॉइंट फलन $T(x)$ है। लेखक त्रिकोणमितीय पहचानों के आधार पर एक प्रतिस्थापन चर (substitution variable) $u$ को परिभाषित करते हैं:
$$u = \frac{\sin(x) - t \cos(x)}{\cos(x) + t \sin(x)}$$
इस $u$ को $\arctan(u)$ की प्रथम $P$ पदों की पावर सीरीज़ में प्रतिस्थापित करके, लेखक इटरेशन फंक्शन का निर्माण करते हैं:
$$T(x) = x - \sum_{k=1}^{P} \frac{(-1)^{k-1}}{2k-1} \left( \frac{\sin(x) - t \cos(x)}{\cos(x) + t \sin(x)} \right)^{2k-1}$$
इटरेशन $x_{n+1} = T(x_n)$ के माध्यम से आगे बढ़ता है।

शोध पत्र इस फलन के गुणों के लिए एक कठोर गणितीय प्रमाण प्रदान करता है:

1. फिक्स्ड पॉइंट गुण (Fixed Point Property): यह सिद्ध किया गया है कि $T(\arctan(t)) = \arctan(t)$। यह इस तथ्य पर निर्भर करता है कि जब $x = \arctan(t)$ होता है, तो प्रतिस्थापन पद शून्य हो जाता है।
2. डेरिवेटिव विश्लेषण (Derivative Analysis): लेखक प्रदर्शित करते हैं कि फिक्स्ड पॉइंट पर, $T(x)$ के प्रथम $2P$ डेरिवेटिव शून्य होते हैं ($T^{(k)}(\arctan(t)) = 0$ जहाँ $0 \le k \le 2P$)। $(2P+1)$-वाँ डेरिवेटिव गैर-शून्य है और इसका मान $(-1)^P (2P)!$ है।
3. अभिसरण प्रमाण (Convergence Proof): बानाच फिक्स्ड पॉइंट थ्योरम (Banach Fixed Point Theorem) का उपयोग करते हुए, शोध पत्र स्थापित करता है कि $\arctan(t)$ के पर्याप्त निकट प्रारंभिक मान $x_0$ के लिए, अनुक्रम (sequence) अभिसरित होता है।
4. अभिसरण का क्रम (Order of Convergence): शून्य डेरिवेटिव्स के आधार पर, त्रुटि पद (error term) का टेलर विस्तार यह दर्शाता है कि विधि का अभिसरण क्रम ठीक $2P + 1$ है।

मुख्य योगदान और परिणाम

• अनिश्चित विषम क्रम (Arbitrary Odd Order): प्राथमिक योगदान एक एकीकृत ढांचा है जो उपयोगकर्ता को एक विशिष्ट विषम क्रम ($3, 5, 7, \dots$) प्राप्त करने के लिए किसी भी प्राकृतिक संख्या $P$ का चयन करने की अनुमति देता है।
• $\pi$ के लिए अनुप्रयोग: इस विधि को $\pi/4$ के लिए मैचिन-जैसे (Machin-like) सूत्रों, विशेष रूप से $\pi/4 = \arctan(1)$ पर लागू किया गया है। शोध पत्र विवरण देता है कि कैसे तर्क $t = 1/q$ के लिए इटरेशन फॉर्मूला सरल हो जाता है।
• उच्च-परिशुद्धता गणना (High-Precision Computation): शोध पत्र $\pi/4$ की एक मिलियन अंकों की गणना करने के लिए एक व्यावहारिक कम्प्यूटेशनल प्रयोग प्रस्तुत करता है।
  • कॉन्फ़िगरेशन: प्रयोग ने $P=2$ का उपयोग किया, जिससे अभिसरण क्रम 5 (quintic convergence) प्राप्त हुआ।
  • रणनीति: प्रदर्शन को अनुकूलित करने के लिए, लेखक ने एक "स्टेप-होल्डिंग" परिशुद्धता रणनीति का उपयोग किया। सभी इटरेशन को पूर्ण परिशुद्धता के साथ गणना करने के बजाय, प्रत्येक चरण में महत्वपूर्ण अंकों की संख्या को ज्यामितीय रूप से ($14 \times 5^n$) बढ़ाया गया, जो सही अंकों की अपेक्षित वृद्धि से मेल खाता है।
  • परिणाम: 14-अंकों के सन्निकटन (approximation) से शुरू होकर, विधि ने केवल 8 इटरेशन में दस लाख से अधिक अंकों की सटीकता प्राप्त की। प्रयोग में देखे गए त्रुटि न्यूनीकरण कारकों (error reduction factors) ने अभिसरण के सैद्धांतिक क्रम की पुष्टि की।

महत्व और दावे
शोध पत्र का दावा है कि यह विधि आर्कटैन्जेंट और, इसके विस्तार के रूप में, $\pi$ की उच्च-परिशुद्धता गणना के लिए एक कुशल विकल्प प्रदान करती है। इसका महत्व इस बात में निहित है कि एल्गोरिदम की मौलिक संरचना को बदले बिना, पैरामीटर $P$ के माध्यम से अभिसरण क्रम को ट्यून किया जा सकता है। लेखक उल्लेख करते हैं कि यह विधि स्व-सुधारात्मक (self-correcting) है और जब इसे तेज़ त्रिकोणमितीय फलन मूल्यांकन (जैसे टेलर सीरीज़ या अन्य एल्गोरिदम के माध्यम से) के साथ जोड़ा जाता है, तो यह पर्याप्त कम्प्यूटेशनल संसाधनों के दिए जाने पर $\arctan(t)$ की "असीमित परिशुद्धता" के साथ गणना करने की अनुमति देती है। यह कार्य दक्षता के प्रदर्शन के रूप में प्रस्तुत किया गया है, जिसे मानक हार्डवेयर और कंप्यूटर अलजेब्रा सॉफ्टवेयर (Maple) का उपयोग करके दस लाख से अधिक दशमलव स्थानों तक $\pi/4$ की सफल गणना द्वारा सत्यापित किया गया है।