a = zeros(100,2); 

for j=1:100
  M=1000000;
  Z=randi([1,6],M,5) ; 
  X=sum(Z,2);
  a(j,1) = mean(X); 
  sigma = std(X);
  a(j,2) = std(X)/sqrt(M);
end