Question

SAGE中的函数solve（）返回我求解方程式的变量的符号值。例如：

sage: s=solve(eqn,y)
sage: s
[y == -1/2*(sqrt(-596*x^8 - 168*x^7 - 67*x^6 + 240*x^5 + 144*x^4 - 60*x - 4) + 8*x^4 + 11*x^3 + 12*x^2)/(15*x + 1), y == 1/2*(sqrt(-596*x^8 - 168*x^7 - 67*x^6 + 240*x^5 + 144*x^4 - 60*x - 4) - 8*x^4 - 11*x^3 - 12*x^2)/(15*x + 1)]

我的问题是我需要在其他计算中使用为y获得的值，但我无法将这些值分配给任何其他变量。有人可以帮帮我吗？

Answer 1

（1）您应该访问Sage用户，专家和开发人员的类似Stack Overflow的ask.sagemath.org论坛！ </plug>

（2）如果你想在某些东西中使用solve（）调用的值，那么使用solution_dict标志可能最简单：

sage: x,y = var("x, y")
sage: eqn = x**4+5*x*y+3*x-y==17
sage: solve(eqn,y)
[y == -(x^4 + 3*x - 17)/(5*x - 1)]
sage: solve(eqn,y,solution_dict=True)
[{y: -(x^4 + 3*x - 17)/(5*x - 1)}]

此选项将解决方案作为字典列表而不是方程列表。我们可以像访问任何其他词典一样访问结果：

sage: sols = solve(eqn,y,solution_dict=True)
sage: sols[0][y]
-(x^4 + 3*x - 17)/(5*x - 1)

如果我们愿意，我们可以将其分配给其他人：

sage: z = sols[0][y]
sage: z
-(x^4 + 3*x - 17)/(5*x - 1)

并替换：

sage: eqn2 = y*(5*x-1)
sage: eqn2.subs(y=z)
-x^4 - 3*x + 17

等等。虽然恕我直言上面更方便，你也可以通过.rhs（）获得相同的结果而无需解决方案：

sage: solve(eqn,y)[0].rhs()
-(x^4 + 3*x - 17)/(5*x - 1)

求解方程并使用在其他计算中获得的值 - SAGE

1 个答案: