function  a = pfijo(fun,x0,tol,maxiter)
% Aproxima por el método del punto fijo una raiz de la ecuacion fun(x)=x
%cercana a x0, tomando como criterio de parada abs(fun(x)-x)<tol o la cota sobre
% el numero de iteraciones dada  por maxiter.
%
% Variables de entrada:
%     fun(x): funcion a iterar, se debe introducir con notación simbolica (eg. 'g')
%     x0: estimación inicial para el proceso de iteración
%     tol: tolerancia en error absoluto para la raiz
%     maxiter: maximo numero de iteraciones permitidas
%
% Variables de salida:
%     a: valor aproximado de la raiz
fprintf(1, 'Metodo del punto fijo \n');
incr=1;
iter=1;
while(incr>tol) & (iter<maxiter)
	a = subs(fun,x0);		% Itera la función g
	incr = abs(a-x0);		% Calcula la variación
fprintf(1, 'iter= %i, a= %x0, incr= %e  \n', iter, a, incr);
	iter = iter + 1;		% Cuenta los pasos
	x0 =a;			% Actualiza x0
end;
% Salida