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