samer@4: function Y=sum1(X)
samer@4: % sum1 - sum over dimension 1 and shiftdim 1
samer@4: %
samer@4: % sum1 :: [[N M]] -> [[M]].
samer@4: 
samer@4: 	Z=sum(X,1);
samer@4: 	S=size(Z);
samer@4: 	Y=reshape(Z,[S(2:end) 1]);