// Enton2.sce
// Version améliorée de Enton1.sce
// On choisit les couleurs (courbe, entonnoir, asymptote)
// On affiche un titre et une légende
// Tracé économique de l'asymptote avec xsegs

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
xtitle("entonnoir pour alpha="+string(alpha));
leg="entonnoir@moyennes arithmétiques de la LFGN";
plot2d([K' K' K'],[BH' SPN' BB'], [5 2 5], "121",leg);
xsegs([0 length(X)],[EX EX],3); // tracé en vert (3) de l'asymptote