function used in method for fast modified band depth (MBD) calculation
combinat(n, p)
n
number of columns in your dataset
p
number of rows in your dataset
Ying Sun and Marc G.Genton