A=[19,-3,2,-1;2,15,4,-1;-2,0,5,-1;4,-6,-10,19]
2 B=[11;21;30;48]
3 n=length(B);
4 L=zeros(n);
5 D=zeros(n);
6 U=zeros(n);
7 for i=1:n
8 for j=1:n
9 if i>j
10 L(i,j)=A(i,j);
11 elseif i<j
12 U(i,j)=A(i,j);
13 else
14 D(i,j)=A(i,j);
15 end;
16 end;
17 end;
18
19 N=inv(D);
20 H=zeros(n);
21 for i=1:n
22 for j=1:n
23 if i!=j
24 H(i,j)=-A(i,j)/A(i,i);
25 end;
26 end;
27 end;
28
29 if(max(sum(abs(H)))<1) || (max(sum(abs(H')))<1)
30 X=zeros(4,1);
31 for m=0:12
32 disp(['Wynik dla iteracji nr ' num2str(m)]);
33 X=-N*(L+U)*X+N*B;
34 disp(X);
35 end;
36 disp("Sumaryczny blad rozwiazania");
37 disp(sum(abs((A*X)-B)));
38 disp(inv(A)*B);
39 else
40 disp("warunek wystarczajacy uzyskania zbieznosci nie zostal spelnionwy");
41 end;