Gauss Yöntemi Mathlab'da

[ALpMontana]

clear all
clc
A=input('x in katsayılarını giriniz=');%[17 -2 -3;-5 21 -2;-5 -5 22];
B=input('sonucları giriniz=');%[500;200;30];
P=input('x in baslangıc degerlerini giriniz=');%[0;0;0];
N=length(B);
x=P;
makiter=10000;
epsilon=0.05;
for k=0:makiter;
k=k+1;
for j=1:N
if j==1
x(1)=(B(1)-A(1,2:N)*P(2:N))/A(1,1);
elseif j==N
x(N)=(B(N)-A(N,1:N-1)*x(1:N-1))/A(N,N);
else
x(j)=(B(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j);
end
end
hata=abs(norm(x-P));
P=x;
if (hata<epsilon)
fprintf('iterasyon maksimum iterasyondan önce bitti\n')
break
end
x;
end
fprintf('bulunan x ler=%g\n',x)
fprintf('epsilon degeri=%g\n',hata)