diff --git a/ABACUS_Usage_Example_LiebLin.cc b/ABACUS_Usage_Example_LiebLin.cc
index fc8447e..377ad2f 100644
--- a/ABACUS_Usage_Example_LiebLin.cc
+++ b/ABACUS_Usage_Example_LiebLin.cc
@@ -22,9 +22,9 @@ int main()
 {
   clock_t StartTime = clock();
 
-  DP c_int = 2.0;
-  DP L = 3.0;
-  int N = 3;
+  DP c_int = 25.0;
+  DP L = 50.0;
+  int N = 10;
   DP nbar_required = 1.0;
   DP kBT = 4.0;
   int Npts = 4*N;
@@ -34,9 +34,9 @@ int main()
 
   if (true) { // State-by-state checks
 
-    DP c_int = 4.0;
-    DP L = 16.0;
-    int N = 16;
+    DP c_int = 25.0;
+    DP L = 50.0;
+    int N = 10;
 
     LiebLin_Bethe_State gstate (c_int, L, N);
     gstate.Compute_All(true);
@@ -71,9 +71,9 @@ int main()
     // 	 << exp(real(ln_Density_ME(gstate, estate))) << "\t" << exp(real(ln_Density_ME(estate, gstate))) << endl;
 
     //Field operator matrix elements:
-    cout << "omega\tiK\t< estate | Psi | gstate > matrix element:" << endl;
-    cout << setprecision(16) << estate.E - gstate.E << "\t" << estate.iK - gstate.iK << "\t"
-	 << exp(ln_Psi_ME(estate, gstate)) << endl;
+    // cout << "omega\tiK\t< estate | Psi | gstate > matrix element:" << endl;
+    // cout << setprecision(16) << estate.E - gstate.E << "\t" << estate.iK - gstate.iK << "\t"
+    // 	 << exp(ln_Psi_ME(estate, gstate)) << endl;
 
     //LiebLin_Bethe_State flipstate = estate;
     //flipstate.Parity_Flip();
diff --git a/Makefile b/Makefile
index 76b84fc..76fe22d 100644
--- a/Makefile
+++ b/Makefile
@@ -24,7 +24,8 @@ BINDIR = $(BASEDIR)bin/
 
 # Set the compiler choice
 #CXX = g++
-CXX = /usr/local/bin/g++
+#CXX = /usr/local/bin/g++
+CXX = /usr/local/Cellar/gcc@7/7.4.0/bin/g++-7
 # On OS X, using a separate installation of llvm tools:
 # CXX = /usr/local/opt/llvm/bin/clang++ -Xpreprocessor -fopenmp -lomp -L/usr/local/opt/llvm/lib:lib/ -I/usr/local/opt/llvm/include
 # CXX = clang++ -Weverything -ferror-limit=1