def probaS3d():
    from matplotlib import pyplot as plt
    UnivS3d = tuple(i + j + k for i in range(1, 7) for j in range(1, 7)
               for k in range(1, 7))
    cardUniversS3d = len(UnivS3d)
    listeEff = [0] * 16
    for k in range(cardUniversS3d):
        listeEff[UnivS3d[k] - 3] += 1
    Proba = [round(listeEff[i] / cardUniversS3d, 4) for i in range(16)]
    tabloProba = dict(zip(range(3, 19), Proba))
    print(tabloProba)
    x = range(3, 19)
    plt.figure()
    plt.bar(x, Proba) ; plt.xticks(x)
    plt.title("DISTRIBUTION DE PROBABILITÉ DE LA VARIABLE S3d")
    plt.xlabel("VALEURS DE S3d SOMME DES TROIS CHIFFRES OBTENUS")
    plt.ylabel("PROBABILITÉ")
    plt.grid()
    plt.show()