print MeanOutliers.nummeandevs(None, 90, 2, 3)
print MeanOutliers.nummeandevs(None, 20, 2, 3)
print MeanOutliers.nummeandevs(None, 21, 2, 3)
print MeanOutliers.nummeandevs(None, 23, 2, 3)
print MeanOutliers.nummeandevs(None, 25, 2, 3)
print MeanOutliers.nummeandevs(None, 15, 2, 3)
print MeanOutliers.nummeandevs(None, 20, 2, 3)
print MeanOutliers.nummeandevs(None, 25, 2, 3)
print MeanOutliers.nummeandevs(None, 25, 2, 3)
print MeanOutliers.nummeandevs(None, 25, 2, 3)
print MeanOutliers.nummeandevs(None, 60, 2, 3)
"""
meansum = 0.0
meandevsum = 0.0
count = 0
std = 40
incount = 0
for i in range(1,10000):
val = random.gauss(0,std)
if abs(val) < 2*std:
incount += 1
meansum += 1.0*val
count += 1
meandevsum += abs(1.0*(val - meansum/count))
MeanOutliers.nummeandevs(None, val, 2, 3)
print meansum/count, meandevsum/count
print "incount", incount
print MeanOutliers.nummeandevs(None, 80, 2, 3)
pass
