A = [1 2 3; 4 5 6]; %data
[xx yy zz] = meshgrid(1:2,3:4,5:6); %column-wise linear indices
B = [A(xx(:)) A(yy(:)) A(zz(:))] %extract
More - generalized to n-dimensions and rows:
A = [1 2 3 7; 4 5 6 9]; %data
[nr nc] = size(A);
C = num2cell(reshape(1:(nr*nc),nr,nc),1);
X = cell(nc,1);
[X{1:nc}] = ndgrid(C{:}); %column-wise linear indices
for ii = nc:-1:1
V(:,ii) = A(reshape(X{ii},[],1));
end
