A fixed-point iteration method for the number Pi with arbitrary odd order of convergence based on the sine function

यहाँ एलोइस शीसल (Alois Schiessl) के शोध पत्र "A fixed-point iteration method for the number Pi with arbitrary odd order of convergence based on the sine function" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र उच्च परिशुद्धता (high precision) के साथ गणितीय स्थिरांक $\pi$ की गणना करने की चुनौती को संबोधित करता है। जबकि कई एल्गोरिदम मौजूद हैं (जैसे, Chudnovsky algorithm, Machin-like formulas), यह कार्य एक नवीन फिक्स्ड-पॉइंट इटरेशन विधि प्रस्तावित करता है जो एक स्वेच्छ विषम क्रम की अभिसरण (arbitrary odd order of convergence) प्राप्त करने में सक्षम है। इसका लक्ष्य साइन (sine) फलन पर आधारित एक सरल पुनरावृत्ति योजना (iterative scheme) का निर्माण करना है, जो पर्याप्त निकट प्रारंभिक अनुमान (initial guess) उपलब्ध होने पर, जटिल श्रेणी विस्तारों (series expansions) या प्रत्येक चरण में मॉड्यूलर अंकगणित की आवश्यकता के बिना $\pi$ की ओर तीव्र अभिसरण की अनुमति देता है।

2. कार्यप्रणाली (Methodology)

मुख्य फलन (The Core Function)

लेखक एक पूर्णांक $P \in \mathbb{N}$ द्वारा पैरामीटराइज्ड एक फिक्स्ड-पॉइंट फलन $S(x)$ को परिभाषित करते हैं। इस फलन का निर्माण $\arcsin(t)$ के पावर सीरीज़ विस्तार को ट्रंकेट (truncate) करके और $t = \sin(x)$ को प्रतिस्थापित करके किया गया है।

फलन को इस प्रकार परिभाषित किया गया है:
$$S(x) = x + \sum_{k=1}^{P} \left( \prod_{\ell=1}^{k-1} \frac{2\ell-1}{2\ell} \right) \frac{\sin(x)^{2k-1}}{2k-1}$$

इटरेशन नियम है:
$$x_{n+1} = S(x_n)$$

सैद्धांतिक आधार (Theoretical Foundation)

यह विधि $\arcsin(t)$ के टेलर सीरीज़ (Taylor series) पर निर्भर करती है:
$$\arcsin(t) = \sum_{k=1}^{\infty} \left( \prod_{\ell=1}^{k-1} \frac{2\ell-1}{2\ell} \right) \frac{t^{2k-1}}{2k-1}$$
$t = \sin(x)$ को प्रतिस्थापित करने पर, अनंत श्रेणी सैद्धांतिक रूप से $x$ के बराबर होती है (जब $x \approx \pi$ हो जहाँ $\sin(x)$ छोटा होता है)। यह शोध पत्र इस श्रेणी को $P$ पदों तक ट्रंकेट करके $\sin(x)$ के एक परिमित बहुपद (finite polynomial) को $x$ में जोड़ने का कार्य करता है।

गणितीय प्रमाण (Mathematical Proof)

शोध पत्र तीन प्रमुख गुणों का कठोर प्रमाण प्रदान करता है:

1. फिक्स्ड पॉइंट (Fixed Point): $S(\pi) = \pi$ क्योंकि $\sin(\pi) = 0$, जिससे सभी योग पद शून्य हो जाते हैं।
2. लुप्त अवकलज (Vanishing Derivatives): लेखक सिद्ध करते हैं कि $x=\pi$ पर $S(x)$ के पहले $2P$ अवकलज (derivatives) शून्य हैं ($S^{(k)}(\pi) = 0$ जब $1 \le k \le 2P$)। यह $\sin(x)^{2P-1}$ पद पर सावधानीपूर्वक अवकलन और लीबनिज नियम (Leibniz rule) को लागू करके प्राप्त किया गया है।
3. अभिसरण का क्रम (Convergence Order): $\pi$ पर $(2P+1)$-वाँ अवकलज गैर-शून्य है। विशेष रूप से:
   $$S^{(2P+1)}(\pi) = \left( \prod_{\ell=1}^{P} (2\ell-1) \right)^2$$
   फिक्स्ड-पॉइंट इटरेशन के सिद्धांत के अनुसार, यदि फिक्स्ड पॉइंट पर पहले $m-1$ अवकलज शून्य हैं और $m$-वाँ अवकलज गैर-शून्य है, तो अभिसरण का क्रम $m$ होता है। अतः, अभिसरण का क्रम सटीक रूप से $2P + 1$ है।

अनुकूलन रणनीतियाँ (Optimization Strategies)

व्यावहारिक उच्च-परिशुद्धता गणना के लिए, लेखक तीन अनुकूलन का सुझाव देते हैं:

