用R编程求二重积分

owenqi 发表于 2019-9-30 11:35
蒙特卡洛方法是这样的

请问如果我想写x是(y,100),y是(0,100)的区间可以这么写吗？
n <- 10000
x<-runif(n);
y<-runif(n);
y1<--log(y);
x1<--log(x);
z=sum((2*exp(-(2*x1+y1))/(x*y))[x1<y1])/n;
我验证了一下结果差点，请问错在哪里了？

nongyanhua 发表于 2013-1-26 10:08
刚开始接触R，实在不知从何下手，请各位大侠帮帮忙，谢谢啦

学习了