function x=jacobi2(A,b,k)
%risolve un sistema A*x=b usando k passi
%del metodo di jacobi
%ad ogni passo, scrive il valore corrente\
%di x_k sullo schermo

n=length(b);
x=b;
for h=1:k
x_vecchio=x;
for i=1:n
    s=0; %accumulatore
    for j=1:n
        if i~=j
            s=s+A(i,j)*x_vecchio(j);
        end
    end
    x(i)=1/A(i,i)*(b(i)-s);    
end
x
end

    
    
    
    