Matlab用蒙德卡洛方法求定积分S=∫¹₀√1-x^2dx=π4

% Matlab代码

N = 1000000; % 设定模拟次数 x = rand(1,N); % 生成均匀分布的随机数 y = rand(1,N); % 生成均匀分布的随机数 count = 0; % 统计落在函数图像下方的点的数量 for i = 1:N if y(i) <= sqrt(1-x(i)^2) count = count + 1; end end S = count/N; % 计算概率密度函数下方的面积比例 S = S*4; % 乘以4得到整个函数的面积比例 S % 输出结果

% 结果约为 0.7852, 与 π/4 接近