function [A,b]=metgauss(A,b)
% Eliminare gauss
[n,n]=size(A);
b=b(:);
for p=1:n
    poz=p;
    piv=A(p,p);
    for j=p:n
        if(piv<A(j,p))
            piv=A(j,p);
            poz=j;
        end
    end
    x(1:n)=A(poz,:);
    A(poz,:)=A(p,:);
    A(p,:)=x(1:n);
    for i=p+1:n
        c=-A(i,p)/A(p,p);
        A(i,p:n)=A(i,p:n)+c*A(p,p:n);
        b(i)=b(i)+c*b(p);
        
    end
end