%
% Calcolo Numerico 1 -  4 luglio 2019 - ESERCIZIO 3
%
clear all
alpha=sqrt(5)

g=@(x) (x.^2+5)./(2*x); 
h=@(x) 1/2*x*(3-x.^2/5); 

x0=3.4;
epsilon=1e-4;
nmax=100;
[xg,itg]=pfisso(g,x0,epsilon,nmax);
[xh,ith]=pfisso(h,x0,epsilon,nmax);


errore_g=abs(xg-alpha);
errore_h=abs(xh-alpha);




rapp_g=errore_g(end)/(errore_g(end-1)).^2;
rapp_h=errore_h(end)/(errore_h(end-1)).^2;

fprintf('itg=%2d errg=%12.6e rapp_g=%12.9f \n',itg,errore_g(itg),rapp_g);
fprintf('ith=%2d errh=%12.6e rapp_h=%12.9f \n',ith,errore_h(ith),rapp_h);

 

gh=@(x) 0.37*g(x)+0.63*h(x);
[xgh,itgh]=pfisso(gh,x0,epsilon,nmax);
errore_gh=abs(xgh-alpha);
 rapp_gh=errore_gh(end)/(errore_gh(end-1)).^2;
fprintf('itgh=%2d errgh=%12.6e rapp_gh=%12.9f \n',itgh,errore_gh(itgh),rapp_gh);
  
%itg= 5 errg=4.894840e-11 rapp_g= 0.223605937 
%ith= 8 errh=4.147793e-13 rapp_h= 0.671195681 
%itgh= 6 errgh=3.737899e-12 rapp_gh= 0.339898008 