first commit

1 files changed, 27 insertions, 0 deletions

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()