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