Метод Гаусса-Зейделя

clear;
clc;
a=[4 -1 1;2 6 -1;1 2 -3]   
b=[4;7;0]
e=0.1;
n=3;
miter=100;
k=0;
l=0;
for i=1:n
    x(i)=0;
end
    b=0;
        for i=1:n
            s=0;
                      for j=1:i-1
                      s=s+a(i,j)*x(j);
                         end
                      for j=i+1:n
                     s=s+a(i,j)*x(j);
                      end
                        x=(b(i)-s)/a(i,i);
                        d=abs(x-x(i));
                    
                        if d>b
                   b=d
                     end
                     x(i)=x;    
        end
        
            if b<e 
                l=1;
            else 
                if k==m
                    l=2;
                else k=k+1;
                end
            end
            
            if l==1
                x(i)
            else disp('Интегрция расходится');
            end