用MATLAB编程求出三位数中全部的水仙花数

代码：

m=100:999;
m1=rem(m,10);  %求个位数
m2=rem(fix(m/10),10); %求十位数
m3=fix(m/100);   %求百位数
k=find(m==m1.^3+m2.^3+m3.^3)   %find(一维向量) 得出一维向量的下标序号
s=m(k)

运行结果如下：

           

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