完美数 Matlab实现

%完美数
x=input('数字：');
A=[]; 
for i=1:x
    if mod(x,i)==0  
       A(end+1)=i; 
    end
end
b=numel(A);  %count
A(1,b)=0;  
c=sum(A);
if c==x
    disp(x)
else
    disp('不是完美数')
end