3、出的结果打印在屏幕上。（4）用三维图像表示在平面区域D上的f(x,y)。（5）每种计算法都要有计算框图，且每种计算法都要编成一个自定义函数.五、程序及其运行结果程序Functiongailvsymsabcden;a=input('输入值a=');b=input('输入值b=');c=input('输入值c=');d=input('输入值d=');e=input('请选择:1用二重积分计算语句计算概率：2等距网格法计算概率；3用分布函数计算概率；4蒙特卡罗法计算概率.5三维图像在D上表示f(x,y)');whilee>0&&
4、e<6ife==1p=erchong(a,b,c,d)endife==2p=wangge(a,b,c,d);endife==3p=fenbu(a,b,c,d);endife==4p=mente(a,b,c,d);endife==5[X,Y]=meshgrid(-3:0.2:3);Z=1/(2*pi)*exp(-1/2*(X.^2+Y.^2));meshz(X,Y,Z);ende=input('请选择:');end%===============================用二重积分计算functionp=erchong(a,b,c,d)symsxy;f0
5、=1/(2*pi)*exp(-1/2*(x^2+y^2));f1=int(f0,x,a,b);%对x积分f1=int(f1,y,c,d);%对y积分p=vpa(f1,9);%================================等距网格法functionp=wangge(a,b,c,d)symsxy;n=100;r1=(b-a)/n;%求步长r2=(d-c)/n;za(1)=a;fori=1:n,za(i+1)=za(i)+r1;end%分块zc(1)=c;forj=1:n,zc(j+1)=zc(j)+r2;endfori=1:nx(i)=unifrn
6、d(za(i),za(i+1));end%随机取点fori=1:ny(i)=unifrnd(zc(i),zc(i+1));ends=0;fori=1:nforj=1:ns=1/(2*pi)*exp(-1/2*(x(i)^2+y(j)^2))+s;%求和endendp=s*r1*r2;p=vpa(p,9)%============================用正态分布的分布函数计算functionp=fenbu(a,b,c,d)symsxy;f0=1/(2*pi)*exp(-1/2*(x.^2+y.^2));%联合密度函数F=int(f0,x);F=int(F
7、,y);%分布函数F=simple(F);F1=subs(F,{x,y},{b,d});%F(b,d)F2=subs(F,{x,y},{a,d});%F(a,d)F3=subs(F,{x,y},{b,c});%F(b,c)F4=subs(F,{x,y},{a,c});%F(a,c)p=F1-F2-F3+F4%P=F(b,d)-F(a,d)-(b,c)+(a,c)%===========================蒙特卡罗法functionp=mente(a,b,c,d)symsxy;N=10000000;%取点数h=0.5;%¸高度x=a+(b-a)*ran
8、d(N,1);%随机生成点y=c+(d-c)*rand(N,1);z=h*rand(N,1);F=1/(2*pi)*exp(-1/2*(x.^2+y.^2));%联合密度函数i=z
