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