MonteCarlo.py

import numpy as np

rolls = 10000
n = 6  
E = [0] * rolls 
     
for i in range(0,rolls):
  flag = False
  ar = []
  while flag == False:
    rnum = np.random.randint(low=1, high=7, size=1)[0]

    if rnum not in ar:
      ar.append(rnum)
      
    if len(ar)==6:
      flag =True 
      
    E[i] +=1

# print(ar, E, len(E))    
value = sum(E) / len(E)
print(value)