%metodo de gauss com pivotagem total Function[x]=gauss_pivot_total(a,b) n=length(b); for k=1:n-1 %pivotagem total max=abs(a(k,k)); t=k; r==k; for i=k+1:n for m=k+1:n if abs(a(i,m))>max max=abs(a(i,m); t=i; r=m; end; end; end; %troca linhas if (t~=k) for j=k:n auxt=a(t,j); a(t,j)=a(k,j); a(k,j)=auxt; end; auxb=b(t); b(t)=b(k); b(k)=auxb; end; % troca colunas if(r~=k) for j=1:n auxr=a(j,r); a(j,r)=a(j,k); a(j,k)=auxr; end; end; %triangularizaçao for i=K+1:n mult=a(i,k)/a(k,k); a(i,k)=0.0; for j=k+1:n a(i,j)=a(i,j)-mult*a(k,j); end; b(i)=b(i)-mult*b(k); end; % substituiçao x=zeros(n,1); x(n)=b(n)/a(n,n); for i=n-1:-1:1 soma=a(i,i+1:n)*x(i+1:n); x(i)=(b(i)-soma)/a(i,i); end; return;