function p = prob_node(bnet, n)
% PROB_NODE Compute Pr(observed self | observed parents)
% p = prob_node(bnet, n)

bnode = bnet.bnodes{bnet.equiv_class(n)};
bnode.parents = parents(bnet.adj_mat, n);
bnode.self = n;
p = prob_bnode(bnode, bnet);
 
