Algoritmo de GaussSeidel

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

function [x] = GaussSeidel(b,A,E,x)
 n = length(b);
 erro = E+1;
 xk = x;
 while(erro>E)
 for i = 1:n;
 soma = 0;
 for j = 1:n
 if j ~= i
 soma = soma + A(i,j)*x(j);
 end
 end
 xk(i) = (b(i)-soma)/A(i,i); 
 end
 erro = norm(A*x-b); 
 end 
endfunction