--- a/src/agf_io.cc
+++ b/src/agf_io.cc
@@ -175,7 +175,14 @@
     mat=read_matrix<real_a, nel_ta>(fs, m, n);
     if (mat==NULL || m<=0 || n<=0) goto fail;
     ave=new real_a[m];
-    for (dim_ta i=0; i<m; i++) ave[i]=0;
+    //printf("read_stats2: constant term:\n");
+    n--;
+    for (dim_ta i=0; i<m; i++) {
+      ave[i]=mat[i][n];
+      printf("%g\n", ave[i]);
+    }
+    //printf("read_stats2: transformation matrix:\n");
+    //print_matrix(stdout, mat, m, n);
   }
 
   delete [] header;