J'espère ne pas avoir ébréché l'orthographe du nom de Mlle COTZAKIOTIS. Je viens d'essayer l'idée qu'elle a eu lors du TP5. Cette idée est tellement astucieuse qu'il ne faut pas la laisser perdre (Jean DEMARTINI).

L'idée est d'utiliser la fonction grand(N,M,'bin',1,p) qui engendre un tableau NxM de valeurs suivant une loi binomiales constituées d'un somme de 1 seul terme. L'idée est d'exploiter le fait qu'une loi binomiale utilise une somme de processus de Bernouilli qui engendre une variable aléatoire valant 0 ou 1 avec des probabilités valant respectivement (1-p) et p.

x = grand(1,1000,'bin',1,0.2);  // engendre 1000 0s ou 1s non équirépartis
                                // probabilité d'apparition des 1s : 0.2
s = sum(x);                     // on vérifie ...
s
    185                         // ... et oui !