
function stirling

% qualita' dell'approssimazione del fattoriale mediante formula di Sterling

err=[];

for n=1000
    fattoriale_n=gamma(n+1);
    stirling_n=sqrt(2*pi*n)*(n/exp(1))^n;
    err_n=abs(fattoriale_n-stirling_n);
    err_n=err_n/abs(fattoriale_n);
    err(end+1)=err_n;
end

% salvataggio su file
fileID=fopen('risultati.txt','w');
fprintf(fileID,'%1.15e \n',err);
fclose(fileID);

% figura (spessore=2)
semilogy(1:100,err,'r-','LineWidth',2)
