Método de Jacobi em SCILAB

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

[ Hits: 13.672 ]

Homepage: http://arielgdc.wordpress.com

Download metodo_jacobi.sce




Solução computacional para sistemas de equações lineares.

  



Esconder código-fonte

//Solução analítica: x = inv(A)*b
A = [3, -1, -1; 1, 5, 1; 1, 1, 7]; //vetor de coeficientes
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;
it=0;
i=0; s=0;
for i=1:n
    for s=1:n
        AA(i,s) = A(i,s)/A(i,i);
    end;
    be(i) = b(i)/A(i, i);
end;

AB = AA;
bc = be;

while p > eps & it<=1000
    x = bc - ((AB-I)*x);
    p = norm((AB*x)-bc);
    it = it+1;
end;

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

Scripts recomendados

Método das Cordas em SCILAB

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

Insertion Sort em Scilab

impressora brother

Programação para sistemas embarcados em Assembly


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner
Linux banner
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts