[Alg]利用蒙特卡洛思想计算圆锥体积

%% 用蒙特卡洛思想算圆锥体积
clear
clc

%% 
R = 0.5;
H = 0.1;
n = 10000000;%R和H越小，模拟次数越大，精确度越好

Pr = unifrnd(-R,R,2,n);
Ph = unifrnd(0,H,1,n);
cnt = 0;

for i = 1:n
    r = sqrt(Pr(1,i)^2 + Pr(2,i)^2);
    h = H-Ph(i);%坐标系不一样，也可以 h=Ph(i)   
    if (r/R)<(h/H)
        cnt = cnt+1;%说明点在圆锥里面
    end
end

V_f = 1/3*pi*R^2*H;%公式计算圆锥体积
V_rect = (2*R)^2*H;%公式计算长方体体积
scale = cnt/n;%圆锥体积占长方体的比例
V_m = scale*V_rect%利用蒙特卡洛思想计算圆锥体积
error_cmp = abs(V_m - V_f)%蒙特卡洛思想算的体积与公式计算的体积的误差