diff --git a/cpp/switchchain_canonical_mixingtime.cpp b/cpp/switchchain_canonical_mixingtime.cpp
index 77d47dd39c6ae922e1b1a0b9f80c7e50052101a0..31e2ef690fde93856529d961936300272108b687 100644
--- a/cpp/switchchain_canonical_mixingtime.cpp
+++ b/cpp/switchchain_canonical_mixingtime.cpp
@@ -13,9 +13,9 @@
 
 int main(int argc, char* argv[]) {
     // Simulation parameters
-    const int numVerticesMin = 2000;
-    const int numVerticesMax = 5000;
-    const int numVerticesStep = 1000;
+    const int numVerticesMin = 10000;
+    const int numVerticesMax = 20000;
+    const int numVerticesStep = 10000;
 
     float tauValues[] = {2.1f, 2.3f, 2.5f, 2.7f, 2.9f};
 
@@ -52,7 +52,7 @@ int main(int argc, char* argv[]) {
     outfile << "tauValues: " << tauValues << std::endl;
     outfile << "Canonical degree sequence.\n";
     outfile << "sample runs: " << sampleRuns << std::endl;
-    outfile << "time stamps: {0.1 n, 0.2 n, ... , 20.0 n}\n";
+    outfile << "time stamps: {0.2 n, 0.4 n, 0.6 n, ... , 40.0 n}\n";
     outfile << "For uniform samples:\n";
     outfile << "mixingTime: 50 * (50 - 5 (tau - 2)) n\n";
     outfile << "measurements: 100000\n";
@@ -103,7 +103,7 @@ int main(int argc, char* argv[]) {
             outfile << '}';
 #else
             std::vector<std::pair<int, Histogram>> samples;
-            for (int i = 1; i <= 200; i++) {
+            for (int i = 2; i <= 400; i += 2) {
                 samples.push_back({(i * numVertices / 10), Histogram()});
             }