wolffd@0: % row sum
wolffd@0: % function Z=rsum(X)
wolffd@0: 
wolffd@0: function Z=rsum(X)
wolffd@0: 
wolffd@0: [N M]=size(X);
wolffd@0: 
wolffd@0: Z=zeros(N,1);
wolffd@0: 
wolffd@0: if M==1,
wolffd@0:   Z=X;
wolffd@0: elseif M<2*N,
wolffd@0:   for m=1:M,
wolffd@0:     Z=Z+X(:,m);
wolffd@0:   end;
wolffd@0: else
wolffd@0:   for n=1:N
wolffd@0:     Z(n)=sum(X(n,:));
wolffd@0:   end;
wolffd@0: end