function prob = prob_gaussian(bnode, bnet)
% PROB_GAUSSIAN Compute the probability of a Gaussian node.
% prob = prob_gaussian(bnode, bnet)

%assert(subset(bnode.self, bnet.observed_nodes));

dist = cond_dist_gaussian(bnode, bnet);
self_val = bnet.values(bnode.self);
prob = gaussian_prob(self_val, dist.mu, dist.Sigma);

