
function demo1

% Approssimazione di "gamma_5=\int_0^{1000} t^{(x-1)} exp(-t) dt" con
% tolleranza 10^(-6) e paragone con "gamma(5)".

f=@(t) (t.^4).*exp(-t); % integranda 
a=0; b=100; % estremi di integrazione
toll=10^(-10); % tolleranza
[gamma5,flag] = integrale_adattativo (f,a,b,toll);

sol=gamma(5);
err=abs(sol-gamma5);

if flag == 0
    fprintf('\n \t terminazione corretta');
else
        fprintf('\n \t terminazione non corretta');
end
fprintf('\n \t gamma5  : %1.15f',gamma5);
fprintf('\n \t sol     : %1.15f',sol);
fprintf('\n \t err.ass.: %1.4e',err);
fprintf('\n \n');