diff --git a/cube.sage b/cube.sage
new file mode 100644
index 0000000000000000000000000000000000000000..2a00e844260157314d362dd0069411c01efef814
--- /dev/null
+++ b/cube.sage
@@ -0,0 +1,25 @@
+# generate a nice picture of subdivided cube
+# which is Z_4 but not Z_2^2 connected
+
+from groupConnectivityNaive import *
+
+G = prettifyGraph(Graph('J?CSZ`GHC@_'))
+
+G.set_pos({ 
+  7: [50,50], 1: [100,50], 8: [150,50], 10: [150,150], 9: [50,150], 2: [50,100],
+  3:[0,0], 4:[200,0], 6: [200,200], 5:[0, 200],
+  0: [175,175]
+})
+
+l = testGroupConnectivityNaive(G, Zkn(2,2), list_all = False)
+
+def labeledGraph(G_, labels):
+  G = G_.copy()
+  for (u, v, i) in G.edges():
+    G.set_edge_label(u, v, labels[i])
+  return G
+
+P = labeledGraph(G, l[2][0]).plot(edge_labels = True)
+
+P.save('cube.pdf')
+