求解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库,我可以高度推荐它.