def S3d_1(nbsim):
    from random import randint
    from matplotlib import pyplot as plt
    serieEffectifs = [0] * 16
    for n in range(nbsim):
        de1 = randint(1, 6)
        de2 = randint(1, 6)
        de3 = randint(1, 6)
        troisDes = de1 + de2 + de3
        serieEffectifs[troisDes - 3] += 1
    tabloEffec = dict(zip(range(3, 19), serieEffectifs))
    serieFreq = [round(serieEffectifs[i] / nbsim, 4) for i in range(16)]
    tabloFreq = dict(zip(range(3, 19), serieFreq))
    print(tabloEffec)
    print(tabloFreq)
    "-----Le graphique en barres du tableau des fréquences simulées-----"
    x = range(3, 19)
    plt.figure()
    plt.bar(x, serieFreq, color = 'green') ; plt.xticks(x)
    plt.title("DISTRIBUTION SIMULÉE DE LA VARIABLE S3d")
    plt.xlabel("VALEURS DE S3d SOMME DES TROIS CHIFFRES OBTENUS")
    plt.ylabel("FRÉQUENCE SIMULÉE")
    plt.grid()
    plt.show()