function T = prod_CPT_and_pi_msgs(CPD, n, ps, msgs, except)% PROD_CPT_AND_PI_MSGS Multiply the CPD and all the pi messages from parents, perhaps excepting one% T = prod_CPY_and_pi_msgs(CPD, n, ps, msgs, except)if nargin < 5, except = -1; enddom = [ps n];%ns = sparse(1, max(dom));ns = zeros(1, max(dom));CPT = CPD_to_CPT(CPD);ns(dom) = mysize(CPT);T = dpot(dom, ns(dom), CPT);for i=1:length(ps)  p = ps(i);  if p ~= except    T = multiply_by_pot(T, dpot(p, ns(p), msgs{n}.pi_from_parent{i}));  endend