diff --git a/LabPSA_3/distrib.py b/LabPSA_3/distrib.py
index 8b13789..c3b4841 100644
--- a/LabPSA_3/distrib.py
+++ b/LabPSA_3/distrib.py
@@ -1 +1,15 @@
 
+import random
+
+estimation = (50*200)/8
+print('estimated number of deer:' + str(estimation))
+prob_deer, deer, n = 1/25, 0, 10000
+
+for i in range(n):
+    flag = 0
+    while flag < 50:
+        deer += 1
+        if random.uniform(0, 1) <= prob_deer:
+            flag += 1
+
+print('number of deer by simulation:' + str(deer/n))