clear all
close all

lambda = .1;
Nmax = 10000;
Mmax = 1000;

y = rand(Nmax,Mmax);

x = -1/lambda*log(1-y);

s(:,1) = x(:,1);
for i=2:Mmax,
    s(:,i)=sum(x(:,1:i)')'/i;
end;

for i=1:10:Mmax,
    m = 1/lambda;
    vars = 1/i*(1/lambda^2);
    
    [h,bins]=hist(s(:,i),50); %calcola l'istogramma e i centri degli intervalli ma non lo disegna 
    
    %hist(s(:,i),50);
    Delta = bins(2)-bins(1);
    Knorm= sum(h)*Delta ; % costante di normalizzazione (=area dell'istogramma)
    bar(bins,h/Knorm); % disegna un ``barplot'' con l'istogramma normalizzato
    
    hold on
    tmin=m-4*sqrt(vars);
    tmax= m+4*sqrt(vars);
    dt=(tmax-tmin)/1000;
    t=[tmin:dt:tmax];
    fx = 1/sqrt(2*pi*vars)*exp(-1/2*(t-m).^2/vars);
    plot(t,fx,':r','LineWidth',2)
    
    drawnow
    hold off
    pause
end;




