// on teste K fois le generateur avec n tirages de la loi p1 // sur {0,...,m} à chaque test par rapport a la loi // uniforme p0 function[nb_rejets]=test_generateur_faux(n,K,m,p1) nb_rejets=0; p0=ones(1,m+1)/(m+1); for i=1:K, // simulation des occurences de 1,...,m lors de n simulations No=grand(1,'mul',n,p1([1:$-1])')'; // calcul des p-valeurs alpha=test_chi2(No,p0); // nombre de rejets à 5% nb_rejets = nb_rejets + (alpha<=0.05); end; endfunction;