Método de Jacobi em SCILAB

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

[ Hits: 18.180 ]

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

Verifica se o link caiu e manda aviso por email - MIkrotik v5.*

Método da Bissecção em SCILAB

SSH NAO NAVEGA

Cálculo de HugePages - Excel

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

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts