function [s] =tirage_multiple(p,n)
  s=0
  for i=1:n
    if (rand()<p) then
      s=s+1
    end
  end
endfunction