predict_bn function

Performs inference over a fitted GBN