function x=jacobi(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);
D=zeros(n);
for i=1:n
    D(i,i)=A(i,i);
end
N=A-D;
x=b;
for i=1:k
    x=D\b - D\(N*x);
    x
end
eig(D\N)
    
    