Método de Jacobi em SCILAB
Publicado por Ariel Galante Dalla Costa (última atualização em 20/01/2014)
[ Hits: 18.739 ]
Homepage: http://arielgdc.wordpress.com
Solução computacional para sistemas de equações lineares.
//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;
MoonScript - Agenda telefônica semifuncional em 101 linhas
Nenhum comentário foi encontrado.
Por que seu __DIR__ falhou ou o "inferno" dos caminhos no PHP
Preparando-se para certificações da LPI através do LPI Lab
Migração de Arch Linux para repositórios CachyOS (Uso de Instruções v3 e v4)
Ativando o modo Quake no XFCE4-Terminal
[Tutorial] Configurando Multimaster no Samba 4 AD (DC02) + Explicação de FSMO Roles
[Resolvido] Google Chrome reclamando de perfil em uso após mudar hostname
Instalando o Tema de Ícones Tela Circle
Copiar Para e Mover Para no menu de contexto do Nautilus e Dolphin
Senhor Einstein tinha razão mesmo! [RESOLVIDO] (9)
Como instalar Warsaw no Gentoo? (1)
Da pra formatar um netbook MGB e usar algum linux nele? É possível usa... (4)









