wolffd@0: function [m, pot] = marginal_nodes(engine, query)
wolffd@0: % MARGINAL_NODES Compute the marginal on the specified set of nodes (global_joint)
wolffd@0: % [m, pot] = marginal_nodes(engine, query)
wolffd@0: 
wolffd@0: pot = marginalize_pot(engine.jpot, query);
wolffd@0: m = pot_to_marginal(pot);
wolffd@0: %[m.T, lik] = normalize(m.T);
wolffd@0: %loglik = log(lik);