Método de Jacobi em SCILAB

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

[ Hits: 15.031 ]

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

Fazendo um ajuste não linear em dados experimentais - FORTRAN 90

Raízes - Fortran

Código "Vírus do Atalho"

Programação para sistemas embarcados em Assembly

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