Scilab - Método de Jacobi

Scilab - Método de Jacobi


DisciplinaCálculo Numérico12.595 materiais251.074 seguidores
Pré-visualização1 página
//A = Matriz dos coeficientes (NxN)
//B = Matriz dos resultados (Nx1)
//Ep = precisão (0<Ep<1)
//it - número máximo de iterações
//Método para resolução de sistemas lineares
function [x]=jacobi(A,B,Ep,it)
 n=length(A(1,:));
 x=zeros(n,1)
 c=ones(n,1)
 k=0
 while k<it &abs(max(x-c))>Ep
 c=x
 k=k+1
 for s=1:n
 for j=[1:n]
 M(j)=A(s,j)*x(j)
 M(s)=0
 end
 soma=sum(M)
 x(s) = (1/A(s,s))*(B(s)-soma);
 end
 end
endfunction