Runge-Kutta 2. Mertebe Yöntemi Hesaplayıcı nedir?
Bu araç, \(y' = F(x, y)\) biçimindeki birinci mertebe adi diferansiyel denklemi (ADD) \([x_0, x_n]\) aralığında, \(y_0 = f(x_0)\) başlangıç koşulundan yola çıkarak sayısal olarak çözer. 2. mertebe Runge-Kutta yöntemini (orta nokta yöntemi) kullanır; \((x, y)\) yaklaşık değerlerinden oluşan bir tablo ile birlikte son değeri \(y_n = f(x_n)\) verir. Herhangi bir ülke veya yasal düzenlemeyle sınırlı olmayan, evrensel bir matematik aracıdır.
Nasıl kullanılır?
Sağ tarafı, yani \(F(x,y)\) ifadesini \(x\) ve \(y\) cinsinden bir matematiksel ifade olarak girin (örneğin 1-y^2, x*y veya sin(x)+y). Başlangıç noktası \(x_0\) ve \(y_0\) değerlerini, aralığın bitiş noktası \(x_n\)'yi ve eşit alt bölüm sayısı \(n\)'i belirleyin. Aralık, \(h = (x_n - x_0)/n\) boyutunda \(n\) adıma bölünür. Daha büyük bir \(n\), daha küçük bir adım ve daha yüksek doğruluk sağlar. Gösterim hassasiyeti seçeneği yalnızca kaç anlamlı basamağın ekranda gösterileceğini belirler, hesaplamayı etkilemez.
Formülün açıklaması
Orta nokta Runge-Kutta şeması çözümü adım adım ilerletir:
$$k_1 = h \cdot F(x_i, y_i)$$
$$k_2 = h \cdot F\left(x_i + \frac{h}{2},\; y_i + \frac{k_1}{2}\right)$$
$$y_{i+1} = y_i + k_2, \quad x_{i+1} = x_i + h$$
Eğim, adımın orta noktasında tahmin edilir; bu da baskın hata terimini ortadan kaldırır. Yerel kesme hatası \(O(h^3)\), genel hata ise \(O(h^2)\) mertebesindedir; dolayısıyla \(h\) değerini yarıya indirmek hatayı kabaca dörtte birine düşürür.
Çözümlü örnek
\(x_0 = 0\), \(y_0 = 0\), \(x_n = 1\), \(n = 50\) (yani \(h = 0.02\)) ile \(y' = 1 - y^2\) denklemini çözelim. Tam çözüm \(y = \tanh(x)\)'tir. 1. adım: $$k_1 = 0.02 \cdot (1-0) = 0.02$$ $$k_2 = 0.02 \cdot (1-0.01^2) = 0.019998$$ $$y_1 = 0.019998$$ 50 adımın tamamı sürdürüldüğünde \(y(1) \approx 0.76159\) elde edilir; bu da beş ondalık basamağa kadar \(\tanh(1) \approx 0.7615942\) ile uyuşur.
Sıkça sorulan sorular
Ne kadar doğru sonuç verir? Genel hata \(h^2\) ile orantılı olduğundan, \(n\) büyüdükçe doğruluk artar. Katı (stiff) denklemlerde ya da çok büyük adımlarda sonuç ıraksayabilir.
\(x_n\), \(x_0\)'dan küçük olabilir mi? Evet. Bu durumda \(h\) negatif olur ve yöntem \(x\) boyunca geriye doğru integral alır; bu da geçerlidir.
Hangi fonksiyonları kullanabilirim? Standart fonksiyonlar: sin, cos, tan, exp, ln, log, sqrt; ayrıca \(+\), \(-\), \(*\), \(/\), \(\hat{}\) işlemleri, parantezler ve \(e\) ile \(\pi\) sabitleri.
