Algoritmo Jacobi – Sistemas Lineares – Método Interativo
Disponível na aula MNC_Sistemas Lineares_Aula7
Orientações: Para executar os arquivos com a extensão “sci ou sce”, salve-os no seu computador e siga os passos abaixo:
-
Abrir o Scilab. Em seguida acessar o menu Arquivo->Abrir Arquivo
-
Selecione a pasta onde você salvou os arquivos, busque o arquivo e clique em abrir. O arquivo será aberto em um ambiente denominado SciNotes.
-
Para alterar o caminho padrão do Scilab. Acesse Arquivo->Alterar o diretório corrente. Em seguida, selecione o diretório onde se encontram os arquivos (Sci/Sce).
-
Para executar o arquivo: No console do Scilab digite exec(“nome do arquivo.sci”)
-
Em seguida digite o comando de entrada de dados referente a cada função.
OBS: caso necessite limpar as variáveis digite o comando clear, ou clc para limpar a janela de comandos
Salve o algoritmo abaixo e execute:
exec(“jacobi.sci”)
-> [x, err , niter ] = jacobi ([10 3 -2; 2 8 -1; 1 1 5], [57 20 -4]’, [0 0 0]’, 0.00001 , 50)
function [x , err , niter ] = jacobi (A ,b ,x0 , tol , maxiter )
err = %inf ;
niter = 1;
D = diag ( diag (A) );
T = A - D;
while err >= tol & niter < maxiter ,
x = D \(b -T * x0 );
err = max ( abs (x - x0 ));
x0 = x;
niter = niter +1;
end
if niter == maxiter then
disp ('Numero maximo de iteracoes atingido ');
disp ('iteracoes ', niter , 'Condicao max_k |x(k) - x0(k)| nao atingida apos ');
break ;
end
endfunction