clear
clf()// Nettoie la fenetre graphique
// Sondages
M=[100,100,100,100]; // on fait 4 series de 100 simulations dun sondage
N=[10,100,400,1000]; // taille commune des sondages dans chaque serie.
Et=[0];
p=input('proportion de 1 dans la population')
for j=1:4, // debut de la grande boucle
  A=rand(1,M(1,j)*N(1,j));
  B=A<p; // T si cest oui, F sinon
  C=bool2s(B); // 1 si cest T, 0 sinon
  D=[];
    for k=1:M(1,j),
      D=[D,sum(C([1],[(k-1)*N(1,j)+1:k*N(1,j)]))/N(1,j)];
    end // D est la liste des proportions calculees au cours des M(1,j) sondages.
  S=D;
//----------------------------------------------------------------------------------------------
// Mediane, quartiles, deciles et centiles de S
  Y=trier(S);
  n=M(1,j); // n est la taille de S
// Mediane
  r=modulo(n,2);
  q=(n-r)/2; 
    if r==0 then 
      med=(Y(q)+Y(q+1))/2; 
    else 
      med=Y(q+1); 
    end
// Quartiles
  Quartiles=[];
  for k=1:3
      Quartiles=[Quartiles,Y(1,ceil(n*k/4))];
  end
// Deciles
  Deciles=[];
  for k=1:9
    Deciles=[Deciles,Y(1,ceil(n*k/10))];
  end
// Centiles
  Centiles=[];
  for k=1:99
    Centiles=[Centiles,Y(1,ceil(n*k/100))];
  end
  m=med                  // Mediane
  Q1=Quartiles(1,1)  // Premier quartile
  Q3=Quartiles(1,3)  // Troisieme quartile
  D1=Deciles(1,1)    // Premier decile
  D9=Deciles(1,9)    // Neuvieme decile
  C1=Centiles(1,1)   // Premier centile
  C99=Centiles(1,99) // Quatre-vingt-dix-neuvieme centile
//----------------------------------------------------------------------------------------------
// Diagramme en boite de S
  u=min(S);
  v=max(S);
  et=v-u;
  Et=[Et,et];
// On dispose maintenant de u, C1, D1, Q1, m, Q3, D3, C99, v et de letendue de la serie
// Il reste a tracer le diagramme en boite de S
  a=-0.2-Et(1,j+1)/3+(j-1)/20+sum(Et(:,[1:j]))/3; // pour eviter les chevauchements de diagrammes
  plot([a,a],[u,v])
  plot([a-et/3,a-et/3,a+et/3,a+et/3,a-et/3,a-et/3,a+et/3,a+et/3],[m,Q1,Q1,m,m,Q3,Q3,m])
  plot([a-et/6,a+et/6],[u,u])
  plot([a-et/6,a+et/6],[v,v])
  plot([a-et/9,a+et/9],[C1,C1])
  plot([a-et/9,a+et/9],[C99,C99])
  plot([a-et/12,a+et/12],[D1,D1])
  plot([a-et/12,a+et/12],[D9,D9])
  plot([a],[u-0.1]) // pour aggrandir un peu le cadre
  plot([a],[v+0.1]) // idem
end