%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-6-7程序
% 蒙特－卡羅隨機測試法求圓周率
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
clear
N=input('N= ')
%x=rand(1,N);y=rand(1,N);
k=0;
for i=1:N
    if rand(1)^2+rand(1)^2<1 k=k+1;end
end
k/N*4