clear
Ms = [1:9];
N = 5e4;

figure
nsV = 3
nsH = 3;
for qM = 1 : length(Ms)
  M = Ms(qM);
  S = zeros( N, 1 );
  for j=1:N
    S(j) = sum( TEFE2016_Ex07_gera_x(M) );
  end
  subplot( nsV, nsH, qM )
  passo = 0.05*sqrt(M);
  hist( S, (-(M+passo/2):passo:(M+passo/2)) )
  xlim( (max(abs(S))+passo)*[-1 1] )
  title( sprintf( 'M=%d, curtose=%.2f', M, kurtosis(S) ) )
  drawnow
end