diff --git a/docs/src/dev/how_it_works.md b/docs/src/dev/how_it_works.md
index fc923d2d..753602e8 100644
--- a/docs/src/dev/how_it_works.md
+++ b/docs/src/dev/how_it_works.md
@@ -102,24 +102,28 @@ two calls to ``f`` to evaluate the gradient:
                x_3 \\
                x_4
            \end{bmatrix}
-
-\vec{x}_{\epsilon} = \begin{bmatrix}
+```
+1st call:
+```math
+\to \vec{x}_{\epsilon} = \begin{bmatrix}
                         x_1 + \epsilon_1 \\
                         x_2 + \epsilon_2 \\
                         x_3 \\
                         x_4
                      \end{bmatrix}
 \to
-f(\vec{x}_{\epsilon}) = f(\vec{x}) + \frac{\delta f(\vec{x})}{\delta x_1} \epsilon_1 + \frac{\delta f(\vec{x})}{\delta x_2} \epsilon_2
-
-\vec{x}_{\epsilon} = \begin{bmatrix}
+f(\vec{x}_{\epsilon}) = f(\vec{x}) + \frac{\delta f(\vec{x})}{\delta x_1} \epsilon_1 + \frac{\delta f(\vec{x})}{\delta x_2} \epsilon_2,
+```
+2nd call:
+```math
+\to \vec{x}_{\epsilon} = \begin{bmatrix}
                         x_1 \\
                         x_2 \\
                         x_3 + \epsilon_1 \\
                         x_4 + \epsilon_2
                      \end{bmatrix}
 \to
-f(\vec{x}_{\epsilon}) = f(\vec{x}) + \frac{\delta f(\vec{x})}{\delta x_3} \epsilon_1 + \frac{\delta f(\vec{x})}{\delta x_4} \epsilon_2
+f(\vec{x}_{\epsilon}) = f(\vec{x}) + \frac{\delta f(\vec{x})}{\delta x_3} \epsilon_1 + \frac{\delta f(\vec{x})}{\delta x_4} \epsilon_2.
 ```
 
 This seeding process is similar for Jacobians, so we won't rehash it here.