Gauss

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

function x=resret (A,b)
 [l,c]=size(A); // dimensão do argumento/tamanha da matriz A.
 for i=l:-1:1
 soma=0;
 for j=i+1:c
 soma=soma+A(i,j)*x(j);
 end
 x(i)=(b(i)-soma)/A(i,i)
 end
endfunction
function [x,E]=Gauss(A,b)//E = matriz aumentada
 [l,c]=size(A);
 E=[A b];
 for i=1:l-1;
 pivo=E(i,i)
 for j=i+1:l
 m=E(j,i)/pivo;
 E(j,:)=E(j,:)-m*E(i,:)
 end
 end
 x=resret(E(:,1:c),E(:,c+1))
endfunction
A= [2 3 -1;4 4 -3;2 -3 1];
b=[5;3;-1]
[x,E]=Gauss(A,b)
E