用Julia解一个2变量方程，任何包

我试图解决：

using SymPy
@syms x y
eq = Equality(x^2 + y^3, 31)
solve(eq)

我得到：

Dict（x=>-sqrt（31-y^3））

Dict（x=>sqrt（31-y^3））

有没有办法得到

x=2，y=3

？不需要和辛皮在一起。我正在寻找任何包，可以有工作完成。

---

PS：丢番图不是一个选项，因为立方体尚未实现。

您正在解决一个未确定的系统，但需要整数解，所以我猜您希望将其作为丢番图问题。我不知道如何从Julia那里做到这一点，但直接使用SymPy（在Python中）有一个函数

diophantine

用于解决此类问题，例如：

In [3]: diophantine(x**2 + y**2 - 25, [x, y])
Out[3]: {(-5, 0), (-4, -3), (-4, 3), (-3, -4), (-3, 4), (0, -5), (0, 5), (3, -4), (3, 4), (4, -3), (4, 3), (5, 0)}

但是，您的示例似乎没有实现：

In [4]: diophantine(x**2 + y**3 - 31)
---------------------------------------------------------------------------
NotImplementedError: No solver has been written for cubic_thue.

---

是的，立体派还没有被用于丢番石。