m   = 7;
n   = 10000;

U   = rand([n, m]) - 0.5;
sim = zeros(m, 3);

for i = 1:n 
    U_sort = sort(U(i,:));
    T      = mean(U_sort([2 3 4 5 6]));
    M      = U_sort(4);
    m      = mean(U_sort);
    
    sim(i, 1) = T;
    sim(i, 2) = M;
    sim(i, 3) = m;
end

hist(abs(sim(:, 2)) - abs(sim(:, 1)))
hist(abs(sim(:, 3)) - abs(sim(:, 1)))