diff --git a/src/main.cpp b/src/main.cpp
index c5f1ab1..b6760cf 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -139,9 +139,8 @@ void simulate_100_particles()
 
     double time = 50.; // microseconds
 
-    //trap.write_simulation_to_dir("output/simulate_100_particles", time, N,
-                                 //"rk4", false);
-    trap.simulate(time, N);
+    trap.write_simulation_to_dir("output/simulate_100_particles", time, N,
+                                 "rk4", false);
 }
 
 /** @brief Simulate 100 particles over 500 \f$ \mu s \f$ using a time
@@ -264,16 +263,16 @@ int main()
 {
     double t0 = omp_get_wtime();
 
-    //simulate_single_particle();
+    simulate_single_particle();
 
-    //simulate_two_particles();
+    simulate_two_particles();
 
-    //simulate_single_particle_with_different_steps();
+    simulate_single_particle_with_different_steps();
 
     simulate_100_particles();
 
-    //simulate_100_particles_with_time_potential_wide_sweep();
-    //simulate_100_particles_with_time_potential_narrow_sweep();
+    simulate_100_particles_with_time_potential_wide_sweep();
+    simulate_100_particles_with_time_potential_narrow_sweep();
 
     double end = omp_get_wtime();