问题描述
即使我将变量约束为实数，例如this answer，SymPy仍为我提供了复杂的解决方案。
下面是查找单位球面与平面X = 1的交点的示例。应该有1个真实的解决方案（1,0）。
>>> import sympy
>>> x,y,z = sympy.symbols('x,z',real=True)
>>> list(sympy.solve([x**2 + y**2 + z**2 - 1,x - 1],[x,z]))
[(1,-I*z,z),(1,I*z,z)]
这些是正确的，但不是真实的。
我不能使用solveset(…,domain=S.Reals)，因为Solveset似乎不支持方程组。
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）