function [ x ] = MetGaussPivPartiala( A, b )
    n = length(b);
    Aext = [A b]; % Concatenare de matrice

    for k = 1:n-1
        p = k;
        max = abs(Aext(k,k));
        for c = k+1:n
            if abs(Aext(c,k)) > max
                max = abs(Aext(c,k));
                p = c;
            end
        end

        if max == 0
            fprintf 'Sistem incompatibil sau sistem compatibil nedeterminat.';
            return;
        end

        if p ~= k
            Aext([p,k],:) = Aext([k,p],:);
        end

        for l = k+1:n
            m(l,k) = Aext(l,k) / Aext(k,k);
            Aext(l, :) = Aext(l, :) - m(l,k) * Aext(k, :);
        end
    end

    if Aext(n,n) == 0
        fprintf 'Sistem incompatibil sau sistem compatibil nedeterminat.';
        return;
    end

    x = SubDesc(Aext(:,1:n), Aext(:, n+1));
end