Составить последовательность наибольшей длины, в которой каждое следующее число делится на предыдущее

clc, clear all
a=randi(10,1,5)
[x1 y1]=meshgrid(a,a);
b=rem(x1,y1)==0-eye(length(a));
c=[x1(b==1)';y1(b==1)']

a =
 
     4    10     5     2    10
 
 
c =
 
     4    10    10    10    10    10    10
     2     5     2    10    10     5     2

clc, clear all
a=randi(10,1,15); b=sort(a);
k=length(b); l1=1; l=1;
while k~=1
    c=b(rem(b,b(l))==0);
    if length(c)<=1
        break
    end
    b=c(2:end);
    x(l1)=c(1);
    k=length(b); l1=l1+1;
    if k==1
        x(l1)=c(end);
    end
end
disp(a)
disp(x)

     5     8     9     8     4     6     1     2     2     7     5     2     5     2     1
 
     1     1     2     2     2     2     4     8     8

a=[3,2,8,4,4,1]; % задание вектора чисел
b=sort(a); % сортировка чисел по возрастанию
[A,B]=meshgrid(b);
C=mod(B,A)==0; % формирование логического массива для определения позиций чисел, делящихся друг на друга без остатка
s=sum(C,2); % суммирование построчно количества чисел, делящихся друг на друга без остатка
D=b(C(s==max(s)),:) % искомая последовательность

D=b(find(b.*C(find((s==max(s))>0),:)>0)) % искомая последовательность

a =
 
     5     8     9     8     4     6     1     2     2     7     5     2     5     2     1
 
>> D1=b(C(s==max(s)),:) % искомая последовательность
Index exceeds matrix dimensions.
 
>> D2=b(find(b.*C(find((s==max(s))>0),:)>0)) % искомая последовательность
Error using  .* 
Matrix dimensions must agree.