二重积分

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问在R语言里怎么计算变上限(或下限或上下限)的二重积分呢?对于求指定区域(例如x^2+y^2=1)的二重积分又该怎么求呢?

烦请告诉一下可用的函数或可参考的文章,拜托了!

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：二重积分 R语言 拜托了 积分

同问，可是怎么没人理呢

monte carlo? 或者别的数值积分法应该也很多吧。。。。。

cubature包 adaptIntegrate函数
看帮助找具体用法吧。

qoiqpwqr 发表于 2010-12-16 11:53
cubature包 adaptIntegrate函数
看帮助找具体用法吧。

恩，好东西，先收藏了。
I.2d <- function(x) {
  x1 = x[1]
  x2 = x[2]
  sin(4*x1+1) * cos(4*x2) * x1 * (x1*(x1*x1)^2 - x2*(x2*x2 - x1) +2)
}

adaptIntegrate(I.2d, rep(-1, 2), rep(1, 2), maxEval=10000)

UsageadaptIntegrate(f, lowerLimit, upperLimit, tol = 1e-05, fDim = 1, maxEval = 0, absError=0)Arguments

f	The function (integrand) to be integrated
lowerLimit	The lower limit of integration, a vector for hypercubes
upperLimit	The upper limit of integration, a vector for hypercubes
tol	The maximum tolerance, default 1e-5.
fDim	The dimension of the integrand, default 1, bears no relation to the dimension of the hypercube
maxEval	The maximum number of function evaluations needed, default 0 implying no limit
absError	The maximum absolute error tolerated

太牛了，崇拜啊

请问在上面那个例子中x1和x2的上下限是多少？

frederic7 发表于 2011-4-4 10:50
恩，好东西，先收藏了。
I.2d

高手！不错，谢谢！