function res=arbre(n,p) // X(m,k) = nombres d'enfants du k-ieme élément de la generation m X=[]; Z=ones(1,n); for m=1:1:(n-1) for k=1:1:Z(m) X(m,k)=grand(1,1,"geom",p); end Z(m+1)=sum((X(m,:))); end res=X; //res=sum(X,2) endfunction