#!/usr/bin/python import math comic_count = 361 data = [float(x) for x in open('numbers.data')] mean = sum(data) / len(data) length = len(data) diff = [(x - mean)**2 for x in data] print 'actual stats', math.sqrt(sum(diff)/length), mean a = 1.0 b = comic_count n = b - a + 1 print 'uniform stats', math.sqrt((n**2 - 1)/12), (a+b)/2