function realis=probadis(n,x,p)
% Generateur aleatoire d une loi discrete a support fini
% n>=1 est le cardinal du support de la loi
% p est un vecteur de n nombres reels positifs tel que sum(p)=1
% x est un vecteur de n nombres reels donnant le support de la loi
% Les realisations successives sont independantes
r=rand;
a=0;
b=p(1);
for i=1:n-1,
   if ((r>=a)&(r<b))
      realis=x(i);
      return;
   end
   a=b;
   b=b+p(i+1);
end
realis=x(n);
return;