compute_xi function

Auxiliary function. Computes combinatorial matrix.