suanfa.txt

一个用于求圆周率pi的matlab算法。题目背景：取一长面积为1平方米的正方形白纸

n=str2num(get(findobj('tag','n'),'string'));
x=rand(n,1);
y=rand(n,1);
t=0:1/100:2*pi;
a=(1+sin(t))/2;b=(1+cos(t))/2;
plot(a,b);
hold on;
d=x.*x+y.*y;
h=d(find(d<1));
no=size(h);
piz=(4*no(1))/n;
set(handles.pi,'string',num2str(piz));