Método de Newton en Matlab

Encuentra la raíz de f iterando un punto p: p=f(p)/f'(p).

function [p]=newton(fun,dfun,a,n,tol)

for h=1:n
d=(feval(fun,a)/feval(dfun,a));
fprintf ('%d: d: %.10f \n',h,d);
if (feval(fun,d))==0
'se ha llegado a la raíz'
break;
end;
if abs(a-d)<tol
'se ha llegado a la tolerancia.'
break;
end;
a=d;
end;

p=d;

if h==n
fprintf ('se ha llegado al limite de iteracciones %d',n);
end;