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