diff --git a/docs/src/harmonics/harmonics.md b/docs/src/harmonics/harmonics.md
index 14280452..fa124fdf 100644
--- a/docs/src/harmonics/harmonics.md
+++ b/docs/src/harmonics/harmonics.md
@@ -30,6 +30,7 @@ sd = EmbeddedDeltaDualComplex2D{Bool, Float64, Point3D}(s);
 subdivide_duals!(sd, Barycenter());
 ```
 
+Compute the Laplacian eigenvectors using [LinearAlgebra.eigen](https://docs.julialang.org/en/v1/stdlib/LinearAlgebra/#LinearAlgebra.eigen). This requires making the sparse Laplacian matrix dense with `collect`. Alternatively, use [Arpack.jl](https://arpack.julialinearalgebra.org/stable/).
 
 ```@example Harmonics
 Δ0 = -Δ(0,sd)