MDISC function

Compute multidimensional discrimination index