S=[70 110 65 30 25 40 60] K=[4 2 3 1 5 7 6]; [~,idx]=sort(K); out=S(idx)

There are a couple of ways to achieve this:

[y1,I] = sort(x(1,:))

y2 = x(2,I)

y = [ y1 ; y2 ]

This basically sort the first row of your data, and save the sorting index in I, and then use those index to get the 'sorted' second row, and then just join them

or

y = sortrows(x')'

The ' operator transposes the matrix, which allows you to use sortrows and then use it again to reshape your output matrix.

A = [1 2 3 ; 4 5 6 ; 7 8 9]; idx = [2 3 1] ; A = A(idx,:)

Md, simply use [B,I] = sort(A,'ascend') B = 1 2 3 4 5 I = 3 4 5 2 1

I guess/hope that you made a mistake when you built your example of sorted |t| (that seems to be sorted according to |c=[2;1;3]|). If so, the solution should be: t_sorted = t(c,:) ;

x = randi([0,20], 5, 5); [~, index] = sort(x, 2, 'descend'); % Sorting index result = zeros(size(x)); % Create output result(sub2ind(size(x), 1:5, index(:, 1).')) = 1; % Mark largest result(sub2ind(size(x), 1:5, index(:, 2).')) = 1; % Mark 2nd largest

>> [a, sort_index]=sort([100;20;30;40]) a = 20 30 40 100 sort_index = 2 3 4 1 This output is correct. It tells you that a(1) == 20 came from element #2 of the input, and that a(2) == 30 came from element 3 in the input, and that a(3) == 40 came from element 4 of the input, and that a(4) == 100 came from element 1 of the input. To get the output you want, use b(sort_index) = 1:length(sort_index);

[m,n]=size(A); A=A(sub2ind([m n],idx,repmat(1:n,m,1)))