diff --git a/nalgebra-lapack/src/generalized_eigenvalues.rs b/nalgebra-lapack/src/generalized_eigenvalues.rs
index 5c273e9b..4a2a293f 100644
--- a/nalgebra-lapack/src/generalized_eigenvalues.rs
+++ b/nalgebra-lapack/src/generalized_eigenvalues.rs
@@ -194,6 +194,11 @@ where
     ///
     ///  u(j)**H * A  = lambda(j) * u(j)**H * B .
     ///  where u(j)**H is the conjugate-transpose of u(j).
+    ///
+    /// What is going on below?
+    ///  If the j-th and (j+1)-th eigenvalues form a complex conjugate pair,
+    ///  then v_l(j) = VSL(:,j)+i*VSL(:,j+1) and v_l(j+1) = VSL(:,j)-i*VSL(:,j+1).
+    ///  and then v_r(j) = VSR(:,j)+i*VSR(:,j+1) and v_r(j+1) = VSR(:,j)-i*VSR(:,j+1).
     pub fn eigenvectors(self) -> (OMatrix<Complex<T>, D, D>, OMatrix<Complex<T>, D, D>)
     where
         DefaultAllocator: