Matlab：水仙花数

>> a=100:1:999;              %定义矩阵
>> b1=rem(a,10);             %求取矩阵a中每个元素的个位，并形成和a一样的同型矩阵
>> b2=rem(fix(a/10),10);     %求取矩阵a中每个元素的十位，并形成和a一样的同型矩阵
>> b3=fix(a/100);            %求取矩阵a中每个元素的百位，并形成和a一样的同型矩阵
>> k=find(a==b1.*b1.*b1+b2.*b2.*b2+b3.*b3.*b3)   %个位矩阵每个元素的三次方 加上 十位矩阵每个元素的三次方 加上 百位矩阵每个元素的三次方，获得一个总矩阵与a矩阵比较。将相等的序号为保存到k中。

k =

    54   271   272   308

>> c=a(k)            %通过k序号位，获得对应数据

c =

   153   370   371   407