round_deanon_prob.py

#!/usr/bin/env python2

def capture_probability( n_full_set, n_capture_set, rounds):
    return (n_capture_set / n_full_set) ** rounds


masternode_count = 4270

if __name__ == "__main__":
    for percentage in range (5, 100, 5):
        for rounds in range(1,9):
            capture_set = masternode_count * float(percentage)/100
            print "%s %s %s %.12f " % ( masternode_count, percentage, rounds, capture_probability(masternode_count, capture_set, rounds) * 100)