matlab 计算m=226;s=0;for i=1:ma=randint(1,2,[1,10^9]);if gcd(a

matlab 计算
m=226;s=0;
for i=1:m
a=randint(1,2,[1,10^9]);
if gcd(a(1),a(2))==1;
s=s+1;
end
end
pi=sqrt(6*m/s)
求结果
a=0;b=1;m=226;
H=1;s=0;
for i=1:m
xi=rand();
yi=H*rand();
if yi