代码:
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)
运行结果如下:
注:水仙花数是指各位数字的立方之和等于该数本身的三位正整数。
代码:
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)
运行结果如下:
注:水仙花数是指各位数字的立方之和等于该数本身的三位正整数。