作者:Peter_831
参考链接:概率统计知识求解以定积分的近似值
来源:CSDN
题目:请使用概率统计知识求解以下积分的近似值,并尝试整理成一种通用的方法:
R语言实现:
f=function(x,y,z) x*exp(-y*y+z*z)+x*y*z x1=3;x2=6 y1=2;y2=4 z1=1;z2=3 arr=array(1:1000) v=(x2-x1)*(y2-y1)*(z2-z1) for(i in 1:1000) { r1=runif(10^5,x1,x2)#产生10^5个区间范围内的随机数 r2=runif(10^5,y1,y2) r3=runif(10^5,z1,z2) s=f(r1,r2,r3) out=v*mean(s) arr[i]=out } print(mean(arr))
输出:
[1] 324.009