function a =newton(fun,x0,tol,maxiter)
% Aproxima por el método de Newton una raiz de la ecuacion fun(x)=0
%cercana a x0, tomando como criterio de parada abs(fun(x))<tol o la cota sobre
%el numero de iteraciones dada  por maxiter. 
%
% Variables de entrada:
%     fun: funcion a calcualr la raiz, se introduce en modo simbolico 'fun'   
%     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 de Newton \n');
dfun=diff(fun);
f=subs(fun,x0); 
df=subs(dfun,x0);
iter=1;
while(abs(f)>tol) & (iter<maxiter)
	a = x0-f/df;		         		% formula de iteracion
	f=subs(fun,a); df=subs(dfun,a);		%Actualiza f y df
     fprintf(1, 'iter= %i, a= %x0,f= %e  \n', iter,a ,f);
	iter = iter + 1;	                           	% Cuenta los pasos
	x0=a;				% actualiza x
end
% Salida