Jakobijanov postupak

% Jakobijanov postupak izracunavanja sistema linearnih jednacina:
function x = jakobi(A,b,x0,itMax,errMax)
[n,m] = size(A);

x = x0;
for it = 1:itMax
    for i = 1:n
        s = 0;
        for j = 1:i-1
            s = s+A(i,j)*x0(j);
        end
        for j = i+1:n
            s = s+A(i,j)*x0(j);
        end
        x(i) = (b(i)-s)/A(i,i);
    end
    err = 0;
    for i=1:n
        err = err + ((x(i) - x0(i))^2);
    end
    if (sqrt(err) < errMax)
        break; % ??
        % errMax = sqrt((x(i) - x0(i))^2);
    end
    x0 = x;

% Poziv:
%  A = [4 1 3; 2 5 6; 1 2 9]
% b = [1; 4; 7]
% errM = 0.01;
% itrM = 100;
% x0 = [1;1;1];
% x = jakobi(A, b, x0, itrM, errM)

% Provjera Gausovom metodom:
% A\b;
end