matlab用蒙特卡罗法估算定积分Word格式文档下载.doc
《matlab用蒙特卡罗法估算定积分Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《matlab用蒙特卡罗法估算定积分Word格式文档下载.doc(3页珍藏版)》请在冰点文库上搜索。
b=1;
m=0;
d=1;
c=[];
fori=1:
n
x=a+rand*(b-a);
y=rand*(maxfx-mimfx)+mimfx;
ify<
=x^2
m=m+1;
end
end
fprintf('
积分值是:
%f\n'
m/n*(maxfx-mimfx)*(b-a))
结果是:
0.333516
0.333116
0.334410
0.333132
0.333022
0.332899
0.333604
0.333535
0.333117
0.333312
0.334002
0.332966
0.332274
0.334320
0.333143
0.333144
0.333840
0.332857
0.333487
0.333332
0.333501
0.333358
0.333129
0.332874
0.333495
0.334282
0.333183
0.333594
0.333894
0.333766
0.333126
0.332993
0.333032
0.333090
0.332987
0.333853
0.333412
0.332419
0.333338
0.333622
0.334497
0.333565
0.333759
0.332973
0.332783
0.332990
0.333409
0.332452
0.334463
0.332927
0.332749
0.332269
0.333389
0.332691
0.333038
0.333918
0.333517
0.333182
0.333919
0.333463
0.332823
0.333679
0.332826
>