// Complément à l'exo surréservation aérienne
// poly ICP M206 et partiel I.S. 2 avril 2010


p=0.9; q=1-p;
// Calcul des coefficients binomiaux pour n=300 (triangle de Pascal)
C=[1]; // initialisation n=0
for i=1:300, C=[C 0]+[0 C];end
//
for n=301:325, 
    C=[C 0]+[0 C];
   //Proba de refuser du monde à l'embarquement P(S_n>300)
   // P_Refus(1) = P(S_{301}>300) 
   P_Refus(n-300)=sum(C(301:(n+1)).*p.^(300:n).*q.^(n-(300:n)));
end
// meilleur n :
n=300+max(find(P_Refus<0.01))
// affichage des valeurs de P(S_j<=300) sur 5 colonnes
P_OK = matrix(1-P_Refus,5,5)';
M=zeros(10,5);
M(2:2:10,:) = P_OK;
M(1:2:10,:)= matrix(301:325,5,5)';
M
// conversion en tableau LaTeX (enfin presque)
//getf("/home/suquet/Enseignement/Scilab/Tableaux_Latex.sci");
//frmat2latex(M,'%8.6f')
//et on fignole à la main avec emacs