diff --git a/mitiq/lre/inference/multivariate_richardson.py b/mitiq/lre/inference/multivariate_richardson.py
index 7e81a10c5..00e3f5670 100644
--- a/mitiq/lre/inference/multivariate_richardson.py
+++ b/mitiq/lre/inference/multivariate_richardson.py
@@ -23,7 +23,8 @@
 def _full_monomial_basis_term_exponents(
     num_layers: int, degree: int
 ) -> list[tuple[int, ...]]:
-    """Finds exponents of monomial terms required to create the sample matrix.
+    """Finds exponents of monomial terms required to create the sample matrix
+    as defined in Section IIB of :cite:`Russo_2024_LRE`.
 
     $Mj(λ_i, d)$ is the basis of monomial terms for $l$-layers in the input
     circuit up to a specific degree $d$. The linear combination defines our