Komputasi untuk Sains dan Teknik - Universitas Indonesia

5.3. ITERASI JACOBI 73Mohon diperhatikan untuk mengganti angka 4 pada statemen for i=1:4 dengan angka yanglain disesuaikan dengan jumlah elemen vektor yang mau dihitung norm2-nya.5.2.2 Perhitungan norm-selisihMisalnya kita punya vektor bernama xlama. Lalu ada vektor lainnya bernama xbaru. Normselisih dari xlama dan xbaru dapat dihitung dengan bantuan script berikut ini1 clear all2 clc3 xlama(1,1)=3;4 xlama(2,1)=-2;5 xlama(3,1)=8;6 xlama(4,1)=5;7 xlama %menampilkan elemen vektor xlama89 xbaru(1,1)=9;10 xbaru(2,1)=4;11 xbaru(3,1)=6;12 xbaru(4,1)=1;13 xbaru %menampilkan elemen vektor xbaru1415 n=4; %jumlah elemen vektor1617 %--------menghitung norm2 selisih -------------18 s=0;19 for i=1:n20 s=s+(xbaru(i,1)-xlama(i,1))^2;21 end22 norm2=sqrt(s)23 %----------------------------------------------Cara perhitungan norm-selisih seperti ini akan diterapkan pada kebanyakan metode iterasi.Jadi tolong diingat baik-baik!!5.3 Iterasi JacobiSekarang kita mulai pembahasan tentang metode iterasi untuk menyelesaikan problem sistempersamaan linear. Metode ini berbeda dengan metode-metode yang telah dijelaskan sebelumnya,dimana ia dimulai dari penentuan nilai awal (initial value) untuk setiap elemen vektor x.Kemudian berdasarkan nilai awal tersebut, dilakukan langkah perhitungan untuk mendapatkanelemen-elemen vektor x yang baru.x (baru) = T x (lama) + c (5.5)