Algoritmo Pronto CN

Prévia do material em texto

Algoritmo de implementação da resolução do sistema de equação linear pelo Método de Gauss
% Resolução do sistema de equações lineares pelo método de Gauss
% 3x1+4x2+7x3+20x4=504
% 20x1+25x2+40x3+50x4=1970
% 10x1+15x2+20x3+22x4=970
% 10x1+8x2+10x3+15x4=601
 
 
a=[ 3 4 7 20 504
 20 25 40 50 1970
 10 15 20 22 970
 10 8 10 15 601];
 
% Aplicação do Metodo
[m,n]=size (a);
 for j=1:m-1
 for z=2:m
 if a(j,j)==0
 t=a(j,:);a(j,:)=a(z,:);
 a(z,:)=t;
 end
 end
 for i=j+1:m
 a(i,:)=a(i,:)-a(j,:)*(a(i,j)/a(j,j));
 end
 end
 x=zeros(1,m);
 for s=m:-1:1
 c=0;
 for k=2:m
 c=c+a(s,k)*x(k);
 end
 x(s)=(a(s,n)-c)/a(s,s);
 end
 disp('Metodo da Eliminacao de Gauss:');
Através do método de eliminação de Gauss, implementado em Matlab, temos a resolução de sistemas lineares de uma forma simples. O método permite transformar qualquer matriz em uma matriz triangular superior por pivoteamento não tendo assim pivô nulo o que tornaria impossível sua resolução obtendo um sistema equivalente de resolução imediata, através de regressões sucessivas. O programa desenvolvido possui uma condição de uso que consiste num determinante diferente de zero.