m=4;n=3;a=1:m*n;
A=reshape(a,m,n)
A=A'
% Comment recuperer A(1,3), A(2,1) et A(3,2) ?
I=[1,2,3];J=[3,1,2]; % A(1,3)=A(I(1),J(1)), ...
K=sub2ind([n,m],I,J); % ou K=sub2ind(size(A),I,J);
b=A(K)  % => b(1)=A(1,3), b(2)=A(2,1) et b(3)=A(3,2)
A(K)=0  % =>  A(1,3)=0, A(2,1)=0 et A(3,2)=0
