function y=Gauss_elimination(a,b)
n=length(b);
% eliminação com pivoteamento parcial
for k=1:n-1
    r=k;
    w=abs(a(k,k));
    for j=k+1:n
        if abs(a(j,k))>=w 
            w=abs(a(j,k));
            r=j;
        end
        
    end
    %trocar linhas k e r
    aux1=a(k,:);
    a(k,:)=a(r,:);
    a(r,:)=aux1;
    aux2=b(k);
    b(k)=b(r);
    b(r)=aux2;
    fprintf('Passo %d: pivô a%d%d=%f',k,k,k,a(k,k));
    a
    b
    fprintf('\n\n')
    for i=(k+1):n
        m=a(i,k)/a(k,k);
        a(i,k)=0;
        b(i)=b(i)-m*b(k);
        for j=k+1:n
            a(i,j)=a(i,j)-m*a(k,j);
        end
    end
end
    fprintf('Passo %d: Resolução do sistema triangular superior',n);
    a
    b
    fprintf('\n\n')
    disp('Solução do Sistema')
y=triangular_superior(a,b);