// on teste K fois le generateur avec n tirages de la loi uniforme // sur {0,...,m} à chaque test. function[nb_rejets]=test_generateur(n,K,m) // K fois n realisations de la loi uniforme sur {0,...,m} nb_rejets=0; p0=ones(1,m+1)/(m+1); for i=1:K, X=grand(1,n,'uin',0,m); // calcul des nombres d'occurences de 1,...,m dans X No=occurences(X,m); // calcul de la p-valeur alpha=test_chi2(No,p0); // nombre de rejets à 5% nb_rejets = nb_rejets + (alpha<=0.05); end; endfunction;