function[x]=jacobi_dejen(A,b,x0,tol,itr)
x=x0;
n=size(x,1);
normVal= Inf; 
disp('       i        x1         x2         x3      x4  ')
while normVal>tol
    xold=x;  
    for i=1:n
        sigma=0;      
        for j=1:n         
            if j~=i
                sigma=sigma+A(i,j)*x(j);
            end           
        end     
        x(i)=(1/A(i,i))*(b(i)-sigma);
    end   
    itr=itr+1;
    normVal=abs(xold-x);
    disp([itr  x']);
end
end