function X = jacobi2(A,b,x0,tol,MAX)
[n, n]=size(A);
X=zeros(n,1);
Xprim=zeros(n);
D=zeros(n);
E=zeros(n);
F=zeros(n);
D=diag(diag(A));
E=-(tril(A)-D);
F=-(triu(A)-D);
G=zeros(3,3);
G=diag(diag(A).^(-1))*(E+F);
for i=1:MAX
    Xprim=X;
    X=G*X+diag(diag(A).^(-1))*b;
    if(norm(X-Xprim)<tol)
        break
    end
end
    
X=x0;