clear

n = 1:20;
n2 = linspace( 1, 20, 20 );

% calcolo il fattoriale di n ( n! )
nf = factorial(n);

% approssimazione di Stirling
Sn = sqrt( 2*pi*n ) .* ( n/exp(1) ).^n;

%   calcolo l'errore assoluto
Eass = abs ( Sn - nf );

%   calcolo l'errore relativo
Erel = Eass ./ nf ; 


%   apro la figura 1 in cui voglio inserire il grafico 
figure(1)

%   cancella eventuali grafici precedenti
clf
%   tiene sulla figura tutti i plot
hold on
%   traccia la linea del grafico
plot( n, Erel, 'k-' )
%   rappresenta con simboli i puntoi calcolati
plot( n, Erel, 'r*' )

%   introduco le etichette degli assi
xlabel( '\itn', 'FontSize', 16)
ylabel( 'sin({\itE_{rel}})', 'FontSize', 16, 'Rotation', 0)

%   introduco titolo
title('{\bfErrore relativo per l''approssimazione di Stirling}')


%   attivo il reticolo di rifermento
grid on

%  NON uso assi monometrici
%axis equal
%   restringe gli assi alla prozione disegnata
axis tight
