blob: f8167cc6ab31d5b42de62ebcda57bbe22aaf1700 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
import random
import matplotlib.pyplot as plt
denominator = 0
numerator = 0
results = []
for N in range(10000, 1000000, 10000):
for i in range(N):
b = random.random() # 0-1
a = random.random() # 0-1
if a < b:
denominator += 1
if a+(b-a) > 1-b and a+(1-b) > b-a and (b-a)+(1-b) > a:
numerator += 1
print(N, numerator*1./denominator)
results.append(numerator*1./denominator)
# print(results)
plt.scatter(range(len(results)), results)
plt.show()
|
