Metodo de Newton Modificado

function newtonmodificado
global fun dfun ddfun
fprintf('metodo de newtonmodificado:\n');
fun=input('ingrese la funcion:\n','s');
x0=input('ingrese el punto inicial:\n');
tol=input('ingrese la tolerancia:\n');
dfun=diff(fun);
ddfun=diff(diff(fun));
it=0;
fprintf('it x0 x1 x0-x1');
while(it<50) it="it+1;" x="x0;" x1="x0-((eval(fun)*eval(dfun))/((eval((dfun)^2))-(eval(fun)*eval(ddfun))));"> 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 newtonmodificado:
ingrese la funcion:
exp(x)-x-1
ingrese el punto inicial:
0.5
ingrese la tolerancia:
10^(-5)

it	x0	x1	x0-x1
1	0.500000000	-0.049299708	0.549299708
2	-0.049299708	-0.000398480	0.048901228
3	-0.000398480	-0.000000026	0.000398453
4	-0.000000026	-0.000000014	0.000000012

el procedimiento se completo satisfactoriamente:

la raiz buscada es = -0.000000014