function Y=sum1(X)
% sum1 - sum over dimension 1 and shiftdim 1
%
% sum1 :: [[N M]] -> [[M]].

	Z=sum(X,1);
	S=size(Z);
	Y=reshape(Z,[S(2:end) 1]);
