compute probabilities of Savage Worlds deaths

gistfile1.py

import random

def d(n):
    r = random.randint(1, n)
    return r if r < n else r + d(n)

def dw(n):
    return max(d(n), d(6))

def old_rules(n):
    for vigor in 4, 6, 8, 10, 12:
        dead = 0
        for i in xrange(n):
            while True:
                r = dw(vigor) - 3
                if r <= 1:
                    dead += 1
                    break
                if r >= 4:
                    break
        print vigor, dead * 100 / n

def new_rules(spirit, n):
    for vigor in 4, 6, 8, 10, 12:
        dead = 0
        for i in xrange(n):
            if dw(spirit) - 3 < 4 or dw(vigor) - 3 < 4:
                dead += 1
        print vigor, dead * 100 / n