function x = metodo(A, b, k)

x = b;
n = length(b);

for h = 1:k %iterazioni del metodo
    for i = 1:n
       s = 0; %accumulatore
       for j = 1:n
           if j ~= i
               s = s + A(i,j)*x(j);
           end
       end
       x(i) = 1/A(i,i) * (b(i) - s);
    end
    x    %scrive l'iterata corrente
end
