diff --git a/verzweigungsvectoren/main.py b/verzweigungsvectoren/main.py
new file mode 100755
index 0000000..4c7a378
--- /dev/null
+++ b/verzweigungsvectoren/main.py
@@ -0,0 +1,31 @@
+#!/usr/bin/env python
+
+from numpy import roots,isreal
+from itertools import combinations_with_replacement
+
+
+solutions = []
+
+s = [1, 2, 3, 4, 5, 6, 7, 8, 9]
+
+for i in range(2, 7):
+
+    for j in combinations_with_replacement(s, i):
+
+        p = [-j.count(k) for k in range(max(j) + 1)]
+        p[0] = 1
+
+
+        r = [float(k) for k in roots(p) if isreal(k) and k > 0.0]
+        # print(j, max(j), p, r)
+
+        solutions.append((j, r[0]))
+
+
+
+solutions.sort(key=lambda e: e[1])
+
+for s in solutions:
+
+    print(f"{s[1]}\t{s[0]}")
+