function Y=prod1(X)
% prod1 - prod over dimension 1 and shiftdim 1
%
% prod1 :: [[N M]] -> [[M]].

	Z=prod(X,1);
	S=size(Z);
	Y=reshape(Z,[S(2:end) 1]);
