Решение систем уравнений в R

Question

Решение уравнения символически может быть достигнуто в R с использованием библиотеки Ryacas. Например

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} 

Так вот как решить систему уравнений с использованием символической алгебры, кроме как через пакет питона.

Хорошая новость заключается в том, что есть R порт SymPy, называемый rsympy, который доступен на CRAN или Google Code, here.

Я никогда не использовал rsympy, кроме загрузки/установки и работы с несколькими простейшими примерами в rsympy Manual. Я использовал оригинальную библиотеку python много в течение последних трех лет, и я могу рекомендовать ее высоко.

Answer 2

Попробуйте это:

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