MATLAB作业题，蒙特卡洛模拟

img

你好，考试题，请自己回答，这样更能学习进步哦


subplot(1,2,1);
% x 0~1
x=0:0.01:1;
y = sqrt(1-x.^2);
plot(x,y);
hold on;
%【输入n】
n=100;
A=unifrnd(0,1,n,2);

plot(A(:,1),A(:,2),'.');

InArae = A(A(:,2)<=sqrt(1-A(:,1).^2),:);
Pinsum = size(InArae,1);
%总点数
Psum = n;

%估计面积
Sc = Sall* Psum/Pinsum;
title('n=100');
xlabel(['估计面积为',num2str(Sc)]);



subplot(1,2,2);
% x 0~1
x=0:0.01:1;
y = sqrt(1-x.^2);
plot(x,y);
hold on;
%【输入n】
n=10000;
A=unifrnd(0,1,n,2);

plot(A(:,1),A(:,2),'.');

InArae = A(A(:,2)<=sqrt(1-A(:,1).^2),:);
Pinsum = size(InArae,1);
%总点数
Psum = n;

%估计面积
Sc = Sall* Psum/Pinsum;
title('n=10000');
xlabel(['估计面积为',num2str(Sc)]);

img