Método de Jacobi em SCILAB

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

[ Hits: 15.423 ]

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

Cálculo de HugePages - Excel

Criação de .deb do Avidemux

Software via GPO no Logon de Usuário - SAMBA e AD

Tudo-Wifi

Método das Cordas em SCILAB


  

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