%求一对亲和数:任一个数的所有因数的和等于另一个数
function sum_prime=func_factor(x)
sum_prime=0;
for i=1:round(x/2)
if(rem(x,i)==0)
sum_prime=sum_prime+i;
end
end
close all;
Factor=zeros(2,100); %构造 2 X 100 的数组
k=1;
for i=1:10000
num_1=func_factor(i); %计算 i 的所有因数的和
if(i==func_factor(num_1)) %若相等,则找到
Factor(1,k)=i;
Factor(2,k)=num_1;
k=k+1;
end
end
%disp(Factor)