求解R中的方程組

Question

象征性地求解方程可以使用Ryacas庫在R中實現。 例如

library(Ryacas)
yacas("Solve(x/(1+x) == a, x)")

給

expression(list(x == a/(1 - a)))

有人知道如何（象征性地）解決方程組嗎？

謝謝。

Answer 1

好吧，我使用優秀的python庫， sympy ，進行符號計算。

使用sympy ，直接求解方程組：

>>> from sympy import *
>>> x,y = symbols('x y')
>>> solve([Eq(x + 5*y, 2), Eq(-3*x + 6*y, 15)], [x, y])
{y: 1, x: -3}

這就是如何使用符號代數來解決方程組，除了通過python包。

好消息是有一個名為rsympy的R端口，可以在CRAN或Google Code上找到 。

除了下載/安裝它並使用rsympy手冊中的幾個最簡單的例子之外，我從未使用過rsympy。 在過去的三年里，我經常使用原來的python庫，我可以高度推薦它。

Answer 2

試試這個：

yacas( "OldSolve({a*x+y==0,x+z==0},{x,y})" )