wolffd@0: function c = center(x)
wolffd@0: if isempty(x)
wolffd@0:     c = [];
wolffd@0: else
wolffd@0:     c = x - repmat(mean(x),[size(x,1),1,1]);
wolffd@0: end