// Question 3 : Variable aléatoire binomiale
///////////////////////////////////////////////////////////////////

n=10;p=0.4; // Paramètres de la loi binomiale
N=10000; // Nombre de réalisations de la variable aléatoire générée

// Génère un tableau y de taille (1,N) de réalisations
// d'une variable aléatoire binomiale de paramètres (n,p)
y=sum(rand(n,N)<p,'r');

xbasc();

// Trace le diagramme en bâton de y
[ind, occ] = dsearch(y, 0:n, "d");
xbasc();plot2d3(0:n, occ/N,rect=[-1,0,n+1,max(occ/N)*1.1],nax=[0,n+3,1,9]);

// Trace en bleu les probabilités exactes
plot2d3((0:n)+0.2,binomial(p,n),style=2);