% Tracé des lois fractales en faisant varier le paramètre
clear;
more off;
nr=1000;
prec=100;
base1=2;
clf;
title('Differents cas pour la fonction de repartition de U');
hold on;
%
for p1=.1:.1:.9,
 sprintf('%d réal. indép. de loi de Bernoulli sur {0,1}...\n',\
  prec*nr)
  coefs=repmat(cumprod(ones(prec,1)/base1),1,nr);
 realis=sum(rbinom([prec,nr],base1-1,p1).*coefs);
 [effectifs,classes]=hist(realis,prec*10);
 plot(classes,cumsum(effectifs)/nr);
end
plot(classes,classes);
xlabel('Valeur');
ylabel('Fonction de repartition');
%
figure;
clf;
title('Differents cas pour la fonction de repartition de U');
hold on;
base2=7;
for p2=.1:.1:.9;
 sprintf('%d réal. indép. de loi binomiales sur {0,1}...\n',\
  prec*nr)
 coefs=repmat(cumprod(ones(prec,1)/base2),1,nr);
 realis=sum(rbinom([prec,nr],base2-1,p2).*coefs);
 [E,C]=hist(realis,prec*10);
 plot(C,cumsum(E)/nr);
end

%legend('Loi uniforme',\
% sprintf('Avec des Bernoulli %1.2f',p1),\
% sprintf('Avec des binomiales (%d, %1.2f)',base2,p2),\
% sprintf('Avec des Bernoulli symetriques'),\
% 4);
