wolffd@0: function marginal = marginal_nodes(engine, query)
wolffd@0: % MARGINAL_NODES Compute the marginal on the specified query nodes (belprop)
wolffd@0: % marginal = marginal_nodes(engine, query)
wolffd@0: 
wolffd@0: assert(length(query)==1);
wolffd@0: marginal = pot_to_marginal(engine.marginal_nodes{query});