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)
}
}
}
