% Métodos Numéricos - Engs. Civil e Electrotécnica 
% exercico 2 da ficha 8
% Carlos Balsa - 2010/2011

% alinea a)
a=-10;
b=10;
x=a:0.1:b;
fx=x.*sin(x)-1;
clf
plot(x,fx); hold on
grid on

% alinea b)
% resolução no papel (ver aula)

% alinea c)
fun=@(x) x*sin(x)-1;          % função a anular
dfun=@(x) sin(x)+x*cos(x);    % derivada da função
itmax=10;
tol=10e-10;
x0=-0.5;        % aproximação inicial à solução determinada à partir do grafico
[x1,f1,k]=nle_newtraph(fun,dfun,x0,itmax,tol) 
x0=0.5;        % aproximação inicial à solução determinada à partir do grafico
[x2,f2,k]=nle_newtraph(fun,dfun,x0,itmax,tol)


plot(x1,f1,'*r'); hold on % inclusão da raiz no gráfico
plot(x2,f2,'*r')          % inclusão da raiz no gráfico
legend('Função','Raízes de menor valor absoluto')
