func normalDistributionDemo() { maxTrials := 100 printEvery := 10 fmt.Printf("\nGenerating %v normal samples. The descriptive statistics ", maxTrials) fmt.Printf("converge to the appropriate values: mean→0, variance→1, skew→0, kurtosis→0. \n") var d stats.Stats for i := 0; i <= maxTrials; i++ { y := rand.NormFloat64() d.Update(y) if i != 0 && i%printEvery == 0 { mean := d.Mean() variance := d.PopulationVariance() skew := d.PopulationSkew() kurtosis := d.PopulationKurtosis() fmt.Printf("itr %v: mean (→0.0) = %0.5f, variance (→1.0) = %0.5f, skew (→0.0) = %0.5f, kurtosis (→0.0) = %0.5f\n", i, mean, variance, skew, kurtosis) } } }