// Le peigne des intervalles de confiance pour 
// l'estimation d'une probabilité inconnue
// 
// Saisie des parametres 
n=input('Taille d''echantillon ?'); 
p=input('Rentrer la probabilité p à estimer:');
Mn=zeros(1:100); // initialisation des moyennes arithmétiques
a=zeros(1:100);b=zeros(1:100);couleur=ones(1:100);
xbasc();
plot2d([0 100],[0 1],[-1 -1],"021");// astuce pour fixer le cadre
xsegs([0 100],[p p],[2]); // horizontale p
for i=1:100,
    X=bool2s(rand(1:n,"uniform")<p);Mn(i)=sum(X)./n;
    a(i)=Mn(i)-1.96.*(4.*n).^(-0.5);
    b(i)=Mn(i)+1.96.*(4.*n).^(-0.5);
    if (p<a(i))|(p>b(i)) then 
         couleur(i)=5; // rouge si p hors intervalle
    else couleur(i)=13; //  vert si p dans intervalle
    end,
    xsegs([i i],[a(i) b(i)],couleur(i))  
end    
rep=find(couleur == 5); // numeros des intervalles rouges  
if rep == [] then
  printf('%s\n', 'Tous les intervalles de confiance contiennent p');
else
  printf('\n%s\n',...
         ' p hors intervalle de confiance pour les échantillons No : ');
  printf('%5i\n',rep');
end