Script para MatLab gauss

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [Tsup]=gauss(Ab)
% Rotina para o método de Gauss
% Entrada: Ab = matriz aumentada
% Saída: Tsup = matriz triangular superior
[nL,nc]=size(Ab);
for k=1:nL-1
 if abs(Ab(k,k))<1.0e-7
 disp('Zero na diagonal principal');
 return
 end
 for i=k+1:nL
 mult=-Ab(i,k)/Ab(k,k);
 disp ([mult]);
 for j=1:nc
 Ab(i,j)=Ab(i,j)+mult*Ab(k,j);
 
 end
 end
 disp ([Ab])
end
if abs(Ab(nL,nL))<1.0e-7
 if abs(Ab(nL,nc))<1.0e-7
 disp('SPI');
 return
 else
 disp('SI');
 return
 end
end
Tsup = Ab;