
function demo_primo

% massimo naturale considerato
nmax=2000;

% vettore_primo: la i-sima componente e' "1" se "i" e' primo, "0" 
% altrimenti.
vettore_primo=[];
for n=1:nmax
    vettore_primo(end+1)=primo(n);
end

% primi_precedenti: la i-sima componente e' il numero di primi piu' piccoli 
% di "i".
primi_precedenti=cumsum(vettore_primo);

% stima_primi: la i-sima componente e' una stima del numero di primi piu'  
% piccoli di "i".

nV=1:nmax;
stima_primi=logint(nV);

% plot dei primi e della stima
erroreV=stima_primi-primi_precedenti;
plot(nV,erroreV,'r-','Linewidth',4);

ultimo_errore=erroreV(end);
fprintf('\n \t ultimo errore: %3.0f \n',ultimo_errore);



% >> demo_primo
% 
%  	 ultimo errore:  12 
%
% >> 
