//  entonnoir encadrant S_k/k pour un échantillon d'une Bernoulli,
//  version avec initialisation du générateur de nombres aléatoires
//  par l'utilisateur
germe=input('Initialisation générateur : taper un entier (9 chiffres max)');
n=input('Rentrer la taille n de l''échantillon (n>10) : '); 
// Attention méchant piège avec l'apostrophe de échantillon!!!!
p=input('Rentrer le paramètre p de la loi de Bernoulli : ');
alpha=input('Rentrer le paramètre alpha de l''entonnoir (alpha>1/2): ');
rand('seed',germe);
// On commence l'affichage a 10
// Il faut des vecteurs colonnes pour plot2d 
t=(10:n)';
X=bool2s(rand(t,'uniform')<p); 
S=cumsum(X);
// correction liée au debut à 10
S9=sum(bool2s(rand(1:9,'uniform')<p)); S=S9+S;
f=S./t; //suite des fréquences
enton=sqrt(alpha.*log(t)./t);
y=p.*ones(t)+enton;
z=p.*ones(t)-enton;
// affichage  f en bleu (2), branches entonnoir y et z en rouge (5)
xbasc(); plot2d([t t t], [f y z], [2 5 5])//
xsegs([0 n],[p p],3); // tracé en vert (3) de l'asymptote