samer@4: % findmd - Multidimensional find
samer@4: %
samer@4: % findmd :: [D:[[1,E]]] -> [[N,E]->natural].
samer@4: %
samer@4: % Take an E dimensional array and return subscripts of nonzero
samer@4: % elements.
samer@4: function S=findmd(A)
samer@4: 
samer@4: E=ndims(A);
samer@4: [SX{1:E}]=ind2sub(size(A),find(A));
samer@4: S=cellcat(2,SX);