From 9aeaa9dc06b5ecd0d3c5f98b0bf80e8af0556e4b Mon Sep 17 00:00:00 2001 From: zhang Date: Sat, 9 Nov 2019 16:30:58 +0800 Subject: first commit --- monte_carlo_triangle.py | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100644 monte_carlo_triangle.py (limited to 'monte_carlo_triangle.py') diff --git a/monte_carlo_triangle.py b/monte_carlo_triangle.py new file mode 100644 index 0000000..f8167cc --- /dev/null +++ b/monte_carlo_triangle.py @@ -0,0 +1,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() -- cgit v1.2.3