x1=5*rand(1,100);
x2=5*rand(1,100);
x=x1+x2;
num_bin=ceil(log2(length(x))+1);
h=10/num_bin;%si può anche usare il max dei dati
%quindi (max-min)/num_bin
limiti=[0:h:10];
conta=histc(x,limiti);
centri=limiti(1:end-1)+h/2;
%anche conta(end-1)=conta(end)+conta(end-1);
figure;bar(centri,conta(1:end-1));
freq=conta/100;
figure;bar(centri,freq(1:end-1));
