sweet.py

from scipy.stats import norm

BOXES_IN_TEST = 5000

def box(sigma):
    return [norm(loc=310 / 12., scale=sigma).rvs() for i in range(12)]


def trials(sigma_range):
    for sigma in sigma_range:
        boxes = map(sum, [box(sigma) for i in range(BOXES_IN_TEST)])
        outliers = sum(map(lambda weight: 1 if weight > 317 or weight < 303 else 0, boxes))
        outliers /= float(BOXES_IN_TEST)
        if outliers < 0.1:
            print 'Stdev {} is acceptable: {} of samples are out of the interval'.format(sigma, outliers)


if __name__ == '__main__':
    test_range = [x / 100. for x in range(100, 300)]
    trials(test_range)