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

• Download metodo_gauss.sce
• Enviar nova versã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;

Adicionar proxy no Internet Explorer na inicialização

Criação de .deb do Avidemux

Programação para sistemas embarcados em Assembly

Tranposta da matriz em Haskell

Método das Cordas em SCILAB

Nenhum comentário foi encontrado.