%
%   Tema d'esame CALCOLO NUMERICO - MILANO 19 settembre 2013 - Esercizio 3
%
clear all;
f=inline('1-x.*exp(1-x)');
x=[-1:0.05:3];
y=f(x);
%plot(x,y,'-.',x,zeros(length(x),1),':');

f1=fcnchk('-exp(1-x).*(1-x)');
x_vera=1;
x0=input('inserisci x0 per Newton e Secanti =');
x1=input('inserisci x1 per Secanti =');

epsilon=input('inserisci tolleranza epsilon=');
nmax=100;

[xN,itN]=newton(f,f1,x0,epsilon,nmax);

[xS,itS]=secanti(f,x0,x1,epsilon,nmax);

fprintf('Newton: it=%2d x(it)=%16.13f \n',itN,xN(end));
fprintf('Secanti:it=%2d x(it)=%16.13f \n',itS,xS(end));

%epsilon=1e-6
%Newton: it=12 x(it)= 0.9992137332227 
%Secanti:it=17 x(it)= 0.9991013685647

%epsilon=1e-8
%Newton: it=15 x(it)= 0.9999016715689 
%Secanti:it=21 x(it)= 0.9998688243401
 
