%wlasna diag dominujaca
%A = [   30, 3, 3, 2 ,4 ;  5, 54, 2, 9, 6; 4, 7, 43, 3, 1; 4, 2, 1, 65, 5,; 7, 5, 3, 1, 45;];

%Wlasna diag niedominujaca
%A = [1, 2, 3, 4, 5;6, 7, 8, 9, 9;7, 6, 5, 4, 3;4, 5, 7, 8, 10;5, 4,3 ,2 ,1];

%b = [1, 2,4,5,6]';


x_0 = zeros(length(A), 1);
max_iter = 600;
tol = 1e-6;
omega = 1.1;

tic
%[x_i, i, norma] = jacobi(A, b, x_0, max_iter, tol);

%[x_i, i, norma] = gaussSeidel(A, b, x_0, max_iter, tol);

[x_i, i, norma] = SOR(A, b, x_0, max_iter, tol, omega);

toc