// Trace pour plusieurs valeurs du paramètre p,
//  de la fdr F_n de U_n=\sum_{i=1}^n 2^{-i}X_i
// les X_i etant iid de loi Bern(p)
// penser à adapter le vecteur de couleurs si on modifie la taille de p
// éviter les couleurs 8 et 34 qui sont celles du fond d'écran
n=13;
p=0.1:0.1:0.9;
q=1-p;
x=0:2^{-n}:1;
F=zeros(length(x),length(p));
for k=1:length(p),
M=[q(k) p(k)];
for i=2:n,
    M=[M.*q(k) M.*p(k)];
end
F(:,k)=[cumsum(M) 1]';
end

xbasc()
plot2d2("onn",x',F, [2 3 4 5 6 5 4 3 2])