function [U,margs] = prodmeas(P,m,n)
U = 1;
margs = zeros(m,n);
for i=1:n
  Pi = getmargX(P,m,n,i);
  U = tensprod1(U,Pi);
  margs(:,i) = Pi;
end
return