function[N]=occurences(X,m)
// X est une matrice dont chaque ligne est la réalisation 
// de variables aléatoires à valeurs
// dans {0,...,m} 
  mx= size(X,'r');
  // taille de l'échantillon 
  N=zeros(mx,m+1);
  for i=1:m+1
	N(:,i)=sum(X==i-1,'c');
  end
endfunction;