Método de Gauss - Matlab

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

% Este programa calcula a matriz triangularizada a partir da matriz aumentada,
% usando eliminação gaussiana.
% Fornece também os valores das variáveis de sistema representados
% na matriz aumentada.
% Saídas: a (matriz triangularizada)
% x (vetor das variáveis do sistema)
aux=0;
[a,r,c]=inputaug(aux); %nenhum argumento é transmitido à função
for i=2:r
 [a,f]=pivot(a,i,r,c);
 if f==0
 for j=i:r
 m=a(j,i-1)/a(i-1,i-1);
 a(j,i-1)=0; %faz elemento igual a zero (eliminacao gaussiana)
 
 for k=i:c
 a(j,k)=a(j,k)-m*a(i-1,k);
 end
 end
 end
end
x(r)=a(r,c)/a(r,c-1); %primeira substituicao
for i=r-1:-1:1
 x(i)=a(i,c);
 for j=c-1:-1:i+1
 x(i)=x(i)-a(i,j)*x(j);
 end
 x(i)=x(i)/a(i,i);
end
disp('Matriz triangularizada:')
a
disp('Valores das variaveis do sistema:')
x