diff --git a/poisson/main.py b/poisson/main.py
new file mode 100644
index 0000000..dda51d7
--- /dev/null
+++ b/poisson/main.py
@@ -0,0 +1,11 @@
+import math
+
+a = list(map(float, input().split(' ')))
+rate = a[0]
+
+sums = 0
+
+for i in range(1,len(a)):
+    sums += round(math.e**(-1*rate) * (rate**a[i]/math.factorial(int(a[i]))), 4)
+
+print(sums)