我有以下循环
syms k m
summ = symsum(k^2,k,1,m)
for j=1:100
summ = subs(symsum((m-j)^2*summ,m,j+1,k),k,m);
end
在R2014a,W7,64b上计算它需要大约600秒。 Mathematica中的等效计算需要60秒。结果是相同的但存在差异。所有和的一个属性(在这种情况下,它们都是一百个)是它可以以半好的方式分解。默认情况下,Mathematica不会将其分解,因为没有告知它这样做。但是symsum通过检查总和来做到这一点,我认为这是减速的根源。
你同意吗?如果是,有没有办法告诉symsum不简化输出?如果不是,有没有办法加快速度呢?我真的很想在Matlab中使用Mathematica。