kuusinoppaa.py

# Kuuden nopean heitto
# Mikä on todennäköisyys sille, että heitolla tulee kaikki numerot 1-6?
# Simuloidaan heittoja
# Juhani Kaukoranta 19.6.2018
import numpy as np
print("Mukissa on 6 arpanoppaa, heitetään ja katsotaan kuinka monessa on kaikki numerot 1-6")
heittoja = int(input("anna heittojen lukumäärä "))

'''heitetän nopat'''
A = np.random.randint(1,7,size=(heittoja,6))
oikeita = 0
'''Katsotaan onko heiton nopissa kaikki numerot 1-6'''
for i in range(0,heittoja):
    if np.all(np.in1d([1,2,3,4,5,6],A[i])):
        oikeita = oikeita + 1
print("oikeita rivejä oli ",oikeita," kpl 1000 rivistä, siis suhde = ",oikeita/heittoja)
print("teoreettinen todennäköisyys sille,että heitolla tulee kaikki numerot 1-6, on 6!/6**6 = 0.0154")