clear
x0 = 2.6;
s_a = 0.15;
s_s = 0.05;
N = 500;
nREP = 15;
figure, hold on
for qR = nREP : -1 : 1
    x = gera_dados( N, x0, s_a, s_s );
    xm = mean( x );
    plot( qR-0.5+.8*(1:N)/N, x, '*b', 'markersize', 2 )
    plot( [qR-0.5 qR+0.5], [xm xm], '-k', 'linewidth', 3 )
end