In [9]:

    

import random

def rollDie():
    """returns a random int between 1 and 6"""
    return random.choice([1,2,3,4,5,6])


    

In [13]:

    

def checkPascal(numTrials = 1000000):
    yes = 0.0
    for i in range(numTrials):
        for j in range(24):
            d1 = rollDie()
            d2 = rollDie()
            if d1 == 6 and d2 == 6:
                yes += 1
                break
    print('Probability of losing = ' + str(1.0 - yes/numTrials))


    

In [14]:

    

checkPascal()


    

    




Probability of losing = 0.50839

In [12]:

    

(35/36)**24


    

    
Out[12]:





0.5085961238690966

In [ ]: