diff --git a/src/sage/quadratic_forms/special_values.py b/src/sage/quadratic_forms/special_values.py
index 3a5dc102cc0..2a3d5097d88 100644
--- a/src/sage/quadratic_forms/special_values.py
+++ b/src/sage/quadratic_forms/special_values.py
@@ -258,11 +258,11 @@ def quadratic_L_function__numerical(n, d, num_terms=1000):
         sage: RR = RealField(100)
         sage: for i in range(5):
         ....:     print("L({}, (-4/.)): {}".format(1+2*i, RR(quadratic_L_function__exact(1+2*i, -4)) - quadratic_L_function__numerical(RR(1+2*i),-4, 10000)))
-        L(1, (-4/.)):  0.000049999999500000024999996962707
-        L(3, (-4/.)):  4.99999970000003...e-13
-        L(5, (-4/.)):  4.99999922759382...e-21
-        L(7, (-4/.)):  ...e-29
-        L(9, (-4/.)):  ...e-29
+        L(1, (-4/.)): 0.000049999999500000024999996962707
+        L(3, (-4/.)): 4.99999970000003...e-13
+        L(5, (-4/.)): 4.99999922759382...e-21
+        L(7, (-4/.)): ...e-29
+        L(9, (-4/.)): ...e-29
 
     This procedure fails for negative special values, as the Dirichlet
     series does not converge here::