#stddev.py
#standard deviation of a bunch of randomly-generated integers

import random
import math

n = int(input("Enter the number of random numbers you want to generate: "))
max_rand = int(input("Enter the max random number: "))

rands = []
for i in range(n):
    rands.append(random.randint(0,max_rand))

display = input("Display the numbers? (y or n): ")
if display == 'y':
    for num in rands:
        print(num,end=" ")
    print()
 
sum_nums = sum(rands)
mean = sum_nums / n    #"mean" is average

sum_sqrd_diffs = 0    #the sum of each datum's difference from the mean, squared 
for num in rands:
    sum_sqrd_diffs += (num-mean)**2

stddev = math.sqrt(sum_sqrd_diffs / n)

print("Sum:",sum_nums,"  Mean:",mean, "  Std.dev.:", stddev)
#randint makes a uniformly-distributed bunch of numbers,
#whose mean is range/2 and whose standard deviation is range/sqrt(12)