% Método de Newton
% Ficha nº6, Ex. nº1 b) e c)
% Carlos Balsa - Métodos Numéricos 2007/08
%
tol=0.5e-10;
x=10;  % Aproximação inicial à raiz
erro=1; % Erro inicial tem de ser maior que tol
k=0;
kmax=50; % numero máximo de iterações permitidas

while(erro>tol & k<kmax)
    k=k+1;
    fx=x.*sin(x)-1; % função
    dfx=sin(x)+x.*cos(x);   % derivada
    h=-fx/dfx;
    x=x+h;
    erro=abs(h);
end
    k      % nº total de iterações
    x      % raiz aproximada
    erro   % erro absoluto