Método de Gauss-Seidel em SCILAB

Descrição

Esconder código-fonte

//Verificar sempre se a diagonal principal é dominante sobre a diagonal secundária
A = [3, -1, -1; 1, 5, 1; 1, 1, 7]; //coeficientes do sistema
I = eye(A);
n = 3;
b = [-2, 14, 24]; //vetor de termos independentes
eps=10^(-2); //Precisão
x0 = [0;0;0]; //Solução inicial
p = 10;
x = x0;
x1 = x;
it=0;
i=0; s=0;

for i=1:n
  bc(i) = b(i)/A(i, i);
end;

big = 0;
num = 0;
while p > eps & it<=1000
  for i=1:n 
    num = 0;
    for j=1:n
      if j <> i then
        num = num + A(i,j)*x(j);
      end;
    end;
    x(i) = (1/A(i,i))*(b(i)-num);
  end;
  p = norm((A*x)-b');
  it = it+1;
end;

if(it>999) then
  printf("Não converge nesse intervalo");
  abort;
else
  printf("Sucesso! Total de iteracoes %f", it);
end;