Determining the Rolle function in Hermite interpolatory approximation by solving an appropriate differential equation

तकनीकी सारांश: हर्मिट इंटरपोलेटरी एप्रोक्सिमेशन में रोल फंक्शन (Rolle Function) का निर्धारण

समस्या विवरण
यह शोध पत्र हर्मिट इंटरपोलेटरी एप्रोक्सिमेशन (Hermite interpolatory approximations) की सटीकता में सुधार करने की चुनौती को संबोधित करता है। जबकि मानक हर्मिट बहुपद $H_{2n+1}(x)$ एक फलन $f(x)$ और उसके अवकलज (derivative) $f'(x)$ को $n+1$ नोड्स पर इंटरपोलेट करता है, पॉइंटवाइज एरर (pointwise error) $\Delta(x | H_{2n+1})$ एक पद द्वारा नियंत्रित होता है जिसमें अज्ञात बिंदु $\xi(x)$ पर $f$ का $(2n+2)$-वाँ अवकलज शामिल होता है, जिसे रोल फंक्शन (Rolle function) कहा जाता है। मुख्य समस्या यह है कि $\xi(x)$ सामान्यतः अज्ञात होता है, जिससे त्रुटि पद (error term) की सीधी गणना या सुधार करना असंभव हो जाता है। लेखक $f(x)$ के अधिक सटीक एप्रोक्सिमेशन की सुविधा के लिए $\xi(x)$ को संख्यात्मक रूप से निर्धारित करने की एक विधि चाहते हैं।

कार्यप्रणाली (Methodology)
प्रस्तावित कार्यप्रणाली लाग्रेंज इंटरपोलेशन (Lagrange interpolation) पर पहले से लागू की गई एक तकनीक को हर्मिट मामले के लिए विस्तारित करती है। इस दृष्टिकोण में निम्नलिखित चरण शामिल हैं:

1. एक अवकल समीकरण (Differential Equation) का व्युत्पन्न: हर्मिट इंटरपोलेशन के मानक एरर फॉर्मूला से शुरू करते हुए, लेखक त्रुटि पद का $x$ के सापेक्ष अवकलन करते हैं। चेन रूल (chain rule) लागू करके और रोल फंक्शन के अवकलज, $\frac{d\xi}{dx}$ को अलग करके, एक प्रथम-क्रम साधारण अवकल समीकरण (first-order ODE) व्युत्पन्न किया जाता है। यह ODE $\frac{d\xi}{dx}$ को ज्ञात फलनों $f(x)$, $H_{2n+1}(x)$ और उनके अवकलजों, साथ ही अज्ञात $\xi(x)$ से जोड़ता है।
2. प्रारंभिक मान समस्या (Initial Value Problem) का निरूपण: ODE को हल करने के लिए, एक नोड के निकट (परंतु उस पर नहीं) एक बिंदु $x_z$ पर एक प्रारंभिक मान $\xi_z = \xi(x_z)$ की आवश्यकता होती है। चूंकि नोड्स पर त्रुटि शून्य होती है, इसलिए नोड्स पर $\xi$ निर्धारित नहीं किया जा सकता है। लेखक एक नोड के पास एक बिंदु $x_z$ चुनते हैं और एक वैध प्रारंभिक $\xi_z$ खोजने के लिए बीजगणितीय त्रुटि समीकरण को संख्यात्मक रूप से हल करते हैं।
3. संख्यात्मक समाधान (Numerical Solution): व्युत्पन्न ODE को सातवें क्रम के रनगे-कुट्टा (Runge-Kutta - RK) मेथड का उपयोग करके संख्यात्मक रूप से हल किया जाता है। शोध पत्र नोट करता है कि प्रारंभिक मान के लिए कई गणितीय समाधान मौजूद हो सकते हैं; सही समाधान का चयन इस प्रतिबंध के आधार पर किया जाता है कि $\xi(x)$ को इंटरपोलेशन नोड्स द्वारा परिभाषित अंतराल ($x_0 < \xi(x) < x_n$) के भीतर ही होना चाहिए।
4. त्रुटि सुधार (Error Correction): एक बार $\xi(x)$ निर्धारित हो जाने के बाद, त्रुटि फॉर्मूला में $f^{(2n+2)}(\xi(x))$ पद ज्ञात हो जाता है। लेखक इस पद को एक बहुपद $H_\xi(x)$ (लीस्ट-स्क्वायर फिटिंग या क्यूबिक स्प्लाइन के माध्यम से) का उपयोग करके अनुमानित करते हैं। इस बहुपद को नोडल फैक्टर $\prod (x-x_k)^2$ से गुणा किया जाता है और त्रुटि बहुपद $E(x)$ बनाने के लिए $(2n+2)!$ से विभाजित किया जाता है।
5. परिष्कृत एप्रोक्सिमेशन (Refined Approximation): अंतिम, बेहतर एप्रोक्सिमेशन को मूल हर्मिट बहुपद में त्रुटि बहुपद जोड़कर बनाया जाता है: $f(x) \approx H_{2n+1}(x) + E(x)$।

