% ch1example3prg2.m
tic                             % 計(jì)時(shí)器啟動(dòng)
n=10000;                        % 每次隨機(jī)落點(diǎn)10000個(gè)
for k=1:1000                    % 重復(fù)試驗(yàn)1000次
    x1=2*rand(n,1)-1;
    y1=2*rand(n,1)-1;           % 隨機(jī)落點(diǎn)產(chǎn)生
    m(k)=sum((x1.*x1+y1.*y1)<1);% 求落入圓中的點(diǎn)數(shù)和
end
S_c=mean(m).*4./n               % 計(jì)算并顯示結(jié)果
time=toc                        % 顯示耗時(shí)