D = size(X,2); % dimension