主题：【讨论】给大家出一道初中数学题 -- 迷途笨狼

n=10000;

l=20;

num=0;

for i=1:n

x=rand(1,2);

if x(2)<sqrt(1-x(1)^2)

if x(2)>(sqrt(4-x(1)^2)-sqrt(2))

num=num+1;

end

end

end

S=4*l^2*num/n;

S算出来就是阴影面积

数值积分或者蒙特卡洛积分比较靠谱

坐标原点画圆形、正方形，以及对应的“枣核型”

x两个数字，对应了一个二维面上的坐标，生成一万次，坐标将会落在原点中心、变长2的盒子里面

第一个条件是圆形之内，第二个条件要我写就是

(x-(-1))^2 +(y-(-1))^2>2^2

即

y>sqrt(-x^2-2*x+3)-1

表示了右上角那个“西瓜牙”，是在左下角那个扇形的外面

这个西瓜牙的面积，除以整个正方形的面积，就等于被记录的num数量除以总体的n

还是你的好

好久没有手算，一个矩阵不满秩已经没有直觉了

应该在最外面再加一个循环，按照百分之几的收敛率判断一下。

这个问题是有积分的解析解的，俺算了好几遍，极坐标，直角坐标都试了，总是算错，唉，老了。

不知道这个误差能不能有显式的分析

def con1(x,y):

——return y<np.sqrt(1-x**2)

def con3(x,y):

——return y>(np.sqrt(4-x**2)-np.sqrt(2))

num=0

total = 10**6

for i in range(total):

——x = 2*np.random.rand(1)-1

——y = 2*np.random.rand(1)-1

——if con1(x,y) and con2(x,y):

———— num+=1

S = 2*num/total*(40**2)

print(S,num)

具体俺不怎么记得了，什么大数定律，中心极限定律什么的。

简单粗暴的方法就是设置收敛率，把收敛率设置小一点就好了。

快过去看看，记得不要躲，乖乖回答问题

mspaint还能干这个！！

少则得，多则惑，老子诚不我欺。

最近越来越迷惑了。

东施效颦啊，估计学这个比学解题有意思。