function bel = compute_bel_tabular(bnode, bnet)

n = bnode.self;
[bel, lik] = normalise(bnet.pearl.msg{n}.pi .* bnet.pearl.msg{n}.lambda);


