clc;
%設圓的半徑r=1；對右上角1/4區域投點
N=50000;?????
for?n=100:N?????%從100開始，計算不同飛鏢數下PI值
????k=0;????????????%落入圓中的飛鏢數，初始為0
????for?i=1:n???????%將n個飛鏢依次投一遍
?????????x=rand（）;???%隨機產生取值在[01]之間的x坐標值
?????????y=rand（）;???%隨機產生取值在[01]之間的y坐標值
????????if?x*x+y*y<1?%統計在圓內的飛鏢數
?????????????k=k+1;?????
????????end
????end
?PI（n）=4*k/n;???%計算PI值
end
t=100:100:N;????%每隔100個點，繪出Pi值
plot（tPI（t）t3.1415926‘r-‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????392??2016-11-12?20:05??pi_value.m

?????文件????????504??2016-11-12?19:04??pi1.m

?????文件????????125??2016-11-12?20:21??pi4.m

?????文件????????447??2016-11-12?20:09??second_pi_value.m

-----------?---------??----------?-----??----

?????????????????1468????????????????????4