n=3; N=100; A=[4.63 -1.21 3.22; -3.07 5.48 2.11; 1.26 3.11 4.57]; B=[2.22; -3.17; 5.11]; X=zeros(n,1); X0=zeros(n,1); tol=0.001;

for i=1:N for j=1:n X(j)=(B(j)-A(j,1:j-1)X(1:j-1)-A(j,j+1:n)X0(j+1:n))/A(j,j); end e=norm((X-X0),inf); if(e<tol) break; end X0=X; end X

n=4;N=100;A=[4 1 -1 1; 1 4 -1 -1; -1 -1 5 1; 1 -1 1 3];B=[-2; -1; 0; 1];X=zeros(n,1);X0=zeros(n,1);tol=0.00001;w=1.2; for i=1:N for j=1:n X(j)=(B(j)-A(j,1:j-1)X(1:j-1)-A(j,j+1:n)X0(j+1:n))/A(j,j); X(j)=(X0(j)(1-w))+(X(j)w); end e=norm((X-X0),inf); if(e<tol) break; else X0=X; end end X0