Toss a coin "ntoss" times and repeat "experiments" times. Tallies up the observed average outcome and the observed standard deviation

experiments = 10;
AveResult = 0;
AveSquared = 0;
Do[
ntoss = 100000; 
nheads = 
          0; 
Do [
If[Random[] < 0.5, 
            nheads = nheads + 1], 
{n, 1, 
            ntoss}
]; 
nhbar = 
          nheads/ntoss  // N; 
AveResult = 
          AveResult + nhbar; 
AveSquared = 
          AveSquared + 
            nhbar*nhbar; 
Print["Probability(heads) = \
", nhbar], 
{m, 1, 
          experiments}
];
AveSquared = AveSquared/experiments;
AveResult = AveResult/experiments;
Print ["Sigma (nhbar) = ",  
      Sqrt[AveSquared - AveResult^2] // N]