mb function

Compute the Markov blanket