本文共 1210 字，大约阅读时间需要 4 分钟。

MATLAB���������

P=rand(100000,2);%���������������������������������x���������������y���������100000������
x=2*P(:,1)-1;%x������������-1���1
y=2*P(:,2);%y������������0���2
II=find(y<=2-x.^2 & y.^3>=x.^2);%������������������������
M=length(II);%������������������������������
S=4*M/100000;%������������������
plot(x(II),y(II),'g');%������������������������

������������

���������������������������100000���������������������������y ��� 2 - x�����y�� ��� x�����������������������������������������������������M���������������������������������������������������������������������������������

#������������������

�����������������������������������������������������������������������������������������������������������������

N=100000;%������������������
x=rand(N,1);%������x������������������-1���1
y=rand(N,1);%������y������������������-1���1
count=0;%������������������
for i=1:N
    if (x(i)^2+y(i)^2 <= 1)
        count=count+1;%%������������������������������
    end
end
PI=4*count/N;%������PI���������

�����������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������