13.07.2015 Views

Matemáticas para la Computación

Matemáticas para la Computación

Matemáticas para la Computación

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Matemáticas</strong> <strong>para</strong> <strong>la</strong> computaciónDept. d’InformàticaProblemas Universitat de València%* Salida: x : Matriz columna de soluciones. *%* conv : Convergencia alcanzada. *%* k : Número de iteraciones. *%****************************************************************function [x,conv,k] = seidel(a,b,Delta,M)% Dimension del sistema de ecuaciones.fi<strong>la</strong>s=size(a,1);% Inicialización de los vectores columna r y x.r=zeros(fi<strong>la</strong>s,1);x=zeros(fi<strong>la</strong>s,1);u=zeros(fi<strong>la</strong>s,1);for k=1:M% Contador de iteracionesfor i=1:fi<strong>la</strong>sx(i)=b(i);for j=1:fi<strong>la</strong>sif (i ~= j)x(i)=x(i)-a(i,j)*x(j);end;end;x(i) = x(i) / a(i,i);end;r=u-x;u=x;% Calculo de <strong>la</strong> convergencia.conv = norm(r) / norm(x);if ( conv

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!