diff --git a/src/beyn.py b/src/beyn.py
index 49c9973721bcadce5ee73882cced5d84ce8bacd9..5b6a9c0c118aa02040d5b88739c64004e8dc5634 100644
--- a/src/beyn.py
+++ b/src/beyn.py
@@ -18,7 +18,7 @@ def simple(operator, origin, radius,
     verbose -- should I be verbose? (optional)
 
     Returns the computed eigenvalues, eigenvectors
-    and the associate *absolute* residual norm
+    and the associate *absolute* residual norms
     """
 
     # Display the parameter used
@@ -72,7 +72,7 @@ def simple(operator, origin, radius,
     A1 = integrate(operator, myPath, 1, vHat)
     B  = V0.H * A1 * W0 * S0Inv
 
-    # Eigenvalues & eigenvectors (with projection onto V0)
+    # Eigenvalues & eigenvectors (by projecting QHat onto V0)
     if(verbose): print "Solving linear EVP..."
     myLambda, QHat = numpy.linalg.eig(B)
     Q = V0 * QHat;