रुंग-कुता विधियाँ

thumb|रुंग -कुता विधियाअँ संख्यात्मक विश्लेषण में रुंगा-कुटा विधियाँ (Runge–Kutta methods) साधारण अवकल समीकरणों को हल करने की पुनरावृत्तिमूलक विधियाँ हैं। इनका विकास जर्मनी के गणितज्ञों सी रुंगा तथा एम डब्ल्यू कुटा ने १९०० ई के आसपास किया था। यह एक विधि नहीं बल्कि इसमें कई विधियाँ हैं। रुंगा-कुटा विधि इसे 'क्लासिकल रुंगा-कुटा विधि' या प्रायः "RK4" भी कहा जाता है। यह चार आर्डर वाली विधि है। मान लीजिए कि एक साधारण अवकल समीकरण (आरम्भिक मान समस्या) निम्नलिखित है: यह विधि निम्नलिखित सूत्र से समझी जा सकती है: n = 0, 1, 2, 3, . . . , तथा Press Teukolsky Vetterling Flannery 2007 908; Süli Mayers 2003 328 h का उचित धनात्मक मान लेकर इस विधि को बारबार प्रयोग करके t के किसी भी मान के लिए y का सन्निकट मान निकाल सकते हैं। उदाहरण माना का हल निकालना है। दिया हुआ है कि इस अवकल समीकरण का ठीक-ठीक हल (इग्जैक्ट सलुशन) है। हम लेकर इसका हल निकालते हैं। इससे पर प्राप्त होता है। इसी तरह हम t के विभिन्न मानों के लिये x का मान प्राप्त करते जाते हैं, जो निम्नांकित सारणी में दिखाये गये हैं- 0.0 1.0 0.1 0.994987426585 0.2 0.979795852198 0.3 0.95393908717 0.4 0.916514893222 0.5 0.866024896597 0.6 0.799998909634 0.7 0.714140165921 0.8 0.599991210485 0.9 0.435832710519 1.0 0.0488018582123 ध्यान दें कि इस अवकल समीकरण के विशुद्ध हल से पर मिलेगा, जबकि इस विधि से प्राप्त हुआ है जो विशुद्ध मान के काफी करीब है। इससे भी अधिक परिशुद्ध मान की गणना के लिये 'स्टेप साइज' को 0.1 के बजाय और कम रखना पड़ेगा। स्पष्‍ट रुंगा-कुटा विधियाँ (Explicit Runge–Kutta methods) ये विधियाँ RK4 के सामान्यीकृत रूप हैं। यह निम्नलिखित रूप में है: जहाँ Press Teukolsky Vetterling Flannery 2007 907 (टिप्पणी: अलग-अलग पुस्तकों में यही समीकरण अलग-अलग तरह से पारिभाषित किया हुआ मिलता है किन्तु वे इसके तुल्य ही होते हैं।) किसी विशेष विधि को प्राप्त करने के लिए पहले s (चरणों की संख्याँ) तय करनी पड़ती है। इसके अनुसार गुणांक aij (1 ≤ j < i ≤ s), bi (i = 1, 2, ..., s) और ci (i = 2, 3, ..., s) तय किए जाते हैं। मैट्रिक्स [aij] को रुंगा-कुटा मैट्रिक्स कहते हैं। bi तथा ci को भार (weights) एवं नोड (nodes) कहते हैं।Iserles 1996 38 ये संख्याएँ निम्नलिखित सारणी के रूप में व्यवस्थित की जातीं हैं जिसे 'बूचर टेबुल' कहते हैं। 0 रुंगा कुटा विधि एकरूप (consistent) होगी यदि यदि हम किसी विधि को p आर्डर वाला बनाना चाहते हैं तो कुछ और आवश्यकताएँ इसमें जुड़ जाती हैं। उदाहरण के लिए 2-चरण तथा 2-आर्डर की विधि प्राप्त करने के लिए निम्नलिखित शर्त आएगी: b1 + b2 = 1, b2c2 = 1/2, and a21 = c2.Iserles 1996 39 सन्दर्भ इन्हें भी देखें अवकल समीकरण साधारण अवकल समीकरण के हल की संख्यात्मक विधियाँ श्रेणी:अवकल समीकरण श्रेणी:संख्यात्मक विश्लेषण