#!/usr/bin/python

import random

def candidates():
    return [random.gauss(100, 15) for x in range(100)]

def average(values):
    return sum(values) / len(values)

def merit(values):
    return sum(values)

for x in range(10):
    sexist_merit = []
    equal_merit = []
    pure_merit = []

    for x in range(1000):
        men = candidates()
        women = candidates()
        people = men + women
        sexist_merit.append(merit(sorted(men)[93:100] + sorted(women)[97:100]))
        equal_merit.append(merit(sorted(men)[95:100] + sorted(women)[95:100]))
        pure_merit.append(merit(sorted(people)[190:200]))

    print "Averages-- pure:%f, equal:%f, sexist:%f" % (average(pure_merit), average(equal_merit), average(sexist_merit))

stephenj@pastinha:~/sharing (master)$ python quotas.py
Averages-- pure:1306.405407, equal:1302.856431, sexist:1296.833211
Averages-- pure:1307.397649, equal:1304.207673, sexist:1298.507569
Averages-- pure:1307.631590, equal:1304.235352, sexist:1298.252491
Averages-- pure:1306.141478, equal:1302.768670, sexist:1297.322975
Averages-- pure:1307.144755, equal:1303.818964, sexist:1297.999458
Averages-- pure:1306.688380, equal:1303.364114, sexist:1297.520212
Averages-- pure:1306.749686, equal:1303.590202, sexist:1298.286172
Averages-- pure:1305.615128, equal:1302.247885, sexist:1296.545001
Averages-- pure:1307.377004, equal:1304.206474, sexist:1298.380275
Averages-- pure:1306.336245, equal:1303.085456, sexist:1297.115914