clear all;
clc;
disp(‘Ödev – Ekrem TÜRKAY 020306053 Bilgisayar Öğrt.-3’)
disp(sprintf(‘——————%s’,date,’-‘));
disp(‘——————-SEIDEL Yöntemi——————-‘);
A=input(‘A matrisi=’);
B=input(‘B matrisi=’);
%A=[1 10 2;2 1 10;10 2 1];
%B=[13;13;13];
N=length(B);genis=[A B];
X=zeros(1,N);
%X=input(‘X”in tahmin değerlerini giriniz’);
maxiter=input(‘Iterasyon son degerini giriniz_:’);
TD=input(‘Tolerans Değerini Giriniz_:’);
for p=1:N
[Y,J]=max(abs(genis(p:N,p)));
C=genis(p,:);
genis(p,:)=genis(J+p-1,:);
genis(J+p-1,:)=C;
if genis(p,p)==0
‘A matrisi tekil olduğu için program durduruluyor’
return
end
if genis(p,p)<=sum(genis(p,[1:p-1,p+1:N]))
‘Yakınsama şartı yok’
return
end
end
A=genis(1:N,1:N);
B=genis(1:N,N+1);
for k=1:maxiter
C=X
for j=1:N
X(j)=(B(j)-sum(A(j,[1:j-1,j+1:N]).*X([1:j-1,j+1:N])))/A(j,j)
end
if (abs(X-C)<=TD)
break
end
end
disp(sprintf(‘Iterasyon sayisi=%d’,k));
for n=1:N
eval([‘X’,int2str(n),’=X(n)’])
end
