% jacobi.m ***** METODO DE JACOBI *****
function x=jacobi(A,b,x0,tol,itermax)
L=tril(A,-1); U=triu(A,1); D=diag(diag(A));
MJ=-D\\\\(L+U); vJ=D\\\\b;
disp(\\\'Los modulos de los valores propios de la matriz de iteración son\\\');
abs(eig(MJ))
iter=1; x=MJ*x0+vJ; xJ=[x];
while (norm(x-x0)/norm(x))>tol & (iter<itermax)
x0=x; iter=iter+1;
x=MJ*x0+vJ; xJ=[xJ,x];
end
disp(\\\'Vector de iteraciones:\\\'); xJ

Responder