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: % query must be a single node
wolffd@0: 
wolffd@0: if length(query)>1
wolffd@0:   error('can only handle single node marginals')
wolffd@0: end
wolffd@0: marginal = engine.bel{query};