利用循环找水仙花数

一个3位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。

shu=[];
for m=100:999
    m1=fix(m/100);
    m2=rem(fix(m/10),10);
    m3=rem(m,10);
    if m==m1.^3+m2.^3+m3.^3
      shu=[shu,m];
    end
end
   disp(shu)