1 files changed, 6 insertions, 0 deletions
diff --git a/README.md b/README.md
index 53ef9ce..cb1f485 100644
--- a/README.md
+++ b/README.md
@@ -56,3 +56,9 @@ axis = axis / np.linalg.norm(axis)
 max_planar_angle = np.pi/3
 print(random_on_cap(axis, max_planar_angle, rng))
 ```
+
+# Citing
+
+If you use this code for your research, please cite the paper [An O(n)
+algorithm for generating uniform random vectors in n-dimensional
+cones](https://arxiv.org/abs/2101.00936).