在 MATLAB 中使用 vpa 处理大数

Question

我正在尝试将数字 2^1000 中的每个数字相加，例如，25346 中每个数字的和 = 2+5+3+4+6 = 20。

我在 Matlab 中编写了如下代码：

clc, clear all, close all,
x=2^1000;
x=vpa(x,400);
sum=0;
while x>0 
    num=mod(x,10);
    sum = sum+num;
    x=floor(x/10); 
end
sum % = 1349, correct answer should be 1366

对于较小的数字，此代码有效。我猜 vpa 命令有问题，但我不知道是什么。

有人有什么想法吗？在 matlab 中是否有更简单快捷的方法？

Answer 1

按照文档中的说明合并这两行：

x=2^1000;
x=vpa(x,400);

到这里

x=vpa(2^1000,400);