Método de Jacobi em SCILAB

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

[ Hits: 18.224 ]

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

Perguntas e respostas com Assembly e NASM

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

Programação para sistemas embarcados em Assembly

Crivo de Eratóstenes Simples em XBase (Clipper)

gitignore para gerenciar dotfiles


  

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