ds_inv_simpson function

Compute Simpson Index