diff --git a/plot.gnuplot b/plot.gnuplot
new file mode 100755
index 0000000..1b3c8d7
--- /dev/null
+++ b/plot.gnuplot
@@ -0,0 +1,7 @@
+#!/usr/bin/env gnuplot
+set terminal pdf enhanced
+
+set output "mutualinformation.pdf"
+set xlabel "t"
+set ylabel "Mutual information"
+plot "mutualinformation" using 1:2 title "" with lines lw 3, 0 title ""
diff --git a/scrambling.cpp b/scrambling.cpp
index e8a62fd..4da688c 100644
--- a/scrambling.cpp
+++ b/scrambling.cpp
@@ -31,16 +31,16 @@ int main()
 	double tPlus = 0;  // time on the left boundary
 	double tMinus = 0; // time on the right boundary
 	alpha = 0.4; // encodes the conformal dimention of local operator h_Psi
-	double y = 3; // endpoint of interval A
-	double L = 15; // length of interval A
+	double y = 1; // endpoint of interval A
+	double L = 5; // length of interval A
 	double epsilon = 0.01; // smearing parameter
 	beta = 10; // inverse temperature
-// 	double tOmega = /*5.4418*/; // thermal state is perturbed by operator inserted at time -tOmega
+	double tOmega = 4; // thermal state is perturbed by operator inserted at time -tOmega
 	double c = 600; // central charge. Must be large in our approximation
 	// End of parameters
 	for(unsigned int i = 0; i < 1000; ++i)
 	{
-		double tOmega = i*2*L/1000;
+		tOmega = i*2*L/1000;
 
 		// Operator insertion points: Left boundary
 		std::complex<double> x1 (0, -epsilon), x4 (0, epsilon), x1bar, x4bar;
@@ -79,7 +79,8 @@ int main()
 		double S_thermal = 2*c/3 * log(sinh(M_PI*L/beta)/cosh(M_PI/beta*(tMinus-tPlus)));
 		double I = S_A + S_B - S_union + S_thermal;
 		std::cout << tOmega << "\t" << I << std::endl;
-}
+	}
+
 	return 0;
 }