function[vm]=mont_v(n)
%求体积的蒙特卡洛模拟值
for i=1:length(n)
x=rand(n(i),1); %点的x坐标
y=rand(n(i),1); %点的y坐标
z=rand(n(i),1);%点的z坐标
%落到区域T内的点的频数
m=sum((x.^2+y.^2+z.^2<=1)&(x>0)&(y>0)&(z>0));
vm(i)=m/n(i); %落到所求立体内的点的频率
end
matlab蒙特卡洛求体积
猜你喜欢
转载自blog.csdn.net/qq_24624539/article/details/81462063
今日推荐
周排行