wolffd@0: % column sum
wolffd@0: % function Z=csum(X)
wolffd@0: 
wolffd@0: function Z=csum(X)
wolffd@0: 
wolffd@0: N=length(X(:,1));
wolffd@0: if (N>1)
wolffd@0:   Z=sum(X);
wolffd@0: else
wolffd@0:   Z=X;
wolffd@0: end;