Jacobi

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [xnew, it, dx, xhist] = Jacobi(A, b, x1, eps, maxit)
 xhist = x1';
 n = length(b);
 x = x1;
 xnew = x;
 
 for it = 1:maxit
 for i = 1:n
 suml = b(i) - A(i,[1:(i-1),(i+1):n])*x([1:(i-1),(i+1):n]);
 xnew(i) = suml/A(i,i);
 end
 dx = abs(xnew-x);
 err = norm(dx);
 relerr = err/(norm(x)+eps);
 
 x = xnew;
 xhist = [xhist;x'];
 
 if (relerr < eps)
 break;
 end
 end
end