A root finding method with arbitrary order of convergence

तकनीकी सारांश: मनचाहे अभिसरण क्रम (Order of Convergence) के साथ एक मूल खोज विधि (Root Finding Method)

लेखक: अलोइस शीस (Alois Schiessl)
दिनांक: 10 फरवरी, 2026

1. समस्या विवरण

यह शोध पत्र एक धनात्मक वास्तविक संख्या $a$ के $M$-वें मूल ($\sqrt[M]{a}$) की गणना करने की मौलिक संख्यात्मक समस्या को संबोधित करता है। जबकि न्यूटन की विधि (Newton’s Method) जैसे पारंपरिक तरीकों का व्यापक रूप से उपयोग किया जाता है, वे एक निश्चित द्विघात अभिसरण क्रम (quadratic order of convergence - क्रम 2) तक सीमित हैं, जिसका अर्थ है कि प्रत्येक पुनरावृत्ति (iteration) के साथ सही अंकों की संख्या लगभग दोगुनी हो जाती है। लेखक एक ऐसी विधि चाहते हैं जो केवल बहुपद मूल्यांकन (polynomial evaluations) तक सीमित हो—जो कि गणनात्मक रूप से कुशल हैं—लेकिन एक पैरामीटर $P$ को समायोजित करके मनचाहे अभिसरण क्रम $(P+1)$ की अनुमति दे सके।

2. कार्यप्रणाली

प्रस्तावित विधि मानक फलन $f(x) = x^M - a$ के बजाय एक विशिष्ट फिक्स्ड-पॉइंट इटरेशन (fixed-point iteration) पर ध्यान केंद्रित करती है, जो फलन $f(t) = 1 - \frac{t^M}{a}$ पर आधारित है।

फिक्स्ड-पॉइंट फलन:
लेखक एक निश्चित समाकलन (definite integral) का उपयोग करके एक विशेष फिक्स्ड-पॉइंट फलन $F(x)$ को परिभाषित करते हैं:
$$F(x) = \left[ \prod_{\ell=1}^{P} \left( 1 + \frac{1}{\ell \cdot M} \right) \right] x \int_{0}^{x} \left( 1 - \frac{t^M}{a} \right)^P dt$$

प्रमुख गणितीय चरण:

• बहुपद न्यूनीकरण (Polynomial Reduction): इंटीग्रेंड (integrand) पर द्विपद प्रमेय (Binomial Theorem) लागू करके और पद-दर-पद समाकलन करके, लेखक सिद्ध करते हैं कि $F(x)$, डिग्री $(P \cdot M + 1)$ वाला एक बहुपद है जिसमें केवल $(P + 1)$ पद हैं।
• गुणांक पूर्व-गणना (Coefficient Pre-calculation): बहुपद के गुणांक $c_k$ इस प्रकार व्युत्पन्न किए गए हैं कि वे इटरेशन प्रक्रिया के दौरान स्थिर रहें, जिससे हॉर्नर स्कीम (Horner’s scheme) के माध्यम से कुशल गणना संभव हो सके।
• अभिसरण विश्लेषण (Convergence Analysis): बनाच फिक्स्ड-पॉइंट थ्योरम (Banach Fixed-Point Theorem) और टेलर श्रेणी विस्तार का उपयोग करते हुए, लेखक सिद्ध करते हैं कि:
  1. $\sqrt[M]{a}$, $F(x)$ का एक फिक्स्ड पॉइंट है।
  2. फिक्स्ड पॉइंट पर $F(x)$ के पहले $P$ अवकलज (derivatives) शून्य हैं।
  3. $(P+1)$-वाँ अवकलज गैर-शून्य है।
  4. फलस्वरूप, अनुक्रम $x_{n+1} = F(x_n)$, $\sqrt[M]{a}$ की ओर ठीक $P+1$ के अभिसरण क्रम के साथ अभिसरित होता है।

3. प्रमुख योगदान

• मनचाहा अभिसरण क्रम (Arbitrary Convergence Order): न्यूटन की विधि के विपरीत, उपयोगकर्ता $P$ को समायोजित करके किसी भी वांछित क्रम (जैसे कि क्यूबिक, क्वार्टिक, क्विंटिक आदि) को प्राप्त करने के लिए चुन सकता है।
• केवल बहुपद दृष्टिकोण (Polynomial-Only Approach): यह विधि इटरेशन चरण के दौरान जटिल ट्रांसेंडेंटल फलनों (transcendental functions) से बचती है, जिससे समस्या सरल बहुपद मूल्यांकन तक सीमित हो जाती है।
• विरलता (Sparsity) के माध्यम से दक्षता: उच्च डिग्री के बावजूद, गैर-शून्य पदों की संख्या को $(P+1)$ तक कम रखा गया है, जो इसे उच्च-परिशुद्धता अंकगणित (high-precision arithmetic) के लिए गणनात्मक रूप से हल्का बनाता है।

4. परिणाम और गणनात्मक सत्यापन

लेखक दो प्राथमिक गणनात्मक प्रयोगों के माध्यम से विधि को सत्यापित करते हैं:

• उदाहरण 1 (द्विघात अभिसरण): $P=1$ (क्रम 2) का उपयोग करके $\sqrt[3]{10}$ की गणना करना। परिणामों ने प्रदर्शित किया कि जैसा कि अपेक्षित था, प्रत्येक चरण के साथ सही अंकों की संख्या लगभग दोगुनी हो गई।
• उदाहरण 2 (उच्च-परिशुद्धता क्वार्टिक अभिसरण): $P=3$ (क्रम 4) का उपयोग करके $\sqrt{2}$ की दस लाख दशमलव अंकों तक गणना करना।
  • प्रदर्शन: गणना MAPLE सॉफ़्टवेयर का उपयोग करके एक मानक उपभोक्ता-ग्रेड पीसी (AMD Ryzen 5 5600G, 32GB RAM) पर लगभग 3 सेकंड में पूरी हुई।
  • सत्यापन: त्रुटि का परिमाण (पुनरावृत्तियों के बीच के अंतर में अग्रणी शून्य की संख्या) प्रत्येक चरण में लगभग चार गुना बढ़ गया, जो क्वार्टिक अभिसरण क्रम की पुष्टि करता है।

5. महत्व

यह शोध पत्र उच्च-परिशुद्धता मूल खोज (root finding) के लिए एक अत्यधिक स्केलेबल एल्गोरिदम प्रस्तुत करता है। इसका महत्व इसकी पूर्वानुमेयता और लचीलेपन में निहित है: कम्प्यूटेशनल गणित या क्रिप्टोग्राफी जैसे अनुप्रयोगों के लिए जहाँ अत्यधिक परिशुद्धता की आवश्यकता होती है, यह विधि एक डेवलपर को बहुपद की जटिलता (जो $P$ द्वारा निर्धारित होती है) और अभिसरण की गति के बीच संतुलन बनाने की अनुमति देती है। दस लाख अंकों की परिशुद्धता तक कुछ ही सेकंड में पहुँचने की क्षमता इसे मानक इटरेटिव सॉल्वर का एक सशक्त विकल्प बनाती है।