Baixe o app para aproveitar ainda mais
Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original
//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
Compartilhar