Metode prediktor–korektor

Metode Prediktor-korektor

metode Trapesium-Euler

Metode-metode yang sudah dibahas pada bagian-bagian sebelumnya yaitu Metode Euler dan Runge-kutta merupakan metode satu langkah untuk menyeleseikan persamaan diferensial biasa.sekarang kita akan membahas metode multi langkah,untuk menghitung y_k dengan menggunakan gradien-gradien f_j,dengan j < k,yang sudah diperoleh sebelumnya.metode ini tidak dapat dimulai dengan sendirinya karena tergantung pada metode-metode satu langkah seperti metode Euler untuk mendapatkan beberapa gradien awal.

metode prediktor-korektor terdiri atas dua bagian:(1) bagian prediktor,yang memprediksi y_k dengan menggunakan gradien-gradien f_j (j < k),dan (2)bagian korektor,yang menggunakan suatu rumus integrasi untuk memperbaiki hampiran.

Metode Trapesium-Euler menggunakan metode Euler sebagai algoritma korektor.jika kita gunakan indek pertama untuk menunjukan interval(langkah)dab indek kedua untuk menunjukan urutan hampiran,maka rumus Euler dapat ditulis sebagai

         yk+1,0=yk,* +hfk,*                                              

dengan aturan rumus dan^'*' berturut-turut menunjukan hampiran awal dan akhir.pada rumus Euler,y_k,* = y_k = y(t_k),dan f_k,* = f(t_k,y_k).

Sebagai persamaan korektor,aturan trapesium dinyatakan sebagai

      yk+1,j = yk,* + h/2(fk,* + fk+1,j-1)

Dengan j adalah penghitung iterasi proses koreksi dan

        fk+,j-1 = f(tk+1,yk+1,j-1)

persamaan korektor yang digunakan sebanyak yang diperlukan untuk mendapatkan keakuratan yang diinginkan. perhatikan bahwa dengan menggunakan persamaan Euler sebagai nilai awal,y_k+,j dapat dihitung untuk j=1,2 ....dengan rumus trapesium.proses koreksi dapat dihentikan setelah iterasi ke-n(ditentukan)atau setelah |y_k+1,j+1-y_k+1,j|<€,untuk suatu nilai € yang ditentukan.

Algoritma(Metode Trapesium Euler)

menghitung hampiran penyeleseian masalah nilai awal y'=f(t,y)dengan y(t₀)=y₀ pada [t₀,b].

INPUT:n,t₀,b,y₀,€ dan fungsi f

OUTPUT:(t_k,y_k),k=1,2,..n

LANGKAH-LANGKAH:

1. Hitung h=(b-t₀)/n
2. FOR k=1,2,..n

      hitung f-=f(tk-1yk-1)

      hitung tk=tk-1+h,z0=yk-1+h*f-

REPEAT

(a) Hitung z=z₀+h/2[f^-+f(t_k,z₀)]

(b) Hitung selisih=z-z₀

(c) simpan z₀=z

UNTIL|selisih|<€

simpan y_k=z₀ 3. SELESEI