以下是一个简单的蒙特卡罗模拟的Matlab代码示例:
num_samples = 10000;
in_circle = 0;
for i = 1:num_samples
x = rand;
y = rand;
if x^2 + y^2 <= 1
in_circle = in_circle + 1;
end
end
pi_estimate = 4 * in_circle / num_samples;
disp(pi_estimate);
在这个例子中,我们生成了10000个随机数对(x, y),判断它们是否在单位圆内,然后计算在