function [xn,erro,k]=Jacobi(A,b,x0,tol,maxiter)
    % n e a dimensao da matriz
    n=length(A);
    % definir um matriz coluna Xn de zeros 
    xn=zeros(n,1);
    %inicialmente o erro tem que ser maior que tol
    erro(1)=tol+1;
    % K e igual a 2 pois o erro so e calculado a partir de xn(2)
    k=2;
    while erro(k-1)>tol & k<maxiter
        for i=1:n
            s=0;
            % soma da primeira parcela da matriz dividida por x0
            for j=1:i-1
                s=s+A(i,j)*x0(j);
            end
            % soma da segunda parcela da matriz dividida por x0
            for j=i+1:n
                s=s+A(i,j)*x0(j);
            end
            % achar o vector xn
            xn(i)=(b(i)-s)/A(i,i);
        end
        % norma infinito
        erro(k)=norm(xn-x0,inf)/norm(xn,inf);
        % o x0 passa a ser o xn
        x0=xn;
        % acrescentar 1 unidade a k
        k=k+1;
    end
return