帮我用MATLAB编程一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。求出全部水仙花数。

以下是使用MATLAB编写的程序，用于找出所有的水仙花数：

for num = 100:999
    digit1 = floor(num/100); % 百位数字
    digit2 = floor(mod(num, 100)/10); % 十位数字
    digit3 = mod(num, 10); % 个位数字
    
    if num == digit1^3 + digit2^3 + digit3^3
        disp(num);
    end
end

运行这段程序，将输出所有的水仙花数：

153
370
371
407