viernes, 4 de diciembre de 2009

Programa 2.1 (Iteracion de punto fijo - pagina 54).

%Programa 2.1 (Iteracion de punto fijo - pagina 54). Aproximacion a una solucion
% de la ecuacion x=g(x) mediante la iteracion p=g(p) realizada a partir
% de una aproximacin inicial p<0>
%
% Datos
% - g es la funcion de iteracion
% - p0 es el punto de partida
% - tol es la tolerancia
% - max1 es el numero maximo de iteraciones
% Resultados
% - k es el numero de iteraciones realizadas
% - p es la aproximacion al punto fijo
% - err es la diferencia entre dos terminos consecutivos
% - P el la sucecion {pn} completa

P(1)= p0;

for k=2:max1
P(k)=feval(g,P(k-1));
err=abs(P(k)-P(k-1));
relerr=err/(abs(P(k))+eps);
p=P(k);
if (err < tol) | (relerr < tol),break;end
end

P=P';

No hay comentarios:

Publicar un comentario