ds_shannon function

Compute Shannon Index