function Y=sumn(X,n)
% sumn - sum over first n dimension and shiftdim n
%
% sumn :: [[N:[[1,L]] M]], L:natural -> [[M]].
Y=shiftdim(sumdims(X,1:n),n);
