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