wolffd@0: function marginal = marginal_nodes(engine, nodes, t, add_ev)
wolffd@0: % MARGINAL_NODES Compute the joint distribution on a set of nodes (smoother_engine)
wolffd@0: % function marginal = marginal_nodes(engine, nodes, t, add_ev)
wolffd@0: 
wolffd@0: if nargin < 4, add_ev = 0; end
wolffd@0: 
wolffd@0: marginal = marginal_nodes(engine.tbn_engine, engine.b{t}, nodes, t, add_ev);