%%
clear
NumDados = 1e3;
x0 = 50;
sigma0 = 5;
InicioHist = 25;
LarguraCanal = 2;
FinalHist = 75;
Bins = ( InicioHist : LarguraCanal : FinalHist);

xi = x0 + sigma0*randn(NumDados,1);

[ n_ocorrencias, x_histograma ] = hist( xi, Bins );
DensidadeFrequencia = (n_ocorrencias/NumDados) / LarguraCanal;
figure
hold on
bar( x_histograma, DensidadeFrequencia, 1, 'g' )
set( gca, 'tickdir', 'out', 'box', 'on', 'fontsize', 14, ...
    'xgrid', 'on', 'ygrid', 'on' )
xlabel( 'X (ux)', 'fontsize', 18 )
ylabel( {'Densidade de Frequencia (ux^{-1})'}, 'fontsize', 18 )
title( sprintf( 'N = %d', NumDados ) )
xlim( [InicioHist FinalHist] )