प्रमुख परिणाम
इस कार्यप्रणाली का परीक्षण फलन $f(x) = e^x \sin x$ पर नोड्स $\{0, 3\pi/2\}$ (जहाँ $n=1$) के लिए किया गया।

• रोल फंक्शन का निर्धारण: संख्यात्मक समाधान ने सफलतापूर्वक एक रोल फंक्शन $\xi(x)$ उत्पन्न किया जो अंतराल संबंधी प्रतिबंधों को संतुष्ट करता था। ODE के संख्यात्मक समाधान को त्रुटि समीकरण के बाएँ-पक्ष और दाएँ-पक्ष की तुलना करके सत्यापित किया गया, जिससे पता चला कि वे अनिवार्य रूप से अविभेद्य (indistinguishable) हैं।
• सटीकता में सुधार:
  • मूल हर्मिट बहुपद $H_3(x)$ ने अधिकतम एप्रोक्सिमेशन एरर 7.04 प्रदर्शित किया।
  • डिग्री 11 के लीस्ट-स्क्वायर फिट के माध्यम से प्राप्त एरर बहुपद $E(x)$ को जोड़ने से, अधिकतम एरर घटकर $7.3 \times 10^{-8}$ रह गया।
  • अवकलज पद के एप्रोक्सिमेशन के लिए क्यूबिक स्प्लाइन का उपयोग करने से अधिकतम एरर और घटकर लगभग $10^{-12}$ हो गया।
• इंटरपोलेटरी गुणों का संरक्षण: शोध पत्र सत्यापित करता है कि सुधारा गया एप्रोक्सिमेशन $H_{2n+1}(x) + E(x)$ मूल हर्मिट इंटरपोलेटरी गुणों को बनाए रखता है, जो नोड्स पर फलन के मान और अवकलज दोनों से मेल खाता है।
• संख्यात्मक समाकलन (Numerical Integration): शोध पत्र संक्षेप में प्रदर्शित करता है कि सुधारा गया बहुपद अंतराल पर समाकलन करने से मूल हर्मिट बहुपद ($\sim 0.7$ एरर) की तुलना में काफी अधिक सटीक परिणाम ($\sim 3 \times 10^{-12}$ एरर) प्राप्त होता है।

महत्व और दावे
शोध पत्र का दावा है कि प्रारंभिक-मान समस्या (initial-value problem) के माध्यम से रोल फंक्शन को निर्धारित करने से हर्मिट इंटरपोलेशन में स्पष्ट पॉइंटवाइज एरर का अनुमान लगाना संभव है। इसकी प्राथमिक महत्ता इस क्षमता में निहित है कि एक नया बहुपद एप्रोक्सिमेशन बनाया जा सकता है जो मानक हर्मिट बहुपद की तुलना में कई गुना अधिक सटीक है, जबकि मूल इंटरपोलेशन स्थितियों को सख्ती से बनाए रखता है।

लेखक नोट करते हैं कि अंतिम परिणाम की सटीकता वर्तमान में $\xi(x)$ के लिए उपयोग किए गए RK मेथड की शुद्धता द्वारा सीमित है, जो यह सुझाव देता है कि भविष्य के कार्य RK सॉल्वर के भीतर एरर कंट्रोल पर ध्यान केंद्रित कर सकते हैं। शोध पत्र निष्कर्ष निकालता है कि यह तकनीक, लाग्रेंज इंटरपोलेशन के पिछले कार्यों पर आधारित, संख्यात्मक एप्रोक्सिमेशन की गुणवत्ता में सुधार करने की पर्याप्त क्षमता प्रदर्शित करती है।