GNU Octave 中符号的任意精度不起作用

Question

我试图检查有趣的计算

x = -80538738812075974
y = 80435758145817515
z = 12602123297335631
x**3+y**3+z**3=42 

与GNU Octave。

我使用了符号包并写道：

>> x = vpa('-80538738812075974')
x = (sym) -80538738812075974.0000000000000
>> y = vpa('80435758145817515')
y = (sym) 80435758145817515.0000000000000
>> z = vpa('12602123297335631')
z = (sym) 12602123297335631.0000000000000
>> vpa(x**3+y**3+z**3)
ans = (sym) -23634890844440363008.0000000000

所以答案不是 42，即使有数字 (50(。使用 wxMaxima，我得到了正确的结果：

x : -80538738812075974;
y : 80435758145817515;
z : 12602123297335631;
x**3+y**3+z**3;
42

所以数字还可以。是否可以在 GNU Octave 中以任意精度工作？

Answer 1

我真的不明白出了什么问题，但只有当所有步骤都以正确的顺序完成(和 sym 而不是 vpa(时，它才明显有效：

setenv PYTHON d:/anaconda3w64/python
pkg load symbolic
syms x y z
x = sym('-80538738812075974')
y = sym('80435758145817515')
z = sym('12602123297335631')
x**3+y**3+z**3
x = (sym) -80538738812075974
y = (sym) 80435758145817515
z = (sym) 12602123297335631
ans = (sym) 42

所以问题终于解决了。