2D indexing of 3D array

I'm wondering if there is a better way to achieve the following without the for loop:

A=rand(20,5);
index=[1 2 3; 10 11 12; 15 16 17]';

output=zeros(size(index,1),size(index,2),size(A,2));
for I=1:size(A,2),
    temp=A(:,I);
    output(:,:,I)=temp(index);
end

Thank you

Answer 1

[the cyclist]

% Your input data
A=rand(20,5);
index=[1 2 3; 10 11 12; 15 16 17]';

% Get the dimensions, for convenience and clarity
a2 = size(A,2);
[i1 i2] = size(index);

% Calculate the output
output = reshape(A(index,1:a2),[i1 i2 a2])

Answer 2

[NA]

Thank you for all the contributions. This was very educational and helpful!