function [m, pot] = marginal_nodes(engine, query)% MARGINAL_NODES Compute the marginal on the specified set of nodes (global_joint)% [m, pot] = marginal_nodes(engine, query)pot = marginalize_pot(engine.jpot, query);m = pot_to_marginal(pot);