function repartition_empirique(X)
  //
  //  trace la fonction de repartition (empirique) de
  //  X un vecteur (ligne ou colonne) contenant k realisations 
  //
  k = length(X)
  Xo = matrix(X,k,1)  // pour etre sur d'avoir un vecteur colonne
                           // de cette facon le code marche dans les 2 cas
  Xo = -gsort(-Xo)    // Xo pour X ordonne
  // on rajoute maintenant deux points a chaque extremite :
  dx = 0.05*(Xo(k) - Xo(1)) 
  Xo = [Xo(1)-dx ; Xo ; Xo(k)+dx]
  // calcul des ordonnees
  y = [0 ; (1:k)'/k ; 1]
  // et le dessin
  plot2d2("onn", Xo , y, 1, "121", "repartition empirique")