wolffd@0: function [marginal, pot] = marginal_family(engine, query)
wolffd@0: % MARGINAL_NODES Compute the marginal on the family of the specified query node (belprop)
wolffd@0: % [marginal, pot] = marginal_family(engine, query)
wolffd@0: 
wolffd@0: pot = engine.marginal_domains{query};
wolffd@0: marginal = pot_to_marginal(pot);