matlab蒙特卡洛求体积

 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