function smallT = marg_tableM(bigT, bigdom, bigsz, onto, maximize)% MARG_TABLE Marginalize a table% smallT = marg_table(bigT, bigdom, bigsz, onto, maximize)  % marg_tableM is the same as marg_table.m, but we are sure  % it is not marg_table.cif nargin < 5, maximize = 0; endsmallT = myreshape(bigT, bigsz); % make sure it is a multi-dim arraysum_over = mysetdiff(bigdom, onto);ndx = find_equiv_posns(sum_over, bigdom);if maximize  for i=1:length(ndx)    smallT = max(smallT, [], ndx(i));  endelse  for i=1:length(ndx)    smallT = sum(smallT, ndx(i));  endendns = zeros(1, max(bigdom));%ns(bigdom) = mysize(bigT); % ignores trailing dimensions of size 1ns(bigdom) = bigsz;smallT = squeeze(smallT); % remove all dimensions of size 1smallT = myreshape(smallT, ns(onto)); % put back relevant dims of size 1