import random

# A population that is clearly NOT normal: a skewed set of values
population = [1, 1, 1, 2, 2, 3, 10, 15, 20, 50]

# Draw many samples of size 30 (with replacement) and record their means
sample_means = []
for _ in range(1000):
    sample = [random.choice(population) for _ in range(30)]
    sample_means.append(sum(sample) / len(sample))

avg_of_means = sum(sample_means) / len(sample_means)
print("Average of the 1000 sample means:", round(avg_of_means, 3))