function x=Gauss_Seidel(A,b)
n=length(A);
x=zeros(n,1);
y=zeros(n,1);
k=0;
while 1
    y=x;
    for i=1:n
        z=b(i);
        for j=1:n
            if j~=i
                z=z-A(i,j)*x(j);
            end
        end
        if abs(A(i,i))<1e-10|k==100
            return;
        end
        z=z/A(i,i);
        x(i)=z;
    end
        if norm(y-x,inf)<eps
            break;
        end
        k=k+1;
    end