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

Bubble Sort em Go

Método das Cordas em SCILAB

Método da Bissecção em SCILAB

Crivo de Eratóstenes Simples em XBase (Clipper)

Nenhum comentário foi encontrado.