Método de Gauss-Seidel em SCILAB

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

[ Hits: 25.437 ]

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

MoonScript - Agenda telefônica semifuncional em 101 linhas

Adicionar proxy no Internet Explorer na inicialização

impressora brother

trocar permissão, dono e grupo de arquivos ou diretórios

Crivo de Eratóstenes Simples em XBase (Clipper)


  

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