matlab求水仙花数代码

以下是MATLAB代码，用于计算3位水仙花数：

for i=100:999
a=fix(i/100); %百位数
b=fix((i-a*100)/10); %十位数
c=i-a*100-b*10; %个位数
if i==a^3+b^3+c^3 %判断是否为水仙花数
disp(i) %输出水仙花数
end
end