function []=frempirique(X)
    // Cette fonction permet de tracer la fonction de 
    // répartition empirique de l'échantillon X
    // arguments d'entrée :
    // X : échantillon
    Y=-gsort(-X);
    n=prod(size(Y));
    Y=[Y(1),Y];
    plot2d2(Y,(0:n)/n);
endfunction

// exemple
scf(0);
X=linspace(0,6,1000);
plot2d(X,1-exp(-X),2);
Y=grand(1,1000,'exp',1);
frempirique(Y);



function []=diagbaton(X,valeurs,probatheo)
    // Cette fonction permet d'afficher un diagramme en
    // batons pour un echantillon et de superposer le diagramme 
    // en baton théorique
    //arguments d'entree
    //X : echantillon
    //valeurs : vecteur des valeurs pour lesquelles
    //on souhaite afficher le diagramme en batons
    //probatheo : vecteurs des proba theoriques associees
    //aux valeurs contenues dans le vecteur valeur
    probaempi = [];
    for i=valeurs
        probaempi = [probaempi sum(X == i)/length(x)];
    end
    bar(valeurs,[probaempi;probatheo]');
endfunction


//exemple
scf(1);
x=grand(1,1000,"poi",1);
diagbaton(x,0:10,exp(-1)*(gamma(1+(0:10))).^(-1));
