metodo de jacobi
malugiru
23 de Abril del 2008
necesiyo urgente el metodo de jacobi
Chris
23 de Abril del 2008
% 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
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