Nan Mihai - Jacobi (Metode numerice)

function [x] = jacobi(A,x0,b,ep,maxiter);
%Metode iterative - Jacobi
    [n,n] = size(A);
    m=n;
    x_nou = zeros(m,1);
    x = x0;
    for k = 1:maxiter
        iter = k;
            for i = 1:n
                    s = 0;
                    for j = 1:n
                        if (j ~= i)
                            s = A(i,j) * x(j) + s;
                        end
                    end
                    x_nou(i) = (b(i) - s) / A(i,i);
            end
            if (norm(x_nou - x) / norm(x) ) < ep
                    iter = k;
                    x = x_nou;
                    return;
            end
                x = x_nou;
        end
end