ऑयलर विधि कैलकुलेटर
dy/dx = f(x, y), एक प्रारंभिक बिंदु और अंत-मान दर्ज करें, और ऑयलर विधि, सुधरी हुई ऑयलर (हून) विधि या चतुर्थ-कोटि रुंगे–कुट्टा से अनुमानित हल पाएँ। इसमें हर प्रवणता सहित पूरी चरण तालिका, सटीक हल (या एक सटीक संदर्भ मान) से त्रुटि, चरण आधा करने पर क्या होता है, और ऊपर खिंचे हल के साथ प्रवणता-क्षेत्र का ग्राफ़ शामिल है। निःशुल्क, पंजीकरण आवश्यक नहीं।
x को स्वतंत्र चर और y को अज्ञात के लिए उपयोग करें, और + − * / ^ तथा कोष्ठकों के साथ: x + y, x*y, y*(1 − y), sin(x) − y। फ़ंक्शन: sin, cos, tan, exp, ln, sqrt और अन्य।
यदि आप सटीक हल जानते हैं, तो उसे x में दर्ज करें ताकि वास्तविक त्रुटि दिखे; अन्यथा एक बहुत सटीक RK4 संदर्भ उपयोग होता है।
परिणाम
3.1874849
y(1) का अनुमान
3.4365637
सटीक y(1)
-0.249
त्रुटि
0.1
चरण आकार h
10
चरण n
विधि
yₙ₊₁ = yₙ + h·f(xₙ, yₙ): हर चरण की शुरुआत की प्रवणता का अनुसरण करें। वैश्विक त्रुटि h के समानुपाती है।
चरण तालिका
| n | xₙ | yₙ | k1 | सटीक y(xₙ) | त्रुटि |
|---|---|---|---|---|---|
| 0 | 0 | 1 | 1 | 1 | 0 |
| 1 | 0.1 | 1.1 | 1.2 | 1.1103418 | -0.0103 |
| 2 | 0.2 | 1.22 | 1.42 | 1.2428055 | -0.0228 |
| 3 | 0.3 | 1.362 | 1.662 | 1.3997176 | -0.0377 |
| 4 | 0.4 | 1.5282 | 1.9282 | 1.5836494 | -0.0554 |
| 5 | 0.5 | 1.72102 | 2.22102 | 1.7974425 | -0.0764 |
| 6 | 0.6 | 1.943122 | 2.543122 | 2.0442376 | -0.101 |
| 7 | 0.7 | 2.1974342 | 2.8974342 | 2.3275054 | -0.13 |
| 8 | 0.8 | 2.4871776 | 3.2871776 | 2.6510819 | -0.164 |
| 9 | 0.9 | 2.8158954 | 3.7158954 | 3.0192062 | -0.203 |
| 10 | 1 | 3.1874849 | 3.4365637 | -0.249 |
एक ही चरण के साथ तीन विधियाँ
| विधि | अंत पर y | त्रुटि |
|---|---|---|
| ऑयलर विधि | 3.1874849 | -0.249 |
| सुधरी हुई ऑयलर (हून) | 3.4281617 | -0.0084 |
| रुंगे–कुट्टा 4 | 3.4365595 | -0.00000417 |
चरण आधा करना (ऑयलर विधि)
| n | h | अंत पर y | त्रुटि |
|---|---|---|---|
| 10 | 0.1 | 3.1874849 | -0.249 |
| 20 | 0.05 | 3.3065954 | -0.13 |
| 40 | 0.025 | 3.3701277 | -0.0664 |
| 80 | 0.0125 | 3.4029699 | -0.0336 |
ऑयलर विधि के लिए h आधा करने पर त्रुटि लगभग आधी होती है; हून के लिए एक-चौथाई; RK4 के लिए एक-सोलहवाँ।
ग्राफ़
- सटीक हल
- ऑयलर विधि
ऑयलर विधि कैलकुलेटर का उपयोग कैसे करें
- 1. समीकरण दर्ज करें। x और y का उपयोग करके dy/dx = f(x, y) का दायाँ पक्ष लिखें:
x + y,y*(1 - y),sin(x) - y। - 2. प्रारंभिक शर्त और अंत दें। x₀, y₀ = y(x₀) और वह x-मान दर्ज करें जहाँ तक हल करना है। अंत x₀ के बाईं ओर भी हो सकता है।
- 3. चरण चुनें। चरण आकार h दर्ज करें, या चरणों की संख्या n पर स्विच करें। यदि h अंतराल को समान रूप से नहीं बाँटता, तो कैलकुलेटर उसे समायोजित करके बताता है।
- 4. विधि चुनें और तालिका पढ़ें। हर पंक्ति में xₙ, yₙ और उपयोग की गई प्रवणताएँ दिखती हैं। यदि आप सटीक हल जानते हैं, तो उसे दर्ज करें ताकि वास्तविक त्रुटि दिखे।
तीन विधियाँ
| विधि | चरण | वैश्विक त्रुटि |
|---|---|---|
| ऑयलर | yₙ₊₁ = yₙ + h·f(xₙ, yₙ) | h के समानुपाती |
| सुधरी हुई ऑयलर (हून) | yₙ₊₁ = yₙ + (h/2)(k₁ + k₂), जहाँ k₁ = f(xₙ, yₙ) और k₂ = f(xₙ + h, yₙ + h·k₁) | h² के समानुपाती |
| रुंगे–कुट्टा 4 | yₙ₊₁ = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄) | h⁴ के समानुपाती |
ऑयलर विधि हर चरण की शुरुआत में स्पर्शरेखा रेखा का अनुसरण करती है। हून विधि चरण की शुरुआत की प्रवणता को अनुमानित अंत की प्रवणता के साथ औसत करती है, इसीलिए इसे स्पष्ट ट्रैपेज़ॉइड नियम भी कहा जाता है। रुंगे–कुट्टा 4 प्रति चरण प्रवणता चार बार नमूना लेती है और उन्हें 1, 2, 2, 1 भार देती है। ये Burden and Faires की Numerical Analysis के अध्याय 5 की विधियाँ हैं।
हल किया गया उदाहरण: y′ = x + y, y(0) = 1
सटीक हल y = 2eˣ − x − 1 है, इसलिए y(1) = 2e − 2 ≈ 3.4365637। h = 0.1 के साथ ऑयलर विधि के पहले चरण ये हैं:
| n | xₙ | yₙ | प्रवणता f(xₙ, yₙ) | सटीक y(xₙ) |
|---|---|---|---|---|
| 0 | 0 | 1 | 1 | 1 |
| 1 | 0.1 | 1.1 | 1.2 | 1.1103418 |
| 2 | 0.2 | 1.22 | 1.42 | 1.2428055 |
| 3 | 0.3 | 1.362 | 1.662 | 1.3997176 |
उदाहरण के लिए, y₁ = 1 + 0.1 × (0 + 1) = 1.1 और y₂ = 1.1 + 0.1 × (0.1 + 1.1) = 1.22। दस चरणों के बाद ऑयलर विधि y(1) ≈ 3.1874849 तक पहुँचती है, यानी −0.2491 की त्रुटि; यह कम आँकती है, क्योंकि सटीक हल उत्तल है और स्पर्शरेखाएँ उसके नीचे रहती हैं।
| चरण n (h) | x = 1 पर ऑयलर त्रुटि | हून त्रुटि | RK4 त्रुटि |
|---|---|---|---|
| 10 (0.1) | −0.2491 | −0.0084 | −0.0000042 |
| 20 (0.05) | −0.1300 | −0.0022 | −0.00000027 |
| 40 (0.025) | −0.0664 | −0.00056 | −0.000000017 |
| 80 (0.0125) | −0.0336 | −0.00014 | −0.0000000011 |
h के आधा होने पर ऑयलर की त्रुटि लगभग आधी, हून की एक-चौथाई और RK4 की एक-सोलहवाँ रह जाती है, जैसा कोटि-तालिका बताती है। n = 10 पर RK4 की त्रुटि के बराबर पहुँचने के लिए ऑयलर विधि को कई लाख चरण चाहिए।
प्रवणता-क्षेत्र और त्रुटि कैसे पढ़ें
ग्राफ़ पर छोटे खंड बिंदुओं के एक ग्रिड पर प्रवणता f(x, y) दिखाते हैं: वास्तविक हल एक ऐसा वक्र है जो हर जगह इन खंडों को स्पर्श करता है। संख्यात्मक हल उसके ऊपर बहुभुज रेखा के रूप में खींचा जाता है; यदि आप सटीक हल देते हैं, तो वह उसके साथ बिंदुकित रेखा में दिखाया जाता है। सटीक हल न होने पर त्रुटि-स्तंभ RK4 के कम से कम 2,000 चरणों वाले संदर्भ मान से तुलना करता है, जो यहाँ चुनी जाने वाली किसी भी विधि से कहीं अधिक सटीक है। कुछ समीकरण परिमित समय में अनंत हो जाते हैं, जैसे x = 1 पर y(0) = 1 के साथ y′ = y²; जहाँ मान अनंत हो जाते हैं, तालिका वहीं रुक जाती है।
अक्सर पूछे जाने वाले प्रश्न
- ऑयलर विधि कैसे काम करती है?
- यह हल को चरण-दर-चरण बनाती है। बिंदु (xₙ, yₙ) से यह अवकल समीकरण के आधार पर प्रवणता f(xₙ, yₙ) निकालती है और उस स्पर्शरेखा रेखा पर क्षैतिज दूरी h तक चलती है, यानी yₙ₊₁ = yₙ + h·f(xₙ, yₙ)। हर नया बिंदु अगले चरण का आरंभिक बिंदु बनता है।
- ऑयलर विधि कितनी सटीक है?
- प्रति चरण इसकी स्थानीय त्रुटि h² के समानुपाती होती है, लेकिन त्रुटियाँ लगभग 1/h चरणों में जमा होती हैं, इसलिए वैश्विक त्रुटि h के समानुपाती है। चरण आधा करने पर त्रुटि भी आधी हो जाती है। यह सीखने के लिए अच्छी विधि है, परिशुद्धता के लिए नहीं।
- सुधरी हुई ऑयलर विधि क्या है?
- सुधरी हुई ऑयलर विधि, या हून की विधि, पहले अंतराल के अंत का अनुमान लगाने के लिए एक ऑयलर चरण लेती है, फिर वहाँ की प्रवणता को आरंभ की प्रवणता के साथ औसत करती है। इस द्वितीय-कोटि सुधार से वैश्विक त्रुटि h के बजाय h² के समानुपाती हो जाती है, जिसकी कीमत प्रति चरण दो प्रवणता-मूल्यांकन है।
- चरण आकार कैसे चुनूँ?
- आधा करने वाली तालिका का उपयोग करें: h तब तक घटाएँ जब तक उत्तर आपके लिए महत्वपूर्ण अंकों में बदलना बंद न कर दे। ऑयलर विधि को छोटे h की ज़रूरत होती है, लेकिन RK4 बहुत बड़े h पर भी उतनी ही सटीकता देती है। कठोर (स्टिफ़) समीकरण, जिनमें कुछ घटक तेज़ी से क्षय होते हैं, इन सभी विधियों पर छोटे h की माँग कर सकते हैं।
- क्या यह द्वितीय-कोटि या समीकरण-निकाय हल कर सकता है?
- नहीं। यह केवल एक प्रथम-कोटि समीकरण dy/dx = f(x, y) हल करता है। द्वितीय-कोटि समीकरण को प्रथम-कोटि निकाय में बदला जा सकता है, लेकिन यह उपकरण निकायों को नहीं संभालता।
- हल बीच में क्यों रुक जाता है?
- मान अनंत हो गए, जो तब होता है जब सटीक हल परिमित समय में अनंत हो जाता है, या किसी अस्थिर समीकरण के लिए चरण बहुत बड़ा हो। तालिका अंतिम परिमित पंक्ति पर समाप्त होती है, और संदेश वह x-मान बताता है जिसके निकट यह हुआ।