实验五 MATLAB蒙特卡罗方法课件.pptVIP

? ? ?;;例5.14 计算 其中D为y= x – 2与y2 = x 所围 ;例5.15 用蒙特卡罗方法计算;实验:蒙特卡罗方法计算体积;function [q,error]=MonteC(L) if nargin==0,L=7;end N=10000; for k=1:L P=rand(N,3); x=2*P(:,1)-1; y=2*P(:,2)-1; z=2*P(:,3); R2=x.^2+y.^2;R=sqrt(R2); II=find(z=Rz=1+sqrt(1-R2)); m=length(II);q(k)=8*m/N; end error=q-pi;;实验任务一：记录L次实验的实验数据及误差;function icecream(m,n) if nargin==0,m=20;n=100;end t=linspace(0,2*pi,n); r=linspace(0,1,m); x=r*cos(t);y=r*sin(t); z1=sqrt(x.^2+y.^2); z2=1+sqrt(1+eps-x.^2-y.^2); X=[x;x];Y=[y;y]; Z=[z1;z2]; mesh(X,Y,Z) view(0,-18) colormap([0 0 1]),axis off;思考题与练习题;5.下面程序绘出二维图