optAdjSet function

Compute the optimal adjustment set