// Enton1.sce
// Ce script trace la ligne polygonale des S_k/k où les
// S_k sont les sommes partielles d'un vecteur (aléatoire) X
// Cette ligne polygonale est encadrée par l'entonnoir centré
// sur la droite y=EX. Il faut fournir la valeur de EX
X=input('Rentrer l''échantillon (vecteur ligne)');
EX=input('Rentrer l''espérance');
alpha=input('Donner le paramètre alpha de l''entonnoir');
K=1:length(X); // vecteur des abscisses entières k=1,...,n
SPN=cumsum(X)./K; // SPN: sommes partielles normalisees
enton=sqrt(alpha.*log(K)./K); // forme de l'entonnoir
BH=EX+enton; // branche haute de l'entonnoir
BB=EX-enton; // branche basse
xbasc(); // nettoyage prealable de la fenetre
plot2d([K' K' K'],[BH' SPN' BB']);