Metodo de Newton Rapson

function newton
global fun dfun
fprintf('metodo de newton:\n');
fun=input('ingrese la funcion:\n','s');
x0=input('ingrese el punto inicial:\n');
tol=input('ingrese la tol:\n');
dfun=diff(fun);
it=0;
fprintf(' it x0 x1 x0-x1');
while(it<50) it="it+1;" x="x0;" x1="x0-(eval(fun)/eval(dfun));">fprintf('el procedimiento se completo satisfactoriamente:\n');
break
end
x0=x1;
end
fprintf('la raiz buscada es=%15.9f\n',x1);
ezplot(fun),grid on

Ejemplo:

metodo de newton:

ingrese la funcion:

exp(-x)-log(x)

ingrese el punto inicial:

1

ingrese la tol:

0.01

el procedimiento se completo satisfactoriamente:
la raiz buscada es= 1.309799389