1 # simulacion tiradas de dos dados equilibrados 2 # se calcula la probabilidad de que la suma sea menor a 11 3 # teoricamente deberia ser ~0.916 (por LAPLACE) 4 5 from random import random 6 from math import modf 7 8 sim_iter=10000 9 sim_vect=[[random(), random()] for i in range(sim_iter)] 10 11 def num_to_face(n): 12 return int(modf((n*6)+1)[1]) 13 14 n_salio=0 15 for i, j in sim_vect: 16 if num_to_face(i)+num_to_face(j) < 11: 17 n_salio += 1 18 19 print(n_salio/sim_iter)
