Bu hesaplayıcı ne işe yarar?
Bu araç, \(y' = F(x, y)\) biçimindeki birinci mertebeden adi diferansiyel denklemi (ADD), verilen bir \(y(x_0) = y_0\) başlangıç koşuluyla, \(x_0\) ile \(x_n\) arasındaki aralıkta sayısal olarak çözer. Bunun için sayısal analizde en yaygın kullanılan ve en güvenilir tek adımlı integratörlerden biri olan klasik dördüncü mertebe Runge-Kutta yöntemini (RK4) kullanır. Sonuç olarak gerçek çözüme yaklaşan \((x_i, y_i)\) noktalarından oluşan bir tablo ve son değer olan \(y(x_n)\) elde edilir. Bu, herhangi bir ülkeye veya birime bağlı olmayan, tamamen matematiksel bir araçtır.
Nasıl kullanılır?
Denklemin sağ tarafı olan \(F(x,y)\) ifadesini x ve y cinsinden girin (örneğin 1-y^2, x+y, x*y veya sin(x)+y). Desteklenen işlemler + - * / ^ operatörleri ve sin, cos, tan, exp, log, ln, sqrt, abs, tanh gibi fonksiyonlar ile pi ve e sabitleridir. Başlangıç noktası \(x_0\)'ı, başlangıç değeri \(y_0\)'ı, bitiş noktası \(x_n\)'i belirleyin ve aralığı kaç eşit parçaya (\(n\)) böleceğinizi seçin. Parça sayısı arttıkça doğruluk da artar; çünkü RK4'ün genel hatası \(h^4\) ile orantılı olarak küçülür.
Formülün açıklaması
Aralık, genişliği \(h = (x_n - x_0)/n\) olan \(n\) eşit adıma bölünür. RK4 her adımda eğimi dört kez örnekler: bir kez başlangıçta (\(k_1\)), iki kez orta noktada (\(k_2, k_3\)) ve bir kez de bitişte (\(k_4\)). Bir sonraki değer ise bunların ağırlıklı ortalamasıdır:
$$y_{n+1} = y_n + \frac{k_1 + 2k_2 + 2k_3 + k_4}{6}$$Bu yaklaşım dördüncü mertebeye kadar olan hata terimlerini ortadan kaldırarak \(O(h^5)\) yerel kesme hatası ve \(O(h^4)\) genel hata sağlar.
Örnek çözüm
\(y' = 1 - y^2\) denklemini \(x_0 = 0\), \(y_0 = 0\), \(x_n = 1\) ve \(n = 10\) (\(h = 0.1\)) için çözelim. Bu denklemin kesin çözümü \(y = \tanh(x)\)'tir. İlk RK4 adımı \(y_1 = 0.0996679\) verir; bu da \(\tanh(0.1) = 0.0996680\) ile örtüşür. On adımın tamamı bittiğinde \(y(1) = 0.7615942\) bulunur ve bu sonuç \(\tanh(1) = 0.7615942\) ile yedi haneye kadar uyumludur.
Sıkça sorulan sorular
RK4 neden Euler yönteminden daha iyidir? Euler yöntemi her adımda tek bir eğim kullanır (hata \(O(h)\)). RK4 ise dört eğim kullanıp bunların ortalamasını alır ve aynı adım boyutuyla \(O(h^4)\) doğruluk sağlar; böylece hedeflenen hassasiyete çok daha az adımda ulaşılır.
Kaç adım seçmeliyim? 50 ile başlayın. Çözüm pürüzsüzse bu genellikle yeterlidir; hızlı değişen veya katılığa yakın (near-stiff) problemlerde adım sayısını 100, 200 veya 500'e çıkarın.
NaN veya Sonsuz (Infinity) sonucu alırsam ne yapmalıyım? Çözüm ıraksamış ya da \(F(x,y)\) geçersiz bir işleme (negatif sayının logaritması veya sıfıra bölme gibi) takılmış olabilir. İfadeyi kontrol edin ve daha küçük bir aralık veya daha fazla adım deneyin.
