Jacobi An
Jacobi An
clc;
% Judul program
disp('METODE iterasi jacobi KEL-7')
disp('===========================')
%input matriksss
A=[1.6 2.5 3 6 1;
3 3 2 4 1;
1 2 3 3 2];
[m,n]=size(A);
disp('matriks A');
disp(A);
%masukkan H
H=[2.1;2.3;2;2;2];
[p,q]=size(H);
% masukkan x
X=[5 5 5 9 2];
maks_iterasi=8;
% maksudnya 1.10^-5=0,00001
Eps=1e-5;
%Langkah pengerjaan
for iterasi=1:maks_iterasi
galat=0;
for ii=1:m
Xb=H(ii);
for j=1:n
if j~=ii
Xb=Xb-A(ii,j)*X(j);
end
end
Xb=Xb/A(ii,ii);
selisih=abs((Xb-X(ii))/Xb);
if selisih > galat;
galat=selisih;
end
X(ii)=Xb;
end
if galat < Eps
disp('Hasil adalah');
Xb
break
end
end
if(iterasi==maks_iterasi);
if(galat>Eps);
disp('Hasil adalah');
X
disp('==========000=========')
end
disp('selisihnyaaaaa adalah');
disp(selisih);
end