Método de Gauss-Seidel em SCILAB

Publicado por Ariel Galante Dalla Costa (última atualização em 20/01/2014)

[ Hits: 25.652 ]

Homepage: http://arielgdc.wordpress.com

Download metodo_gauss.sce




Solução computacional para o método de solução de sistemas de equações lineares.

  



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;

Scripts recomendados

Insertion Sort em Scilab

Verifica se o link caiu e manda aviso por email - MIkrotik v5.*

Programação para sistemas embarcados em Assembly

Teste de tempo em LUA

Bubble Sort em Go


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts