diff --git a/R/multitaperHelper.R b/R/multitaperHelper.R
index c2b7b93..bb3ad83 100644
--- a/R/multitaperHelper.R
+++ b/R/multitaperHelper.R
@@ -68,11 +68,11 @@
                     var, dt_, ev, evp=(1-ev),
                     tol=.03, maxadaptiveiteration=100) {
 
-    ## this is equation (5.3) and (5.4) form 
+    ## this is equation (5.2) and (5.3) from 
     ##   Thomson, D.J. Spectrum Estimation and Harmonic Analysis,
     ##   Proceedings of the IEEE, 1982.
-    ## note that the weights are squared, they are |d_k(f)^2 from equation
-    ## (5.4)
+    ## note that the weights are squared, they are |d_k(f)|^2 from equation
+    ## (5.3)
 
     out <- .Fortran("mw2wta", as.double(sa),
                     wt=matrix(as.double(0), nfreq, nord),