Bu hesap makinesi ne işe yarar?
Bu araç, \(y'' = F(x, y, y')\) biçiminde yazılan ikinci mertebeden bir adi diferansiyel denklemi (ODE), \([x_0, x_n]\) aralığı üzerinde klasik dördüncü mertebe Runge-Kutta (RK4) yöntemiyle sayısal olarak çözer. Sağ tarafı oluşturan \(F\) ifadesini \(x\), \(y\) ve \(p\) cinsinden bir matematiksel ifade olarak (burada \(p\), \(y'\) yerine geçer) girersiniz; ardından \(y(x_0)\) ve \(y'(x_0)\) başlangıç koşullarını, aralığın uç noktalarını ve adım sayısını belirtirsiniz. Hesap makinesi, \(x\), \(y\) ve \(y'\) değerlerinden oluşan eksiksiz bir tablo ile uç nokta değerlerini döndürür. Tamamen sayısal analizdir ve her yerde geçerlidir.
Nasıl kullanılır?
\(F(x, y, p)\) ifadesini girin — örneğin \(y'' + 4y' + 4y = 0\) için -4*p - 4*y. \(x_0\), \(y_0 = f(x_0)\) ve \(p_0 = y'(x_0)\) değerlerini ayarlayın. Bitiş noktası \(x_n\)'i belirleyin ve alt aralık sayısı \(n\)'i seçin (daha fazla adım, daha küçük bir \(h = (x_n - x_0)/n\) adım boyu ve daha yüksek doğruluk demektir). Kaç anlamlı basamak gösterileceğini seçin. RK4'ün hatası genel olarak \(O(h^4)\) mertebesindedir; bu yüzden \(n\)'i iki katına çıkarmak hatayı yaklaşık 16 kat azaltır.
Formülün açıklaması
İkinci mertebeden bir ODE, \(p = y'\) tanımı yapılarak iki adet birinci mertebeden denklemden oluşan bir sisteme indirgenir: bu durumda \(y' = p\) ve \(p' = F(x, y, p)\) olur. RK4, her adımda her iki bilinmeyeni de dört ağırlıklı eğim hesabıyla ilerletir (\(y\) için \(k\), \(p\) için \(j\)) ve bunları \((k_1 + 2k_2 + 2k_3 + k_4)/6\) ile \((j_1 + 2j_2 + 2j_3 + j_4)/6\) şeklinde birleştirir.
Çözümlü örnek
\(y'' = -4p - 4y\) denklemini ele alalım; \(x_0 = 0\), \(y_0 = 0\), \(p_0 = 1\), \(x_n = 1\) ve tek bir kaba adım (\(n = 1\), \(h = 1\)). Dört aşama şu sonuçları verir: \(k_1=1\), \(k_2=-1\), \(k_3=2\), \(k_4=-5\) ve \(j_1=-4\), \(j_2=2\), \(j_3=-6\), \(j_4=12\). Buradan $$y(1) = \frac{1 - 2 + 4 - 5}{6} = -\frac{1}{3} = -0.3333$$ ve $$p(1) = 1 + \frac{-4 + 4 - 12 + 12}{6} = 1$$ bulunur. Kesin çözüm \(y = x e^{-2x}\) olduğundan \(y(1) = e^{-2} = 0.1353\)'tür; \(h = 1\)'lik tek bir dev adım çok kabadır. \(n = 50\) ya da \(100\) kullanmak sonucu tam olarak kesin değere yakınsatır.
Sıkça sorulan sorular
Geriye doğru integral alabilir miyim? Evet — \(x_n\) değerini \(x_0\)'dan küçük yapın; bu durumda \(h\) adım boyu negatif olur ve \(x_0\)'dan \(x_n\)'e doğru aşağı yönde integral alınır.
Az adımda sonucum neden yanlış görünüyor? RK4'ün doğruluğu küçük bir adım boyuna bağlıdır. Ardışık sonuçlar değişmeyi bırakana kadar \(n\)'i artırın.
F içinde hangi fonksiyonları kullanabilirim? + - * / ^, parantezler ve sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, exp, log/ln, log10, sqrt, abs ile birlikte \(\pi\) ve \(e\) sabitleri.
