wolffd@0: function marginal = marginal_domain(engine, i)
wolffd@0: % MARGINAL_DOMAIN Return the marginal on the specified domain (belprop)
wolffd@0: % marginal = marginal_domain(engine, i)
wolffd@0: 
wolffd@0: marginal = pot_to_marginal(engine.marginal_domains{i});