利用MATLAB求自守数，方法一：循环

什么是自守数？
一个数的平方的末尾几位数字等于它本身，如55=25,66=36…

tic
f=zeros(0,100);
index=0;
for p=5:100000
    w=1+floor(log10(p));  %求解这个数的位数
%为什么用floor不用ceil，请自己运行ceil（log10（10））便知
    m=mod(p^2,10^w);
    if p==m;
        index=index+1;
        f(index)=p;
    end
end
an=f(1:index)
toc