clear all;
% settings
n = 50;
A = diag([3:(n+2)])+diag((-1).^[1:(n-1)],1)+diag((-1).^[3:(n+1)],-1);
b=ones(n,1);

% Jacobi
tol=10^(-8);
max_it=20000;
x0=zeros(n,1);
[x_jac, err_jac, iter_jac]  = jacobi(A, x0, b, max_it, tol);

% Gauss-Seidel
[x_gs, err_gs, iter_gs]  = gauss_seidel(A, x0, b, max_it, tol);

omega_ott = 1.045454545454545;

% SOR
[x_sor,err_sor,iter_sor]=sor(A,x0,b,omega_ott,max_it,tol);

fprintf('\n\t Le iterate del metodo di Jacobi:       %3d',iter_jac)
fprintf('\n\t Le iterate del metodo di Gauss-Seidel: %3d',iter_gs)
fprintf('\n\t Le iterate del metodo SOR:             %3d\n',iter_sor)