diff --git a/python/cusps/maximal_cusp_area_matrix.py b/python/cusps/maximal_cusp_area_matrix.py
index d463fe39b..12ca431ae 100644
--- a/python/cusps/maximal_cusp_area_matrix.py
+++ b/python/cusps/maximal_cusp_area_matrix.py
@@ -23,6 +23,9 @@ def maximal_cusp_area_matrix(manifold, bits_prec, verified):
     [14.3590180492058 88.4588035788533 9.67661682098102 11.4136568679705]
     [11.4136568679715 9.67661682098102 88.4588035788541 14.3590180492042]
     [9.67661682098105 11.4136568679705 14.3590180492042 88.4588035788038]
+
+    sage: from snappy import Manifold
+    sage: M = Manifold("o9_44206")
     sage: maximal_cusp_area_matrix(M, bits_prec=80, verified=True) # doctest: +NUMERIC6
     [ 88.458803578854197094? 14.3590180492058335371? 11.4136568679715291317?  9.6766168209810445566?]
     [14.3590180492058335371?          88.4588035789?      9.676616820981044?          11.4136568680?]