Matlab找出1-10000范围内的完数

在Matlab中，我们可以创建一个函数来找出1到10000范围内的所有完全数。完全数是指能够被其除自身以外的所有正因子之和所整除的自然数。以下是编写一个简单函数的步骤：

function [perfectNumbers] = findPerfectNumbers(n)
    perfectNumbers = [];
    for num = 1:n
        sumFactors = 0;
        for i = 1:sqrt(num) % 只需检查到num的平方根即可
            if mod(num, i) == 0
                factor = i;
                if factor ~= num / factor % 排除重复的因子
                    sumFactors = sumFactors + factor;
                    if factor * factor ~= num % 如果因子不是另一个因子的两倍，才加上
                        sumFactors = sumFactors + (num / factor);
                end
            end
        end
        if sumFactors == num % 如果sumFactors等于原数，它是一个完全数
            perfectNumbers = [perfectNumbers; num];
        end
    end
end

% 调用函数并显示结果
perfectNumbers = findPerfectNumbers(10000);
disp(perfectNumbers);