ds_blau function

Compute Blau's Index