Python 有界符号变量

我想通过将多元优化问题转换为非线性方程组，并使用Symphy的“solve”函数进行求解，从而解决多元优化问题，如下所示：

xopt = sympy.solve(gradf, x, set=True)

问题是，这个特殊的方程有无穷多个解，调用solve会冻结我的计算机

如果我可以为符号变量设置上下限，即引入一组约束：

l1