%%% Teste la loi du premier chiffre sur des échantillons 
%%% de lois fractales sur [0,1] aléatoirement choisies
%%% Marche pas ! il faut prendre des lois sur [0,1] + générales
%%% dont la loi du xieme chiffre est aléatoire et iid ?
clear;
R=zeros(1,10);
nechant=50;
nlois=50;
for i=1:nlois,
  PDF=rdirichlet(1+rint(9));
  E=rfrac(PDF,20);
  EP=find(E>0);
  R(1)=R(1)+size(find(E==0),1);
  J=(firstdigit(E(EP),10,1))';
  R(J+1)=R(J+1)+1;
end
clf
bar(R)