1. उच्च-परिशुद्धता प्रारंभिक अनुमान (High-Precision Initial Guess): $x_0$ से शुरुआत करें जो पहले से ही $\approx 18$ अंकों तक सटीक हो।
2. परिवर्तनीय परिशुद्धता अंकगणित (Variable Precision Arithmetic): चूंकि यह विधि स्व-सुधारात्मक (self-correcting) है, इसलिए शुरुआती चरणों में पूर्ण परिशुद्धता की आवश्यकता नहीं होती है। परिणाम की बढ़ती सटीकता से मेल खाने के लिए कार्यशील परिशुद्धता (अंकों की संख्या) प्रत्येक चरण में ज्यामितीय रूप से ($18 \times (2P+1)^n$) बढ़ाई जाती है, जिससे गणना का समय बचता है।
3. गुणांक पूर्व-गणना (Coefficient Pre-calculation): बहुपद के गुणांक स्थिर होते हैं और उन्हें पूर्व-गणना करके संग्रहीत किया जा सकता है।

3. प्रमुख योगदान (Key Contributions)

• स्वेच्छ विषम अभिसरण (Arbitrary Odd Convergence): यह शोध पत्र एक एकीकृत ढांचा प्रस्तुत करता है जहाँ पैरामीटर $P$ को बढ़ाकर अभिसरण क्रम को किसी भी विषम पूर्णांक ($3, 5, 7, \dots$) में ट्यून किया जा सकता है।
• कार्यान्वयन की सरलता (Simplicity of Implementation): अन्य उच्च-क्रम विधियों के विपरीत, जिन्हें जटिल तर्क अनुमानों (rational approximations) या कई त्रिकोणमितीय मूल्यांकनों की आवश्यकता होती है, यह विधि केवल साइन फलन और सरल बहुपद अंकगणित पर निर्भर करती है।
• कठोर व्युत्पत्ति (Rigorous Derivation): शोध पत्र टेलर विस्तार और लीबनिज नियम का उपयोग करके अभिसरण क्रम का पूर्ण विश्लेषणात्मक प्रमाण प्रदान करता है, जो इसे केवल अनुमानित संख्यात्मक विधियों से अलग करता है।
• व्यावहारिक दक्षता (Practical Efficiency): $P=4$ (9वें क्रम का अभिसरण) का उपयोग करके केवल 6 इटरेशन में $\pi$ के 10 लाख अंक प्राप्त करने का प्रदर्शन इस विधि की अत्यधिक दक्षता को उजागर करता है।

4. परिणाम (Results)

लेखक ने एक मानक उपभोक्ता पीसी (AMD Ryzen 5 5600G) पर Maple 2025.2 का उपयोग करके कम्प्यूटेशनल परीक्षण किए।

• मामला 1 ($P=1$, 3rd Order):

  • $x_0 = 3$ से शुरुआत।
  • 7 इटरेशन के बाद, परिणाम ने >900 सही दशमलव अंक प्राप्त किए।
  • त्रुटि में कमी (error reduction) क्यूबिक पैटर्न ($10^{-11} \to 10^{-34} \to 10^{-100}$) का पालन करती है।

• मामला 2 ($P=4$, 9th Order):

  • लक्ष्य: 1,000,000 अंक।
  • $x_0 \approx 3.14159...$ (18 अंक) से शुरुआत।
  • परिवर्तनीय परिशुद्धता का उपयोग किया (162 अंकों से शुरू होकर, स्केल करते हुए)।
  • परिणाम: केवल 6 इटरेशन में 1,059,070 से अधिक सही अंक प्राप्त किए।
  • त्रुटि घातांक (error exponent) प्रत्येक चरण में लगभग 9 के कारक से बढ़ा, जो सैद्धांतिक 9वें क्रम के अभिसरण की पुष्टि करता है।

5. महत्व (Significance)

यह शोध पत्र उच्च-परिशुद्धता $\pi$ गणना के लिए एक गणितीय रूप से सुंदर और कम्प्यूटेशनल रूप से कुशल विधि प्रस्तुत करता है।

• सैद्धांतिक मूल्य: यह दर्शाता है कि कैसे $\arcsin$ श्रृंखला को ट्रंकेट करने से एक शक्तिशाली फिक्स्ड-पॉइंट इटरेशन बनता है जिसका अभिसरण क्रम नियंत्रित किया जा सकता है।
• व्यावहारिक मूल्य: यह विधि किसी भी उच्च-परिशुद्धता अंकगणित लाइब्रेरी के लिए अत्यधिक उपयुक्त है। "डबलिंग" (या इस मामले में, $(2P+1)$ से गुणा करने वाली) परिशुद्धता रणनीतियों के साथ एक उच्च-क्रम योजना को जोड़कर, यह महंगे उच्च-परिशुद्धता साइन मूल्यांकन की लागत को कम करती है।
• सुलभता: अन्य अल्ट्रा-फास्ट एल्गोरिदम (जैसे Chudnovsky algorithm जिसमें जटिल मॉड्यूलर एक्सपोनेंटिएशन शामिल है) की तुलना में इस एल्गोरिदम को लागू करना अपेक्षाकृत आसान है, जो इसे शैक्षिक उद्देश्यों या विशिष्ट उच्च-परिशुद्धता अनुप्रयोगों के लिए एक मजबूत उम्मीदवार बनाता है।

नोट: शोध पत्र की तिथि 14 मार्च, 2026 (Pi Day) है, जो यह संकेत देती है कि यह एक भविष्य-दिनांकित या काल्पनिक प्रकाशन है जिसे "Pi Day" उत्सव के लिए बनाया गया है, हालांकि इसकी गणितीय सामग्री एक कठोर व्युत्पत्ति के रूप में प्रस्तुत की गई है।