
L'issue d'un lancer de dé cubique peut être simulée à l'aide de la fonction \texttt{randint(1,6)} du module \texttt{random}.

De même, l'issue d'un lancer de dé tétraédrique peut être simulée à l'aide de la fonction \texttt{randint(1,4)}.

Pour effectuer $n$ simulations, on doit faire une boucle itérative \og for \fg. Cela donne par exemple:
\begin{python}
from random import randint

def simulation(n):
    somme = 0
    for k in range(n):
        a = randint(1,6)
        b = randint(1,4)
        somme = somme + a + b # somme précédente + résultat de l'expérience (a+b)

    return somme / n
\end{python}

L'idée ici est de définir une variable \texttt{somme} qui contiendra la somme des résultats de toutes les expériences.
Pour avoir la somme moyenne, il suffit de diviser cette somme totale par le nombre de simulations.

On obtient:
\begin{console}
>>> for i in range(5):
		print( simulation(1000) )

5.77
6.025
6.033
5.934
6.006
\end{console}

Sur 5 simulations de \numprint{1000} expériences, on s'aperçoit que la somme moyenne est d'environ~6.
