compute_U function

Compute U Statistic