diff --git a/examples/multiple-meshes.py b/examples/multiple-meshes.py
new file mode 100644
index 0000000000000000000000000000000000000000..84918e1779a34932ed6269a8884e30e63bb846d0
--- /dev/null
+++ b/examples/multiple-meshes.py
@@ -0,0 +1,26 @@
+from __future__ import division
+
+import numpy as np  # noqa
+
+
+order = 4
+
+
+def main():
+    from meshmode.mesh.generation import (  # noqa
+            make_curve_mesh, starfish)
+    mesh1 = make_curve_mesh(starfish, np.linspace(0, 1, 20), 4)
+
+    from meshmode.mesh.processing import affine_map, merge_disjoint_meshes
+    mesh2 = affine_map(mesh1, b=np.array([2, 3]))
+
+    mesh = merge_disjoint_meshes((mesh1, mesh2))
+
+    from meshmode.mesh.visualization import draw_2d_mesh
+    draw_2d_mesh(mesh, set_bounding_box=True)
+
+    import matplotlib.pyplot as pt
+    pt.show()
+
+if __name__ == "__main__":
+    main()