R中方程组的求解

使用Ryacas库可以在R中实现方程的符号求解。 比如说

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

给予

有人知道如何（象征性地）解方程组吗

---

谢谢。

好吧，我使用优秀的python库sympy进行符号计算

使用Symphy，直接求解方程组：

>>> 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包

好消息是sympy有一个R端口，名为rsympy，可在CRAN或Google代码上获得

除了下载/安装rsympy并在rsympy手册中使用几个最简单的示例之外，我从未使用过rsympy。在过去的三年中，我经常使用原始python库，我可以高度推荐它。

试试以下方法：

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

---

你对以下几点有想法吗？如果我运行

yacas（“OldSolve（{x+5*y==2，-3*x+6*y==15}，{x，y}）”）

这很好，但我不知道为什么x没有计算到最后得到-3。有可能让R计算出精确的结果吗？我目前没有安装带有Yacas的设备来测试我的解决方案，但在查看文档后，我会尝试：

Yacas（“OldSolve（{a*x+y==0，x+z==0}，x）”）

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