xbasc();
clear 

n=1000;
Y=rand(n,1); // Renvoie un vecteur colonne dont les entrées sont n
                // réalisations pseudo-aléatoires  selon la loi uniforme
                // sur [0,1] et pseudo-indépendantes 

moyenne=sum(Y)/length(Y) //  Moyenne arithmétique =
                         //  (Y(1)+...+Y(n))/n
						 
mediane=median(Y) // Médiane
												 
plot2d([1:n], Y) // Tracé du vecteur Y. Dans ce cas [1:n] est facultatif

xset('window',1)						 
histplot(20,Y) // Histogramme à 20 classes