实验五维维安尼体.pptVIP

function [S,error]=MonteC(L) if nargin==0,L=7;end for k=1:L S(k)=viviani; end error=S-(8/3*pi-32/9); hist(S,7);实验任务一：L次实验的实验数据及误差;[V,er]=montec(128) V = 4.8226 er = 6.0640e-004;function vivian [X,Y,Z]=sphere(36); colormap([0 0 1]) mesh(2*X,2*Y,2*Z),axis off hold on,view(150,24) [x,y,z]=cylinder([1,1],36); x=x+1;z=4.2*z-2; II=31:37; for k=30:-1:1 II=[k,II]; mesh(x(:,II),y(:,II),z(:,II)),pause(.5) end;思考与练习