用rSymPy计算泰勒级数

我一直在试用R接口rSymPy到CAS Symphy，它工作得很好。但是，我找不到使用某些更复杂功能的正确语法，例如查找泰勒级数。 例如，我尝试了以下方法：

library(rSymPy)
sympy("var('p')")
#
##### Cannot make this work ???
#
sympy("from sympy.mpmath import *")
xt <- sympy("p=taylor(exp, 0, 10)")

库（rSymPy）
sympy（“var（'p'）”）
#
#####不能让这个工作？？？
#
sympy（“从sympy.mpmath导入*”）
xt似乎没有明确的泰勒级数可用，但级数函数可用。以下代码起作用：

library(rSymPy)
sympy("var('p')")
sympy("var('x')") # or sympy("x = Symbol('x', real=True)")
#
xt <- sympy("p=series(exp(x), x, 0, 10)") # expand about 0 to 10th order

我们可以通过将代码修改为：

library(rSymPy)
sympy("var('p')")
sympy("var('x')") # or sympy("x = Symbol('x', real=True)")
#
xt <- sympy("p=series(exp(x), x, 0, 10)") # expand about 0 to 10th order
# Remove order information
xt0 <- sympy("p.removeO()")
# Test results
x <- 1/3
T1 <- eval(parse(text=xt0)) # Evaluate the result, xt0
T2 <- exp(x)                # The correct value
print(T1-T2)                # Print the error

我希望这对任何希望使用R包rSymPy的人都有帮助

[1] "1 + x + x**2/2 + x**3/6 + x**4/24 + x**5/120 + x**6/720 + x**7/5040 + x**8/40320 + x**9/362880 + O(x**10)"

library(rSymPy)
sympy("var('p')")
sympy("var('x')") # or sympy("x = Symbol('x', real=True)")
#
xt <- sympy("p=series(exp(x), x, 0, 10)") # expand about 0 to 10th order
# Remove order information
xt0 <- sympy("p.removeO()")
# Test results
x <- 1/3
T1 <- eval(parse(text=xt0)) # Evaluate the result, xt0
T2 <- exp(x)                # The correct value
print(T1-T2)                # Print the error

[1] -4.811929e-12