diff --git a/doc/examples/plot_swirl.py b/doc/examples/plot_swirl.py
index f05246b0..98226a20 100644
--- a/doc/examples/plot_swirl.py
+++ b/doc/examples/plot_swirl.py
@@ -42,9 +42,9 @@ mapping for the swirl transformation first computes, relative to a center
 
 .. math::
 
-    \theta = \arctan(y/x)
+    \\theta = \\arctan(y/x)
 
-    \rho = \sqrt{(x - x_0)^2 + (y - y_0)^2},
+    \\rho = \sqrt{(x - x_0)^2 + (y - y_0)^2},
 
 and then transforms them according to
 
@@ -56,12 +56,12 @@ and then transforms them according to
 
     s = \mathtt{strength}
 
-    \theta' = \phi + s \, e^{-\rho / r + \theta}
+    \\theta' = \phi + s \, e^{-\\rho / r + \\theta}
 
 where ``strength`` is a parameter for the amount of swirl, ``radius`` indicates
 the swirl extent in pixels, and ``rotation`` adds a rotation angle.  The
 transformation of ``radius`` into :math:`r` is to ensure that the
-transformation decays to :math:`\approx 1/1000^{\mathsf{th}}` within the
+transformation decays to :math:`\\approx 1/1000^{\mathsf{th}}` within the
 specified radius.
 
 """