计算R中的积分

Question

我正在自学R，并且在浏览谷歌示例时，我一直看到提到积分或计算R中的整合。我正在使用的三本书（使用R，R's Cookbook，Beginning R发现统计数据）没有提及计算这些。

我在这里发现了几个问题，但它们与这个问题的格式/方向不同。此外，this问题也没有提供任何信息。

所以我的问题是，你如何计算一个：

我猜h(x)在单位间隔内是统一的吗？

鉴于我只找到了“文本”而没有实际的实现示例，我可能会解雇它并继续前进，如果你们不能验证任何东西。提前谢谢。

Answer 1

我不确定你是否在这里问这个非常简单的问题，但是

f <- function(x) x*(x^5-1)
integrate(f,0,1)
## -0.3571429 with absolute error < 4e-15

似乎工作正常。分析检查：

f_indef <- function(x) x^7/7-x^2/2
f_indef(1)-f_indef(0)
## [1] -0.3571429

或懒惰：

library(Ryacas)
(intres <- yacas("Integrate (x) x*(x^5-1)",addSemi=FALSE))
## x^7/7-x^2/2;

除了Eval.yacas之外，我无法让NULL给我任何东西，但是......