Mavzu: Birinchi tartibli differsial tenglamalarni yechish usullari

K₁⁽ⁱ⁾=hf_i(x_i,y_i)

K₂⁽ⁱ⁾=hf_i(x_i +h/2, y_i+K₁⁽ⁱ⁾/2)

K₃⁽ⁱ⁾=hf_i(x_i +h/2, y_i+K₂⁽ⁱ⁾/2) (7.5.1)

K₄⁽ⁱ⁾=hf_i(x_i +h, y_i+K₃⁽ⁱ⁾)
Funktsiyaning orttirmasi y_i ni quyidagi formuladan topiladi
y_i=(K₁⁽ⁱ⁾+2 K₂⁽ⁱ⁾+2 K₃⁽ⁱ⁾+ K₄⁽ⁱ⁾) / 6 (7.5.2)

Bu erda h=(b-a)/n – integrallash qadami. i ni har bir qiymati uchun (7.5.1) va (7.5.2) dagi amallarni bajaramiz va noma’lum funktsiya “u” ni qiymatlarini (tenglamaning yechimini) quyidagi formuladan topamiz.

y_i+1=y_i+ y_i , (i=0,1,2, ...n) 

Runge – Kutta usuli bilan differensial tenglamani yechishda jadval tuzilsa hisoblash jarayoni birmuncha osonlashadi. Jadvalni tuzish tartibi quyidagicha:

1. (2) va (3) ustunlarga x va u ning kerakli bo’lgan qiymatlari yoziladi.

   
   

2. “x” va “u” larning qiymatlarini ((2)-va (3)-ustunlardan) u^’=f(x,y) tenglamani o’ng tarafiga qo’yiladi va natijalarni (4) ustunga (satrlari mos ravishda) qo’yiladi.

   
   

3. Topilgan f(x,y) qiymatlarini integrallash qadami “h” ga ko’paytiriladi va natijalar (5) ustunga yoziladi.

   
   
4. 
   
   K₁⁽⁰⁾ ni 1 ga, K₂⁽⁰⁾ va K₃⁽⁰⁾ larni 2 ga, K₄⁽⁰⁾ ni 1 ga ko’paytirib ularni (6) ustunga yozamiz.
   

I-IV jarayonni Kⁱ ni (i=0,1,2, ...n) har bir qiymati uchun takrorlaymiz. (6)-ustunni qiymatlarining yig’indisini hisoblab, natijani 6 ga bo’lamiz va u=(1/6) (K₁⁽ⁱ⁾+2 K₂⁽ⁱ⁾+2 K₃⁽ⁱ⁾+ K₄⁽ⁱ⁾) ni

topamiz. Va nihoyat y_i+1=y_i+ y_i topiladi. YUqorida keltirilgan hisoblash tartibini [a,b] kesmani barcha nuqtalari uchun takrorlaymiz.