ds-graph-based-phasespace-exercise.py

from pylab import *
import networkx as nx

g = nx.DiGraph()

for x in range(100):
    g.add_edge(x, x ** x % 100)

ccs = list(nx.connected_components(g.to_undirected()))
n = len(ccs)
w = ceil(sqrt(n))
h = ceil(n / w)
for i in range(n):
    subplot(h, w, i + 1)
    nx.draw(nx.subgraph(g, ccs[i]), with_labels = True)

show()