diff --git a/src/rmn.cpp b/src/rmn.cpp
index 5f9649f..92a499a 100644
--- a/src/rmn.cpp
+++ b/src/rmn.cpp
@@ -57,6 +57,7 @@ void estimateZ(
   for (i = 1; i < L + 1; ++i, ++F_it, ++B_it){
 
     if (i >= w){
+      std::cout << i << '\t' << *forlnEmPN_it << std::endl;
       *F_it = logSumExp(0, *forlnEmPN_it++ + lnTPN + *(F_it - w) - *(F_it - 1));
       *B_it = logSumExp(0, *revlnEmPN_it++ + lnTPN + *(B_it - w) - *(B_it - 1));
 
@@ -69,6 +70,8 @@ void estimateZ(
     *F_it += *(F_it - 1);
     *B_it += *(B_it - 1);
   }
+
+  std::cout << "lnL = " << F[L] << std::endl;
   lnL += F[L];
 
   // calculate posterior;