पीexp( z))पी
पीपी( z)) = ∑ज = १पीzजेजे !
| पी5( मैं ε ) | > 1εपी5( z))
पी5( z))
पी5( z))पीपी5( z))
अंत में, उच्च-क्रम रन-कुट्टा विधियों के लिए काल्पनिक स्थिरता अंतराल की सीमा निर्धारित करते समय गलतियां करना आसान है। ऐसा इसलिए है क्योंकि ऐसे तरीकों के लिए स्थिरता क्षेत्र की सीमा काल्पनिक अक्ष के बेहद करीब है । इसलिए, राउंडऑफ़ त्रुटियां गलत निष्कर्ष पर ले जा सकती हैं; केवल सटीक गणना का उपयोग किया जाना चाहिए (निश्चित रूप से, इन परिस्थितियों में व्यावहारिक उद्देश्यों के लिए स्थिरता क्षेत्र की सीमा की प्रासंगिकता पर निश्चित रूप से बहस की जा सकती है)।
उदाहरण के लिए, यहां पांचवीं क्रम विधि के स्थिरता क्षेत्र का एक भाग है, जो फेहेलबर्ग 5 (4) युग्म से लिया गया है:
काल्पनिक स्थिरता अंतराल खाली है, लेकिन आप इस संकल्प पर चित्र से नहीं बता सकते हैं! ध्यान दें कि इस क्षेत्र में स्पष्ट रूप से काल्पनिक अक्ष का हिस्सा शामिल है, लेकिन उत्पत्ति के बारे में कोई अंतराल नहीं है।
इस बीच, डोरमंड-प्रिंस 5 (4) जोड़ी से पांचवें क्रम की विधि के लिए यहां है:
[ - १ , १ ]
पीपी( z))
आपको NodePy पैकेज में भी रुचि हो सकती है , जो ऊपर दिए गए भूखंडों का उत्पादन करता है और जिनका उपयोग किसी विधि की काल्पनिक स्थिरता अंतराल (अस्वीकरण: मैंने NodePy बनाया) जैसी चीजों को सटीक रूप से निर्धारित करने के लिए किया जा सकता है।