从 Ryacas 计算表达式中提取 x^k 的系数

extract coefficient of x^k from Ryacas calculated expression

我得到了这样一个Ryacas表达式,我想提取x^4,11的系数。我应该如何编写代码才能得到这个,谢谢。

希望对您有所帮助:

x <- Sym("x")
P <- (x+1)*(x+2)*(x+3)
Expand(P)
# expression(x^3 + 6 * x^2 + 11 * x + 6)
# coefficient of x^2:
yacas("Coef((x+1)*(x+2)*(x+3), x, 2)")
# expression(6)
# or:
yacas(paste0("Coef(", P, ",x , 2)"))
# expression(6)

编辑 2020-02

已针对 Ryacas 的新版本更新:

library(Ryacas)
x <- ysym("x")
yac_assign((x+1)*(x+2)*(x+3), "P")
### expand the polynomial
yac_str("Expand(P)")
# [1] "x^3+6*x^2+11*x+6"
### extract coefficient of x^2:
yac_str("Coef(P, x, 2)")
# [1] "6"