diff --git a/README.md b/README.md
index afdc88f9e..096a0c216 100644
--- a/README.md
+++ b/README.md
@@ -35,3 +35,16 @@ u0 = (1.0, 2.0) # brackets
 prob = IntervalNonlinearProblem(f, u0)
 sol = solve(prob)
 ```
+
+## Citation
+
+If you found this library to be useful in academic work, then please cite:
+
+```bibtex
+@article{pal2024nonlinearsolve,
+  title={NonlinearSolve. jl: High-Performance and Robust Solvers for Systems of Nonlinear Equations in Julia},
+  author={Pal, Avik and Holtorf, Flemming and Larsson, Axel and Loman, Torkel and Schaefer, Frank and Qu, Qingyu and Edelman, Alan and Rackauckas, Chris and others},
+  journal={arXiv preprint arXiv:2403.16341},
+  year={2024}
+}
+```