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