>> a=[100:900];
m1=rem(a,10);
m2=rem(fix(a/10),10);
m3=fix(a/100);
f=find(a==m1.*m1.*m1+m2.*m2.*m2+m3.*m3.*m3);
k=[a(f)]
k =
153 370 371 407
>> a=[100:900];
m1=rem(a,10);
m2=rem(fix(a/10),10);
m3=fix(a/100);
f=find(a==m1.*m1.*m1+m2.*m2.*m2+m3.*m3.*m3);
k=[a(f)]
k =
153 370 371 407