检索函数的积分表达式

Question

Base R 提供了一个函数D()，它输出给定函数的导数表达式。

例如：

f <- expression(x^3+2*x)
D(f, "x") # with respect to x
# 3 * x^2 + 2

是否有类似的函数可以产生函数积分的表达式？

我知道stats::integrate() 会计算给定区间内函数的积分，但据我所知，它不会输出积分的表达式。

Answer 1

正如cmets中提到的，包Ryacas可以做不定积分：

library(Ryacas)
x = Sym('x')
f = expression(x^3 + 2*x)
Integrate(f, x)
## expression(x^4/4 + x^2)

因为differentiation is mechanics, integration is art.