
function jacobi_script

% definizione matrice A
A=zeros(10);
for k=1:10
    A(k,k)=4;
    if k > 1
        A(k,k-1)=-1;
        A(k-1,k)=-1;
    end
end

% definizione termine noto
b=ones(10,1);

% iterazioni del metodo di Jacobi
x0=zeros(size(b)); toll=10^(-8); nmax=10000;
[xv,step,flag]=jacobi(A,b,x0,toll,nmax);

semilogy(1:length(step),step,'k-');
saveas(gcf, 'figura', 'jpg');