题主的提供公式的计算，可以用matlab的for循环语句来实现。
Xit=[。。。]; %n个数据
Xjt=[。。。]; %m个数据
for i=1:n
for j=1:m
syms k
Xi=Xit（i）,Xj=Xjt（j）
Div=symsum(abs(Xi^k/Xi-Xj^k/Xj),1,19)
end
注：题主你的公式的下标似乎一样，请确认。
