Numberphile guesser

guesser.py

import random

n = 10000

def randomize():
    return random.normalvariate(0, 100)

def guess():
    a = randomize()
    b = randomize()
    assert a != b

    # Choice artefact
    k = randomize()

    if a>k:
        # We'll suppose A is the largest
        return a > b
    elif a<k:
        # We'll suppose A is the smallest
        return a < b

success = 0
count = 0
for i in range(n):
    try:
        if guess():
            success += 1
        count += 1
    except:
        continue

print "guesses : %s" % count
print "successes : %s (%s%%)" % (success, (success*100./n))