用MATLAB写个蒙特卡罗求圆周率的值 要有MATLAB画图的代码

好的，以下是MATLAB代码实现：

% 设置投点数量
N = 100000;

% 生成随机点并判断是否落在圆内
x = rand(1, N);
y = rand(1, N);
in_circle = (x.^2 + y.^2 <= 1);

% 计算圆周率值
pi_estimate = 4*sum(in_circle)/N;

% 绘制结果
theta = linspace(0, 2*pi, 100);
x_circle = cos(theta);
y_circle = sin(theta);
plot(x(~in_circle), y(~in_circle), '.', 'MarkerSize', 5, 'Color', [0.5, 0.5, 0.5])
hold on
plot(x(in_circle), y(in_circle), '.', 'MarkerSize', 5, 'Color', [0.8500, 0.3250, 0.0980])
hold on
plot(x_circle, y_circle, 'r')
axis equal
title(['估算值为: ', num2str(pi_estimate)])

运行以上代码，即可得到蒙特卡罗估算的圆周率值，并且用图像呈现结果。