function melange_carlo(n,N)
    z=zeros(1,N);
    for i=1:1:N
        z(1,i)=temps_melange(n)
    end
    clf()
    // histogramme empirique
    // histplot(30,z);
    // fonction de repartion empirique
    t=tabul(z,'i');
    t(:,2)=t(:,2)/N;
    plot2d2(cumsum(t(:,2)));
    plot([n*log(n),n*log(n)],[0,1],'r')
endfunction