Mi blogpost

(Mi blog) entren a mi blog
http://eleinformatica.blogspot.com/

programas jacobi y seidel

function jacobi(A,B,Es)
x=[0 0 0]
u=triu(A,1);
l=tril(A,-1);
d=diag(A);
e=1;
count=0;
while e>Es
x1=(B-(l+u)*x)./d;
e=max(abs(x1-x));
x=x1;
end
disp('la iteracion es:');
disp(count);
disp(x1);
disp('el error es:');
disp(e);


function seidel(A,B,Es)
x=[0 0 0];
D=diag(diag(A));
u=tril(A,1);
l=triu(A,-1);
count=0;
e=1;
while e>Es
x1=(l+D)\(B-u*x);
e=max(abs(x1-x));
count=count+1;
x=x1;
end
disp('la iteraciones:');
disp(count);
disp(x);
disp('el error es:');
disp(e);

soluciones de la practica