%metodo de gauss com pivotagem total
Function[x]=gauss_pivot_total(a,b)
n=length(b);
for k=1:n-1
%pivotagem total
max=abs(a(k,k));
t=k;
r==k;
for i=k+1:n
    for m=k+1:n
        if abs(a(i,m))>max
            max=abs(a(i,m);
            t=i;
            r=m;
        end;
    end;
end;

%troca linhas
if (t~=k)
    for j=k:n
        auxt=a(t,j);
        a(t,j)=a(k,j);
        a(k,j)=auxt;
    end;
    auxb=b(t);
    b(t)=b(k);
    b(k)=auxb;
end;

% troca colunas
if(r~=k)
    for j=1:n
        auxr=a(j,r);
        a(j,r)=a(j,k);
        a(j,k)=auxr;
    end;
end;
%triangularizaçao

    for i=K+1:n
        mult=a(i,k)/a(k,k);
        a(i,k)=0.0;
        for j=k+1:n
            a(i,j)=a(i,j)-mult*a(k,j);
        end;
        b(i)=b(i)-mult*b(k);
    end;
    
% substituiçao

x=zeros(n,1);
x(n)=b(n)/a(n,n);
for i=n-1:-1:1
    soma=a(i,i+1:n)*x(i+1:n);
    x(i)=(b(i)-soma)/a(i,i);
end;
return;