你好,可以这样做,具体每行都有解释
function s = method3(n) q = factor(n);%因式分解 uq = unique(q);%得到唯一因子 k = diff(find([1,diff(q),1]~=0));%每个因子的个数 s = prod((uq.^(k+1)-1)./(uq-1));%公式计算因数之和 end