MATLAB求100到999之间的水仙花数

水仙花数是指各位数字的立方之和等于该数本身的三位正整数.

x=[];                                                      %空向量

for n=100:999                               

        a=fix(n/100);                                 %求每次循环时n的百位数

        b=rem(fix(n/10),10);                     %求每次循环时n的十位数

        c=rem(n,10);                                %求每次循环时n的个位数

        if(n==a^3+b^3+c^3)                     %判断n值是否等于各位数字立方之和

               x=[x,n];                                  %是水仙花数时放入行向量里

         end

end

disp(x);                                                  %结束循环检测，输出水仙花数