diff --git a/src/WallGo/equationOfMotion.py b/src/WallGo/equationOfMotion.py
index bc669eda..4e217e8f 100644
--- a/src/WallGo/equationOfMotion.py
+++ b/src/WallGo/equationOfMotion.py
@@ -1562,17 +1562,19 @@ def estimateTanhError(
         r"""
         Estimates the EOM error due to the tanh ansatz. It is estimated by the integral
 
-        .. math:: \sqrt{\Delta[\mathrm{EOM}^2]/|\mathrm{EOM}^2|},
+        .. math::
+            \sqrt{\Delta[\mathrm{EOM}^2]/|\mathrm{EOM}^2|},
 
         with
 
-        .. math:: \\Delta[\\mathrm{EOM}^2]=\\int\\! dz\\, (-\\partial_z^2 \\phi+ 
-            \\partial V_{\\mathrm{eq}}/ \\partial \\phi+ \\partial V_{\\mathrm{out}}/ \\partial \\phi )^2
+        .. math::
+            Delta[\mathrm{EOM}^2]=\int\! dz\, (-\partial_z^2 \phi+ 
+            \partial V_{\mathrm{eq}}/ \partial \phi+ \partial V_{\mathrm{out}}/ \partial \phi )^2
 
         and
 
-        .. math:: |\\mathrm{EOM}^2|=\\int\\! dz\\, [(\\partial_z^2 \\phi)^2+ 
-            (\\partial V_{\\mathrm{eq}}/ \\partial \\phi)^2+ (\\partial V_{\\mathrm{out}}/ \\partial \\phi)^2].
+        .. math:: |\mathrm{EOM}^2|=\int\! dz\, [(\partial_z^2 \phi)^2+ 
+            (\partial V_{\mathrm{eq}}/ \partial \phi)^2+ (\partial V_{\mathrm{out}}/ \partial \phi)^2].
 
         """
         Tminus = hydroResults.temperatureMinus