function m=prodc(x);
% get product of each column
m=prod(x);
if size(m,2)>1;
   m=m';
end;
return;