diff --git a/e030.py b/e030.py
new file mode 100644
index 0000000..f3fa3a6
--- /dev/null
+++ b/e030.py
@@ -0,0 +1,16 @@
+def power_sums(n):
+    numbers = []
+    i = 9
+    while i <= 10**(n+1):
+        i = i + 1
+        if int(max(str(i)))**n > i:
+            continue
+        s = sum([int(c)**n for c in str(i)])
+        if i == s:
+            numbers.append(i)
+    return numbers
+
+p = power_sums(4)
+print 'power_sums(4)', p, sum(p)
+p = power_sums(5)
+print 'power_sums(5)', p, sum(p)