Add dimensions to a matrix

My RHSvec is a 51X21 matrix. kdpolind is 11X51X21. Doing the following:

[RHSval,kprimeind] = max(RHSvec,[],2);

gives me a 51X1 RHSval and a 51X1 kprimeind.

if kprimeind is as follows:

16

20

20

16

20

16

16

then I want to store in kprimeind in kdpolind as

kdpolind(act,1,16)

kdpolind(act,2,20)

kdpolind(act,3,20)

kdpolind(act,4,16)...

I am unable to do this due to dimensions mismatch. Is there a simple way of doing this?

Thanks!

Answer 1

[Matt J]

[RHSval,kprimeind] = max(RHSvec,[],2,'linear');

kdpolind(act,kprimeind)= vector_that_is_51x1