Add example

src/sage/geometry/polyhedron/base.py

@@ -10441,6 +10441,16 @@ def affine_hull_manifold(self, name=None, latex_name=None, start_index=0, ambien
             E^3 --> A
                (x, y, z) |--> (x0, x1) = (-x + y + 1, -1/2*x - 1/2*y + z + 1/2)
 
+        Arrangement of affine hull of facets::
+
+            sage: D = polytopes.dodecahedron()
+            sage: E3 = EuclideanSpace(3)
+            sage: submanifolds = [
+            ....:     F.as_polyhedron().affine_hull_manifold(name=f'F{i}', orthogonal=True, ambient_space=E3)
+            ....:     for i, F in enumerate(D.facets())]
+            sage: sum(FM.plot({}, srange(-2, 2, 0.1), srange(-2, 2, 0.1), opacity=0.2)
+            ....:     for FM in submanifolds)  # not tested
+
         """
         if ambient_space is None:
             from sage.manifolds.differentiable.examples.euclidean import EuclideanSpace