SAS水仙花数

data mywork.waT;
do i = 100 to 999;
a100 = int(i/100);
a10 = mod(int(i/10),10);
a1 = mod(i,10);
sum = a100**3+a10**3+a1**3;
if i = sum then output;
end;
run;
proc print;
run;