diff --git a/bfmatrix.cpp b/bfmatrix.cpp
index 64c542fa4c2bf2dcad486ea77a3289fe713b8ec4..d0ef3a0e9a125773cf4918b3b1f1e47287582182 100644
--- a/bfmatrix.cpp
+++ b/bfmatrix.cpp
@@ -5,7 +5,7 @@
 //
 //    Copyright (C) 2007 University of Oxford 
 //
-/*  CCOPYRIGHT  */
+
 #include <iostream>
 #include <iomanip>
 #include <fstream>
@@ -25,9 +25,9 @@ NEWMAT::Matrix BFMatrix::SubMatrix(unsigned int fr, unsigned int lr, unsigned in
 {
   if (fr<1 || fc<1 || lr>Nrows() || lc>Ncols() || fr>lr || fc>lc) throw BFMatrixException("BFMatrix::SubMatrix: index out of range");
   NEWMAT::Matrix omat(lr-fr+1,lc-fc+1);
-  for (unsigned int r=fr; r<=lr; r++) {
-    for (unsigned int c=fc; c<lc; c++) {
-      omat(r,c) = this->Peek(r,c);
+  for (unsigned int r=fr, ri=1; r<=lr; r++, ri++) {
+    for (unsigned int c=fc, ci=1; c<=lc; c++, ci++) {
+      omat(ri,ci) = this->Peek(r,c);
     }
   }
   return(omat);