 n=input('Taille des échantillons ?');
 R=rand(n,6);             % n réal. indép. de loi unif sur [0,1]
 T1=-2*log(R(:,1));       % n réal. indép. de T4 
 T2=-2*log(R(:,2));       % n réal. indép. de T4 
 T3=-2*log(R(:,3));       % n réal. indép. de T4 
 T4=-3*log(R(:,4));       % n réal. indép. de T4 
 T7=T4;                   % n réal. indép. de T7 
 T5=3*R(:,5)+R(:,6)+1/2;  % n réal. indép. de T5 
 T6=3*R(:,5)-R(:,6)+3/2;  % n réal. indép. de T6 
 % n réal. indép. de T
 T=max(T7,min(max(T1,max(T2,T3)),max(T6,min(T4,T5)))); 
 % moyenne et écart type empiriques de T 
 M=[mean(T),std(T)]; 
 % histogramme de l'échantillon
 hist(T,fix(sqrt(n))); 