计算物理第七题作业报告
PB13203248 陈召远
一、题目
用 Monte Carlo 方法计算如下定积分,并讨论有效数字位数。
二、算法
1、 利用 方法产生区间的随机数
2、 计算函数值
3、 计算函数值的平均值 1/n*(
)
4、 计算积分
=(1-0)*
5、 同理,产生 5 个随机数,分别为[0,7/10],[0,4/5],[0,9/10],[0,1],[0,11/10]区间
6、 计算函数值 6-x^2-y^2-z^2-u^2-v^2
7、 计算函数值的平均值 1/n*(
)
8、 计算积分为 7/10*4/5*9/10*1*11/10*1/n*(
)
9、 输出结果
三、结果及分析
N 取 1000000 时,
=1.045464
误差为 1/
=0.001
所以有效数字为 4 位