%
% Tema d'esame 20 febbraio 2019 - Esercizio n.1 - Matematica
% 
clear all
f=@(x) x-exp(-x);
h=@(x) exp(-x);
h1=@(x) -h(x);
 
g=@(x) h(x)-h1(x)./(1-h1(x)).*f(x);

x0=1;
epsilon=1e-8;
nmax=100;

[x,it]=pfisso_mod(g,f,x0,epsilon,nmax);
g1_alpha=(g(x(end))-g(x(end-1)))/(x(end)-x(end-1));
f_alpha=f(x(end));
fprintf('N=%2d a=%12.8f G1=%12.8f f=%12.6e\n',it,x(end),g1_alpha,f_alpha);

%N= 4 a=  0.56714329 G1=  0.00002828 f=-6.927809e-09