应用统计作业(1)——概率统计知识求解以定积分的近似值

作者: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


 

原文地址:https://www.cnblogs.com/weitter/p/13764620.html