Simuler pour vérifier

On veut vérifier par simulation le résultat de l'exercice 6, où quatre feux sont verts avec probabilité 0,4 chacun.

import random

def un_trajet(p, n):
    """Compte les feux verts sur un trajet de n feux."""
    verts = 0
    for _ in range(n):
        if random.random() < p:
            verts = verts + 1
    return verts

def simulation(nb_trajets, p, n):
    compte = [0, 0, 0, 0, 0]
    for _ in range(nb_trajets):
        compte[un_trajet(p, n)] = compte[un_trajet(p, n)] + 1
    return [c / nb_trajets for c in compte]

print(simulation(100000, 0.4, 4))

• Le programme contient une erreur qui fausse complètement le résultat. Laquelle ?
• Corriger la ligne fautive.
• Une exécution correcte affiche [0.1290, 0.3462, 0.3452, 0.1538, 0.0258]. Comparer aux probabilités exactes.
• Pourquoi les valeurs affichées ne sont-elles pas exactement les probabilités théoriques ? Que se passerait-il avec 100 trajets seulement ?