Método de Newton

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [xk, num_inter] =newton(f, df, x0, to1, max_Iter)
 k = 1;
 erro = 10;
 xkm1 = x0;
 %printf(" k = %d, xk = %f\n", k-1, x0)
 while (max_Iter>k && erro> to1)
 xk = xkm1 - f(xkm1)/df(xkm1);
 erro = abs(f(xk));
 printf(" k = %d, xk = %f\n", k, xk)
 k=k+1;
 xkm1 = xk;
 endwhile
 num_inter = k
endfunction