判断3位数是否为水仙花数应该是100到999.所以如下：
set talk off
clear
for I=100 to 999
a=int(i/100)
b=int((i-100*a)/10)
c=i-int(i/10)*10
if i=a^3+b^3+c^3
?i
endif
set talk on
return