answered 8 years ago
You can do the disjoint union and then identify some vertices:
sage: I = g.disjoint_union(h) sage: I.merge_vertices([(0, 0),(1, 0)]) sage: I.plot()