求解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})" )