www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第3章/li3_19.m

    for m=100:999
     a1=fix(m/100);  %求a的百位数字
     a2=rem(fix(m/10),10);  %求a的十位数字
     a3=rem(m,10);   %求a的个位数字
     if a==a1*a1*a1+a2*a2+a2+a3*a3*a3
         disp(a)
     end
 end