用R如何求积分？？？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

integral(f(x)*dx, a, b)    这个函数在哪个软件包里呢 ？？
  运行的时候提示没这个函数    谢谢！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：Integral Inte int 软件包 如何 软件包

是integrate

应该是这个： grDevices

integarate可以求    grDevices好像不对  都谢谢了

自己写程序求定积分：
area <- function(f, a, b, eps = 1.0e-06, lim = 10) {  #### f是被积函数；a，b是积分端点；eps是积分精度要求；lim见下
fun1 <- function(f, a, b, fa, fb, a0, eps, lim, fun) {
d <- (a + b)/2; h <- (b - a)/4; fd <- f(d)
a1 <- h * (fa + fd); a2 <- h * (fd + fb)
if(abs(a0 - a1 - a2) < eps || lim == 0)
return(a1 + a2)
else {
return(fun(f, a, d, fa, fd, a1, eps, lim - 1, fun)
+ fun(f, d, b, fd, fb, a2, eps, lim - 1, fun))
}
}
fa <- f(a); fb <- f(b); a0 <- ((fa + fb) * (b - a))/2
fun1(f, a, b, fa, fb, a0, eps, lim, fun1)
}

#### 是对分区间的上限，默认值为10，即被积分区间最多被等分为2的10次方个子区间。
使用实例：
f <- function(x) 1/x
quad<-area(f,1,5)