matlab蒙特卡洛法求pi动图

好的，这里是matlab蒙特卡洛法求pi的动图代码：

N = 10000; % 投点次数
x = rand(1,N); % 在[0,1]区间内随机生成x坐标
y = rand(1,N); % 在[0,1]区间内随机生成y坐标
r = sqrt(x.^2 + y.^2); % 计算每个点到原点的距离
count = sum(r <= 1); % 统计落在圆内的点数
pi_estimate = 4 * count / N; % 估算pi值
for i = 1:N
    if r(i) <= 1 % 落在圆内的点用红色表示
        plot(x(i),y(i),'r.');
    else % 落在圆外的点用蓝色表示
        plot(x(i),y(i),'b.');
    end
    axis equal; % 设置坐标轴比例为1:1
    title(['\pi = ' num2str(pi_estimate)]); % 显示当前估算的pi值
    drawnow; % 实时更新图像
end