% Artur Wiśniewski, opcja 1

A=[20,-5,4,-4; 3,10,-3,-1; 2,0,4,0; 2,-1,-2,8]
b=[12; 20; 28; 52f]

n=length(b);
L=zeros(n);
D=zeros(n);
U=zeros(n);

for i=1:n
 for j=1:n
 if i>j
 L(i,j)=A(i,j);
 elseif i<j
 U(i,j)=A(i,j);
 else
 D(i,j)=A(i,j);
 end;
 end
end

N=inv(D);
H=zeros(n);

for i=1:n
 for j=1:n
 if i!=j
 H(i,j)=-A(i,j)/A(i,i);
 end
 end
end;

if (max(sum(abs(H)))<1) || (max(sum(abs(H')))<1)
 x=zeros(4,1);
 for t=0:12
 disp(['Wynik dla iteracji nr ',num2str(t)]);
 x=-N*(L+U)*x+N*b;
 disp(x);
 end
disp("Sumaryczny blad rozwiazania");
 disp(sum(abs((A*x)-b)));
 disp("Sprawdzenie - wynik dokladny");
 disp(inv(A)*b);
else
 disp("Warunek wystarczający uzyskania zbieżności nie został spełniony");
end;