用MATLAB编一个求10万以内的所有亲和数对的程序

以下是用 MATLAB 编写的求 10 万以内所有亲和数对的程序：

function amicable_numbers(n)
    for i = 1:n
        sum1 = 0;
        for j = 1:i-1
            if mod(i,j) == 0
                sum1 = sum1 + j;
            end
        end
        sum2 = 0;
        for j = 1:sum1-1
            if mod(sum1,j) == 0
                sum2 = sum2 + j;
            end
        end
        if i == sum2 && i < sum1
            fprintf('%d and %d are amicable numbers\n', i, sum1);
        end
    end
end

你可以调用这个函数并输入 100000 来得到 10 万以内的所有亲和数对。